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Abstract 

Present data on the partonic content of the photon from LEP, TRISTAN 
and HERA accelerators are reviewed and the essential aspects of the underlying 
ideas and methods are pointed out. Results of the unpolarized photon structure 
function F2 from DlSg^ experiments and on large px jet production processes in 
the resolved 77 collisions are presented for both real and virtual photons. The 
results of analysis of the hadronic final state accompanying the DlSe^ measure- 
ments, showing some discrepancies with the Monte Carlo models, are collected 
together and presented as a separate issue. Also results on the DISe7 with 
leptonic final states are shown. The results from resolved real and virtual pho- 
ton processes at HERA collider based on the single and double jet events, also 
charged particles and prompt photons, are presented. In the context of virtual 
photon processes the data for forward jet and forward particle production are 
included. In addition a short presentation of the recent data on the heavy quark 
content of the photon is given. Related topics - the polarized (spin dependent) 
structure functions for the real and virtual photon, the structure function of the 
electron and the photonic content of the proton are also shortly mentioned. 

PACS numbers: 13.60. Fz, 12.38. Bx, 14.70. Bh 
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1 Introduction 



1.1 The concept of the "structure of the photon" 

The photon is one of the basic elementary particles. This "horrendous-looking hypoth- 
esis of light-quanta" introduced in 1905 by Einstein has been based on the Planck's 
concept of quanta of energy dating from 1900. The photon properties, as follows from 
Quantum Electrodynamics (QED), are well known: the photon is a massless, charge- 
less object with a pointlike coupling to the charged fundamental particles. As such, it 
was used in the past and is still used at present as an ideal tool for probing structure 
of more complicated objects, for example hadrons, acting as a microscope with the 
resolution given by its wavelength. 

The concept of "the structure of the light-quanta" which has appeared later and 
which is the topic of this work, makes people confused even today. However it is very 
natural in the quantum field theory, since a photon, as any elementary particle, can 
fluctuate into various states consisting of leptons, quarks, bosons, hadrons.... A 
transition to a particular observed state can be realized as a physical process due to 
the interaction with another object (s), for example, ... "through an interaction with a 
Coulomb field the photon could materialize as a pair of electrons, 7 — >■ e^e~ . Although 
not usually thought of in these terms, this phenomenon was the earliest manifestation 
of photon structure" 

Similarly the process ^ qq leads to the "hadronic (partonic) structure of the 
photon" in the perturbative Quantum Chromodynamics (pQCD). This notion offers 
an effective way to describe hard and semihard processes involving photons. The 
hadronic properties of the photon become apparent also in the soft photon-hadron 
interactions, where the similarity between photon and vector mesons p,uj,(j) interaction 
is observed. Probing the "structure" of the photon can be performed in an analogous 
way as testing the structure of other objects - once a hard probe is at our disposal, 
the partonic structure can be revealed. Often the probe happens to be just another 
photon (..."an ideal tool"). Sometimes this leads to additional complications and even 
confusion when the partonic content of the photon is discussed. 

Contrary to the nucleon case, the Quark-Parton Model (QPM or PM) gives the 
definite predictions for the hadronic structure functions F12 for the real photon, and 
for the related parton densities in the photon. They exhibit a logarithmic dependence 
on the scale. The leading QCD corrections to these quantities can be summed 
up using the evolution equations, which, due to existence of the inhomogeneous (PM) 
term, can be solved in principle without any input, leading to the so called asymptotic 
solution. Despite this fact evolution equations are solved usually in an analogous way 
as for the nucleon, i.e. by introducing some experimental or model input, to cure the 
singular behaviour present in the asymptotic solution. The gross features of the parton 
densities in the photon present already in the PM, as the scaling violation and large 
quark densities at large xsj, are unique for the photon (in general for gauge bosons), 
therefore they are treated as one of the basic tests of QCD. This holds to even larger 
extent for the structure functions of the virtual photon. Here the existence of the 
second scale, the "mass" of the virtual photon, makes it possible, at least in principle, 
to obtain the QCD predictions in the input-independent way. The QCD description 
of all basic processes induced by the real and virtual photons is now available with the 
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NLO accuracy. 

Basic ideas related to the considered topic are discussed in review articles, e.g. 



1.2 Recent measurements 

There are two basic types of inclusive processes where the structure of photon is tested 
in existing experiments: 

• the deep inelastic scattering, DlSg^ and DlSg-y*, where the structure functions 
for a real and virtual photon are measured, 

• the large pt jet (and particles) production in 7p and 77 collisions, where indi- 
vidual quark and gluon densities in the photon may be probed (the so called 
resolved photon processes). 

The early experiments of the first type were performed at PETRA and PEP e^e~ 
colliders almost twenty years ago, starting with the first measurement by the PLUTO 
Collaboration in 1981. Recent results on the structure function F2, for a real 7, 
based on few years runs at LEPl, i.e. at the CM energy ~ M^, are now available 
in published or preliminary form. Also new data taken at LEP 1.5 (CM energy 130- 
136 GeV) and at LEP2 (161-172, 183, 189 and 202 GeV) have been or are being 
analysed. Final analyses appeared from the TRISTAN collider (CM energy ~ 60 
GeV). [] Altogether the measurements cover a wide range of the (average) from 
0.2 to ~ 390 — 400 GeV^ allowing to test the scaling violation pattern. The range 
of the XBj variable extends from ~ 0.001, where the possible rise of the structure 
function of the real photon can be observed, to 0.98. In earlier analyses the charm 
(sometimes also bottom) quark contribution to F2 was subtracted, later on it was 
included, although with a different treatment than for the light quarks (only the QPM 
term). Recently the first dedicated measurement of F2(. has been performed at LEP2 
for 0.0014< XBj <0.87 and 5< <ioo GeV^. 

The measurement of the resolved real photon processes, i.e. production of large 
Pt jets (also particles), started a decade later. First evidence for the resolved photon 



process was reported in 1992 by the AMY collaboration at TRISTAN Data have 
been and are still being taken in 77 collisions at the abovementioned e+e~ machines 
and in the photoproduction processes at the e^p collider HERA (with the CM energy 
~ 300 - 320 GeV). The latter one gave first results on the effective parton density (the 
light quark distributions combined with the gluon distributions) and on the gluonic 
content of the photon, showing the rise at small x, a part of the momentum of the 
photon taken by a parton. The probed range of x is from 0.05-0.1 to 1, where x 
equal to 1 corresponds to the direct interaction of the photon. The hard scale in such 
processes, Q, is usually given by pr, and the range of the measured p^ (= Q"^) extends 
from ~ 20 up to 1000 GeV^. There appeared also first data on the prompt photon 
production, and on the c and h quark content of the photon, both from e^e~ and e^p 
collisions. 



^In principle the measurement could be performed also at the SLC collider (at energy ~ 91 GeV), 
but the luminosity is too low for a precise analysis of two-photon processes. 
•^The first high results were obtained by the AMY Collaboration |^ 
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An impressive progress has been observed during the last year in the measurements 
of the "structure of the virtual photon" . Fifteen years after the first measurement of 
the F2 by PLUTO Collaboration ^j, a new measurement has just been performed at 
LEP. A large amount of data concerning the resolved virtual photon comes from the 
ep collider HERA, with the virtualities of the photon spanning the range from 0.1 to 
85 GeV^. Here an effective parton distribution for a virtual photon has been extracted 
from the dijet production. Here also a new type of analysis related to the forward jet 
production in events with large (DISep-type events) is being performed, with the 
aim to test the underlying dynamics of the QCD radiation processes at the small x bj 
in the proton. The resolved virtual photon subprocesses contribute here as well and 
they should be included when describing the forward jet production. 

The discrepancies between the data and Monte Carlo (MC) models in the hadronic 
final states in the photon-induced processes are observed by many collaborations. They 
are especially visible in the DISe-y experiments, where there is a need of unfolding, i. e. 
reconstruction of basic kinematical variable xbj from the visible energy of hadronic 
system W. The improvement in the unfolding has been obtained recently, still the 
dependence on the chosen MC used in the analysis cannot be avoided [0. Similar 
situation, i.e. the lack of proper description of hadronic final state, occurs also at 
the ep collider HERA. It was found there that implementation of the modified trans- 
verse momenta distribution of the partons in the MC programs helps to describe the 
data. The same modified MC models have been used recently in analysis of the e~^e~ 
data from the LEP and TRISTAN colliders. The observed discrepancies may also be 
a signal of a new production mechanism not included in the standard perturbative 
QCD analysis, as for example the multiple interaction. This mechanism, when added, 
usually improves description of the data, both at the e~^e~ and ep colliders. The situ- 
ation is still not clear and it enforces the advanced study of various aspects of hadron 
production in the DlSe^ setup (see e.g. p, |^, |TT|, ^) and in the photon- induced jet 
production [ff^ . 



1.3 The aim of the survey 

The aim of this survey is to present our knowledge on the structure of the photon, 
as it stands AD 2000. We chose to use the data as a guide through the field. We 
collected main results from precision measurements for the unpolarized real and virtual 
photon^, performed during last few years at LEP, TRISTAN and HERA collider^. 
In order to show the status and also the progress achieved in this rapidly expanding 
field of research in high energy physics, we present the very recent (mainly published) 
results, and also some of the previous measurements playing the role of the milestones 
in the field.Q 

The survey is a collection of short presentations of the various experimental data 
without much detail on the experimental setup. Only basic information on the type 
of measurements for each entry, used Monte Carlo generators and, wherever possible, 
comparison with the theoretical predictions are given. This allows to compare different 

■^So far there is no data for the polarized photon structure functions nor for the individual polarized 
parton densities. 

^ We refer to data taken up to 1998, and papers that appeared till June 2000. 

comprehensive collection of the earlier data can be found in jl^; for the new results see |2)h, j. 
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results in an easy way, and also to show various attempts to improve the Monte Carlo 
description of the data. Although the collection of the experimental results is the main 
body of the survey, we found it important to point out the essential aspects of the 
analyses and to mention also some new theoretical ideas and future measurements. 

This survey is an extension of our first attempt p|i, which covered the data which 
appeared up to early 1998 when the first measurements of some quantities related to 
the structure of photon have appeared, and also some problems in describing data 
have been recognized. The present overview contains much more data, and it reports 
new types of analyses undertaken to solve the observed problems ( already with some 
positive results). 

1.4 The content of the survey 

We start with the results of measurements of the hadronic structure function of the 
unpolarized real photons in DlSe-y experiments. The results of analysis of the hadronic 
final state accompanying these measurements are collected together and presented as 
a separate issue. We discuss also results on double-tag hadronic events (DlSe-y*)) 
single and double-tag leptonic events. The results from resolved real and virtual photon 
processes at HERA collider based on the single and double jet events, and in addition 
on charged particles and prompt photons, are presented. In the context of virtual 
photon processes the HERA data for forward jet and forward particle production are 
included. In addition a short presentation of the very recent data on the heavy quark 
content of the photon is given to signal the new type of measurements. Related topics 
- the polarized photon structure functions, the structure function of the electron and 
the photonic content of the proton - are also shortly mentioned. In the Appendix 
we give a short description of the existing parton parametrizations for the real and 
the virtual photon. Bibliography contains, besides references to the data we describe, 
references to MC generators and the comprehensive list of theoretical papers divided 
into groups according to the topics. 

We will not discuss the total cross sections, either for real or virtual photons, nor 
diffraction phenomena. The data from the dedicated studies of the structure of jets 
and remnant jets are omitted, a compromise between the completeness of the picture 
and the size of the presentation. 

Each of the topics presented in the survey is preceded by a short review of the 
basic theoretical ideas, the description of the specific - for the discussed topic - type 
of analysis, used MC programs, the quality of the description in terms of the MC and 
theoretical predictions, and finally the notation. (In different types of measurements 
and in analyses made by different groups, different notation is used. We try to keep 
a consistent notation throughout the survey. The double role played by the virtual 
photon - serving as a probe of the structure of the target, and in some cases being a 
target itself, which may be resolved by final large pt hadrons/jets, causes additional 
problems, e.g. versus notation for DISep events at HERA.) In each entry we 
present individual measurements by the short description of results together with some 
representative figures, as well as comments/conclusions mainly quoted from original 
publications. (Standard names of the Monte Carlo generators and the acronyms of 
the parton parametrizations are used in the descriptions, see references and Appendix, 
respectively.) We also give tables with the measured values of the hadronic structure 
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function F2 and for parton densities for the real and virtual photon. 

We denote each entry by the name of the collaboration in alphabetic order, the 
publication year (trying to keep the chronological order) and the reference number 
listed in the bibliography, and in addition by the name of the collider. The conference 
presentations or submitted papers are labelled with conf in the header of the entry. 
By adding a small letter we distinguish among the entries published or presented by 
the same collaboration in the same year, in the alphabetic order corresponding to their 
appearance in the text, e.g. DELPHI 96a. 



7 



2 Partonic content of the real photon 



The deep inelastic scattering on the 7 , DISe7, with the hadronic final state produced in 
the single-tag events, corresponding to the 7*7 collision at e^e~ colliders, is the main 
source of information on the partonic content of the real photon. For the unpolarized 
e^e~ case, the hadronic structure functions Fi^2 can be introduced to describe this 
process, in analogous way as for the proton in the DISep experiments. 

Contrary to the proton case, the (Quark) Parton Model gives definite predictions 
for F^2i based on the fully calculable QED process 7*7 — > qq. These functions, and re- 
lated quark densities in the photon, depend logarithmically on Q^. The Parton Model 
contribution together with QCD corrections constitute the inhomogeneous evolution 
equations which can be solved without any extra input, leading to the definite QCD 
predictions both for xsj and dependence. This was why the photon structure 
functions were treated at the beginning as a unique test of QCD, allowing also, in 
principle, for an extraction of Aqcd- The pointlike (called also asymptotic) solution 
obtained in this way suffers, however, from the singularities at small xbj, calling for 
an additional contribution of nonperturbative (hadronic) origin. Even after proper 
implementation of the hadronic contribution the specific properties of F2 , namely the 
large quark density at large xbj and logarithmic rise with Q"^ remain, providing a very 
well defined (a "must") test of QCD. 

First measurements of F2 were performed at PETRA and PEP in early eighties 
ll. Recently new data from LEP (LEPl with the CM energy y/s ^ 90 GeV, LEP 1.5 



with 130-136 GeV and LEP2 with 161-202 GeV) and TRISTAN (with the CM energy 
58 GeV) e'^e~ colliders have appeared. They reach very small xbj ~ 10~^ (at small 
Q^), and also very high < >~ 400 GeV^ regions. In sec. |2.2.2| we collect these 
results for the hadronic structure function F2 of the real photon. 

In extracting F2 the unfolding of the true kinematical variables is a real challenge. 
Specific problems which occur here, special data analyses and related experimental 
results are discussed in the separate section (sec. The contributions due to the 

resolved, both target and probe, photons happen to be also of relevance for proper 
unfolding of F2, see below. On top of it the target is always slightly virtual, which 
infiuences the extraction of F2 ■ The problem how to correct for the offshellness of the 
target photon is easier to handle for a leptonic structure function p^^^^^^ ^ based on 
the process 7*7 —>■ , see corresponding results in sec. ^. 

The description of the heavy quark content of the photon deserves a comment. 
At very high scale (~ Q^) the c and b-quark, or even t-quark play a role of another 
(massless) parton in the photon. At the present relatively low scales it is often enough 
to use just the (Q)PM description of the massive c (and b) - quark contribution to 
F2 . In early measurements heavy quark contributions were even subtracted from data, 
as they are not given by the pointlike solution of the QCD evolution equation. It is 
only recently that the charm part of F2, F2^^, was measured at LEP in the dedicated 
experiment. These and other very recent results on the heavy quark content of the 
photon are shortly commented in sec. ^ 

The other source, of growing importance, of information on the partonic content 
of the real photon is the hard particle(s) or jet production in 77 collision at the 
e+e~ colliders or •jp collision at ep HERA collider {^/s ~ 300-320 GeV). These hard 
processes were observed first at TRISTAN and HERA. Here, as in a typical hard 



8 



hadron-hadron process, the individual partons of parent particles {e.g. two photons 
in e'^e~ or photon and proton in the photoproduction at ep collision) interact and 
the corresponding parton densities in the photon may be extracted from data. The 
collection of the results for these resolved real photon processes is presented in sec. 
27^ . For corresponding results for a virtual photon see sec. |^. 

Neither the structure functions of the polarized real photon nor the corresponding 
resolved polarized photon processes have been measured so far; sec. ^]T] contains a 
short note on this topic. 



2.1 Source of the real photons 

In the present experiments there are no beams of high energy real photons. Therefore 
testing of the "structure" of the real photon is performed at the e~^e~ collisions, both 
in the DlSg-y and the resolved 77 (with both photons real) processes, or in photopro- 
duction at e^p colliders. In all these cases the flux of 'initial real' photons arises from 
the electron or positron beams. The flux of these (almost) real photons in present 
experimental setups for unpolarized electron (positron) beams can be approximated 
by the Weizsacker - Williams (WW, sometimes called also the Equivalent Photon 
Approximation, EPA) formula [T^, |^ : 



mm 



mm 



where z denotes the fraction of the electron energy taken by the photon. P, 
mgZ^/(l — z) with me denoting the electron mass, and P^^x stands for the maximal 
value of the photon virtuality appropriate for the experimental setupQ. Eq. |I] gives the 
typical soft bremsstrahlung spectrum. Note that in 67 or 77 options planned for the 
future Linear Colliders, where the real photons would be produced in the Compton 
backscattering process, the corresponding flux is expected to be much harder |]16| . 



Here also a flux of photons with definite polarization can be obtained in an easy way. 



2.2 DISe7 experiments 

In this section we consider the standard DlSe-y measurements in unpolarized 67 colli- 
sions based on the process (fig. 0): 

e(A;)7(p) — * e(A;') hadrons, (2) 

at e~^e~ colliders (single tagged events) 0. Here the target is a real (in practice almost 

real, with p"^ ~ 0) photon, radiated by the initial electron^. 

In the DISe7 Process (|^) the real photon with four-momentum p, and the virtuality 
= —p^ = 0, is probed by the highly virtual photon with four-momentum q = k — k', 

where k and k' are the four-momenta of the initial and final electron, and a large 

virtuality = -g^ >0.g 

^The critical discussion of this approach for the 77 processes in e+e~ colhsions at LEP 2 can be 
found in [T5[ |. 

^For simphcity we will use the symbol e for both electron and positron. 

^The positive square of the four-momentum of the exchanged particle will be called its virtuality. 
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y(pio) 

Figure 1: Deep inelastic scattering on the (almost) real photon 67 (P^ = — 0) ^ eX. 
The photon probe has a large virtuality = —q^ >0. 

The differential cross section for tlie process (|^) is given by: 

'-[{l-y)F] + XB,y^F2] (3) 



dxBjdy 

\{l + {l-y?)F,^-y'F2], (4) 



witli tlie relation among tlie transverse Fj, (= F^) and longitudinal F2 structure 
functions: = F^ + 2xBjFi [^. The standard DIS variables are defined by 

p-q 

^Bj = , y = r, (5) 

2p ■ q p ■ k 

where xbj is a standard Bjorken variable, while y, in case of a massive target [e.g. 
proton) denotes the scaled energy of the exchanged photon in the target rest frame. 
In practice, in DlS^-y measurements the variable y is small and the cross section (|^, ||) 
is effectively saturated by F2. Note that in e+e~ collisions there is always a small off- 
shellness of the photon target, i. e. = —P^ 7^ 0, and there may appear in addition the 



third structure function , |T^. After integration over the azimuthal angle between 
two scattering planes of the scattered and e~ it disappears however from the cross 
section in the single-tag measurements, since it enters the cross section as the term 
F3'cos(20)). See sec. |] for a discussion on virtual photon measurements where is 
kept far from zero. Discussion on the additional structure functions measured for a 
leptonic final state can be found in sec. 



2.2.1 Theoretical description 

In the (Quark) Parton Model (denoted as QPM[]or PM) one assumes that the hadronic 
final state in eq. (H) is due to the production of quark pairs: qi and (ji {i=l,2...Nf, 
with Nf being the number of quark flavours) with the fractional charge Qi. The F2 
is obtained by the integration over the four-momentum of the outgoing quarks with 
respect to the target photon direction. The full (QPM or the lowest order) expression 

^Usually in QPM quarks are treated as massive. 
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for F2', keeping the terms with quark mass rriq., is given by the Bethe-Heitler formula 



^ i=l ^ 



+[x% + (1 - XBjf + XBj{l - 3XB,)^ - 4,^] In [^], (6) 
where the quark velocity f3, in the 7*7 CM system, is given by 



/5 = \ 1 



4m^ 

Hi 



1 - -^2l£i5L_. (7) 



s \ Q\l-XBj)' 



The squared energy of the 7*7 collision, s, equal to the square of the invariant mass 
of the hadronic system, W"^, is given by 



s = W' = ^{l-XB,). (8) 

XBj 

In the limit of s well above threshold, i.e. for /3 ~ 1, which for fixed XBj (not too 
small and not too close to 1) corresponds to the Bjorken limit, one gets 

In^ ^~ln ^ 9 

1-/5 "^l^Bj 

The structure function F2 can be approximated in such case by 
a 

F^ = -N, J2 Qt^BjiK, + (1 - ^BjY] In — + SxBjil - XBj) - 1] (10) 

^ i=l ''^qi 

and it can be used to define the quark densities in the photon: 



2N 



F2 = XBjY.QhnxB,,Q'). (11) 



In the above formulae Nc=3 denotes the number of colors and Qi - the fractional 
quark charge. In the last equation a (natural) assumption, that quark and antiquark 
distributions in photon are the same, has been introduced. Note that 

• F2 is calculable in the (Q)PM, in contrast to the structure function of hadrons, 
e.g. the nucleon structure function F^ . 

• F2 is proportional to a, the fine structure coupling constant (~ 1/137). There 
is an overall logarithmic dependence on the energy scale squared (or Q^, see 
eq. (|)). A large value of F2 , or in other words a large quark density, is predicted 
for large XBj- 

• F2 is a sum of the quark contributions, each containing the factor Qf, so the 
individual quark density is proportional to (see eqs. (|T0| , [TT|) ). 
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The longitudinal structure function Fl is not zero in PM, in contrast to the corre- 
sponding function for hadrons, and it is scale invariant: 

Nf 

PI = -iV,5:QfxB,[4xB,(l - XB,)]. (12) 

i=l 

Note that also the third structure function is scale invariant: 

Nf 

Fl = -NcY.Qt^BA-^l,]- (13) 
^ i=l 

In the leading logarithmic approximation (LLA) the Callan-Gross relation F2 = 
2xBjFi holds as in the case of hadrons. 

In this approximation the PM formula for the quark density is given by 

qUxb,, Q')\^pis = ^NM^% + (1 - xs,f] In (14) 

Since the LL contribution corresponds to the on-shell quarks one can treat XBj in the 
above formula as equal to the part of four-momentum of the initial photon carried 
by the quark; this latter variable is usually denoted by (or simply x). Note also 
that in the above formula instead of a quark mass there appears the QCD scale Aqqd- 
Therefore this expression allows to describe all the light quark (light as compared 
to the scale W"^) contributions to F2 in a universal way. Heavy quark contributions 
should be treated separately, according to the QPM (eqs. ^ and |^), as long as is 
not considerably larger than [Q. 

An additional dependence will appear in F2 and further in qj due to the QCD 
corrections. Here also the gluonic content of the photon, C^, appears. Large logarith- 
mic corrections ~ InQ^ can be described by the Dokshitzer-Gribov-Lipatov-Altarelli- 
Parisi (DGLAP) type of equations |5|, |l|] or by other techniques [|, |l||, |0|, 0, |2 



m 



the leading logarithm approximation LLA, in the next-to-leading logarithm approxi- 
mation NLLA, or with a higher accuracy (NNLLA,...). 

The inhomogeneous DGLAP equations for the (massless) parton densities in the 
real photon can be represented in the following way (below is used for both quarks 
and antiquarks): 



dq]ix,Q') a 



with the standard splitting functions Pgq, Pqg, Pgq, Pgg and in addition with the function 

Pg^ = N,[x^ + {l-x)\ (16) 
describing the splitting of the photon into quarks (from eq. 14). 
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The dominant LL (Q^) contribution in a general solution of the eqs. ( P3[ ) (corre- 
sponding to the one-loop expressions for and for the splitting functions) is coming 
from the strong ordering of the transverse momenta of radiated gluons along the chain 
between the target and the probe particles. 

Note that it is possible to solve the above equations, based on the pointlike cou- 
plings only, without the initial conditions, assuming the LL (Q^) or NLL (Q^) be- 



haviour of the solution [|1^. Obtained in this way the, so called, asymptotic solutions. 



which give definite predictions for both the xbj and Q dependence of the q^, G"' and 
F2. The logarithmic dependence of this solutions on the QCD parameter Aqcd would, 
in principle, allow to extract this basic parameter from F2 data. However these point- 



like solutions have a singular behaviour at small x |2^, and in practice while solving 
the equations (14) the initial conditions have to be assumed from a model {e.g. VMD) 
or taken from measurements, at some (low) scale [^. This way nonperturbative 
contributions enter, which regularize somehow the unphysical behaviour of the purely 
pointlike solution. 

At this stage the parton parametrizations for the photon, can be constructed. For 
light quarks and a gluon the procedure is based on eqs. (p!5D, while treatment of the 
heavy quarks depends on the scale. At present energies (or Q"^ scales) for c, b quark 
contributions the QPM formula (eq. ^), with possible QCD corrections, is applied 
close to the threshold, see also sec. ^. In the Appendix some details of existing parton 
parametrizations for a real photon are given. 

In pre-QCD times, the hadronic structure of the photon has been solely attributed 
to the vector meson component (p, uj, (p) in the real photon. Therefore for the matrix 
element of the photon between e.g. the nucleon states, the following representation 



by the corresponding matrix elements for the p meson current was assumed 

Tfl^ 1 

<N\J^\N>= <N\p^\N> +... (17) 

9p r - 

That was the basic assumption of the Vector Meson Dominance (VMD) model or, if 
higher than (f) vector meson states were included, of the Generalized VMD (GVMD)|3 
21, n]. Once a photon "becomes" a hadron, the cross section for 77 hadrons 



is dominated by soft exchange processes, as for any other hadron-hadron collisions 
{soft VMD contribution). This mechanism is also present in DlSe^ case, leading to 
nonperturbative (or hadronic) contribution to the measured structure function F2 , as 
well as to or densities. 

It is clear that the VMD not only models the soft photon- induced processes, but, 
being built into the parton densities, it is also important in the description of the hard 
large pt processes involving resolved photons. For the general discussion of various 
components of the photon, with less conventional description of the photon-induced 
processes and the photon structure functions, see ref. and the Appendix. 

Parton distributions of the photon do not fulfill the, standard for hadrons, "mo- 
mentum" sum rule formulated for quarks and gluons. The reason for this is that 7 
can interact directly, being a parton itself. Various attempts to introduce a useful! 



momentum sum rule for the photon were undertaken, see p6|, ||162|| and ||169 |. 
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denoted also VDM and GVDM 
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Finally we would like to mention specific problems when counting the order of 
the perturbation (leading logarithms) in the higher order calculations for the photon- 
induced processes. This is already seen in eq. where the QPM formula contains 
both LL and NLL terms, both being, however, of the QED origin. The presence of the 
inhomogeneous term in the evolution equations for the quark densities in the photon 
( [T^) leads to further complications in the interpretation of what is meant by the LO 
and the NLO QCD analysis. For a recent discussion, see [|^, |135|| , and also comments 
in sees. 2.4.1 and 2.4.5[ 



For future discussion (sec. |3.4|) let us mention other than the DGLAP approaches 
to the deep inelastic scattering processes (for any target). They may become impor- 
tant if the region of small xbj is tested in DIS experiments. In such case large loga- 
rithms ~ ln{l/ XBj) appear, and summing them to all orders becomes necessary. The 
corresponding evolution equation, Balitsky-Fadin-Kuraev-Lipatov equation (BFKL), 
leading to LL [I/xej) terms does not rely on the strong ordering in the transverse 
momenta along the chain, allowing for the parton emission with large transverse mo- 
menta [^0- The combined approach was also proposed (Catani-Ciafaloni-Fiorani- 



Marchesini equation) . Note that BFKL and CCFM approaches should apply for 
any target as the DGLAP does, since they correspond to the QCD chain to be attached 
to the target particle. 

2.2.2 Measurements of F2 

Measurements of the F2 {xbj, <5^) for the real photon base on the single tagged events 
at e~^e~ colliders {i.e. DlSe-y events, see eq. (H) and fig. |l]). The bulk of the hadrons is 
produced at small angles to the incoming e~^e~ beams, and remains undetected. Due 
to the particle losses in the experiment not xbj but the quantity Xyis is measured, 

Q2 

^ Q' + wL' ^^^^ 

where Wms is the invariant mass of the visible hadronic system. Since the target 
photon has always the nonzero virtuality P^, the measured quantity is in fact 

Q2 

= Q2 + WL+P^- ^^^^ 

In order to reconstruct xbj from a measured value for Xms, an unfolding procedure is 
needed. To this end, in some analyses the related variables, Wrec, Xrec, are also used, see 
e.g. L3 98a, OPAL 2000. They are reconstructed from the measured four- momenta 
of electrons with the constraints from the transverse momentum conservation. Even 
using the above technique, the value of Wrec is still generally smaller than the true 
value due to e.g. energy loss in the forward regions. Therefore new variables Wcor 
and formed to correct for this fact, see OPAL 2000. The very important 

improvement in the description of the data has been achieved due to the recently 
proposed two-dimensional unfolding, with the energy in the forward cone used as the 
second variable 0. It was used in ALEPH 99b,conf, and OPAL 2000 leading to 
much smaller errors 0, ^ . The proper inclusion of the nonzero virtuality of the target 



photon^ should lead to the further improvement (see sec. ^ where this effect for the 

^^The NLO BFKL calculations demonstrate large corrections to the LO BFKL results. 
^^The procedure "is not well defined if the photon is not pointlike" |lC|]. 
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leptonic structure functions is discussed). The unfolding can be performed using the 
traditional method based on the linear scale or using the approach with the logarithmic 
scale, especially useful for extracting F2 in the small Xsj region 

In order to perform reliable unfolding of the structure function F2 from the cross 
section for 7*7 — »• hadrons, the corresponding events are grouped into classes of similar 
topologies (see discussion in [§, § and e.g. TOPAZ 94, OPAL 97a, DELPHI 
96b,conf), which roughly coincide with the contributions from QPM, VMD and RPC 
(RPC - Resolved Photon Contribution). First analyses applied the FKP approach, 
with the pointlike QCD contribution assumed for light quarks, QPM description of 
the heavy quarks, and VMD contribution for the soft hadronic production. The 
parameter separates the pointlike and the hadronic (p-type) configuration of the quarks 
in the target photon (at the first vertex level). The additional parameter p™*", in the 
FKP approach of the same order as p^, is introduced as a cutoff for the perturbative 
2 — > 2 subprocesses QCD calculations. More details on the modelling of the hadronic 
final state can be found in sec. |2.3| . 

Nowadays the Monte Carlo generators HERWIG and PYTHIA, adapted for DlSe^ 
in 1995 (see e.g. OPAL 97a), can describe these contributions (using the cutoff 
p^*"-) with any of the existing parton parametrizations for the photon. It is no longer 
necessary for the resolved 7 term to fit an empirical pj, parameter to the data before 
unfolding, as in the parton density in the photon both nonperturbative {e.g. a la VMD) 
and the perturbative contributions are included. The PYTHIA and modified HERWIG 
("HERWIG + power law-p^") programs with broader pt distribution dpl/p^ + pIq 
originally introduced to improve description of 7p processes (ZEUS 95c), seem to 
describe 77 data better than the default versions {e.g. |^). On the other hand, the 
TWOGAM generator, based on QPM + VMD + RPC contributions is successfully 
used to describe the newest F2 data although with a large p^*" parameter, ~ 2.3 
- 3.5 GeV, instead of the default value 1.8 GeV (L3 99a). The general purpose 
generator PHOJET based on the Dual Parton Model is also used successfully in the 
F2 measurements (with p™" = 2.5 GeV). The large dependence on the chosen Monte 
Carlo model leads to large errors in the F2 data points, or even to two sets of data 
points from the same event sample (L3 98a) . The multi-dimensional unfolding method 
reduces the dependence on MC models, as it was found in very recent analyses (OPAL 
2000 and ALEPH 99b,conf), see sec. 

The F2 measurements at small xsj are presented in DELPHI 96a, OPAL 97b, 
OPAL 2000 and L3 98a, with the lowest value xbj (center of bin value) 0.001 
(0.0022). The low Q"^ measurements at LEP (< >~ 2 GeV^) were reported in 
OPAL 97b, L3 98a, OPAL 2000 while the largest < >=400 GeV^ data are 
discussed in DELPHI 98,conf. Note also that in the earlier measurements the c 
(and h) quark contribution was usually subtracted from the structure function F2 
{e.g. OPAL 94). For the dedicated measurement of F2C reported in 1999 by OPAL 
collaboration, see sec. |^. 

The general features of F2, as far as xbj and logQ^ dependences are concerned, 
agree with the theoretical expectations, although the precision of the data does not 
allow in many cases to distinguish between existing parton parametrizations and clarify 

^•^The parameter pto ~ 0.66 GeV should not be confused with separating the soft and hard 
hadronic processes within the FKP approach. 
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the small xbj behaviour of F2 . The parton densities which give, at present, the best 
description of the F2 data are GRV and SaSlD [0. 

As we have already pointed out, the study of the hadronic final state together 
with the possibly resolved (real and virtual) photon contribution became part of the 
measurements of the structure function for the photon at e~^e~ colliders (details of the 
studies of the hadronic final state will be presented in the next section). Here we would 
only like to stress that some discrepancies have been found for certain distributions, 
e.g. the hadronic energy flow. This fact is included in the estimation of the uncertainty 
of the measured function F2 . 

We start the presentation of the F2 data from the LEP collider, then the TRISTAN 
data are discussed. Figures and tables with corresponding numbers for the F2 as a 
function of xbj at fixed < >, and also for F2 averaged over relevant xbj ranges 
are given. Collective figures of F2 versus xbj and F2 versus containing all existing 
results are presented at the end of this section. 
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DATA 



•ALEPH 99a |]3§ (LEP 1) 

Data on F2 for xsj from 0.005 to 0.97 and between 6 and 3000 GeV^, in three 
bins, were collected in the period 1991-95. An analysis of the hadronic final state 
was performed using the QPM+VMD model and the "HERWIG 5.9 + power law 
involving GRV LO parton parametrizations (see next section for details). 

The unfolded results for F2 as a function of xbj are presented in table ^ and in 
figs. H, ^, where the comparison with predictions of the parton parametrizations: 
AFG HO, LAC 1 and GRV LO is made. The values of F2 averaged over XBj (for 0.1 
< XBj < 0.6) for three values of < > are given in table g 



Table 1: 





XBj 




[GeV^] 




{stat. + syst.) 


9.9 


0.005 - 0.08 


0.30 ±0.02 ±0.02 




0.08-0.2 


0.40 ±0.03 ±0.07 




0.2 - 0.4 


0.41 ±0.05 ±0.09 




0.4-0.8 


0.27 ±0.13 ±0.09 


20.7 


0.009 - 0.120 


0.36 ±0.02 ±0.05 




0.12 - 0.27 


0.34 ±0.03 ±0.11 




0.27-0.50 


0.56 ±0.05 ±0.10 




0.50-0.89 


0.45 ±0.11 ±0.05 


284 


0.03-0.35 


0.65 ±0.1 ±0.09 




0.35-0.65 


0.70 ±0.16 ±0.19 




0.65-0.97 


1.28 ±0.26 ±0.26 



ALEPH (Q^) = 9.9 GeV 




Figure 2: The unfolded structure function F2 /a as a function of XBj for < >=9.9 
GeV^. Comparison with predictions of different parton parametrizations: AFG HO (solid 
hne), LAC 1 (dashed hne), GRV LO (dot-dashed line) (from 
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Figure 3: The same as in fig. I for < Q2 >=20.7 GeV^ (from 
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Figure 4: The same as in fig. | for < >=284 GeV^ (from j^). 



Table 2: 



<Q'> 
[GeV^] 


< F^/a > 
(tot) 


9.9 


0.38 ±0.05 


20.7 


0.50 ±0.05 


284 


0.68 ±0.12 



The dependence of F2 (from the 0.1 < xbj < 0.6 range) has also been measured 
(not shown) and a good agreement with the OPAL 97c measurements was found. 

The F2 data were compared to different parametrizations and the test was 
performed. LACl and 2 and WHIT 4,5 and 6, all containing a large gluon density at 
low X, show large x^- 

Comment: "Comparisons to parametrized parton density functions show that those 
containing a large gluon content are inconsistent with the data. " 

•ALEPH 99b,conf g (LEP 2) 

The measurement of F2 at A/i=183 GeV based on the 1997 data was performed using 
a two-dimensional unfolding method (the principle of maximum entropy) 0. The two 
g2 ranges: 7-24 GeV^ and 17-200 GeV^ were studied. In the analysis "HERWIG 5.9 ± 
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power law pt" with the GRV LO and SaSlD parton parametrizations and in addition 
the PHOJET model were used. 

The hadronic final state analysis was performed (see next section for details) and 
it was found that the unfolding procedure reduces the dependence on Monte Carlo 
models, and that "PHOJET is in poor agreement with the data for the higher 
range" . 

The unfolded results for F2 as a function of xsj are presented in fig. |^ for < Q'^ >= 
13.7 GeV^ and 56.5 GeV^ and in table | (only for lower Q'^). 

PRELIMINARY 



] <Q->=13.7 (GeW/cf 






ALEPHdata 




GRVLO 




SaSld 




LACl 





0.1 0.2 0.3 0.4 0.5 0.6 0.7 O.S 0.9 



1 <Q->=56.5 {GeW/cf 









0.1 0.2 0.3 0.4 0.5 0.6 0.7 O.S 0.9 



Figure 5: The unfolded structure function F2 /a as a function of x Bj for < Q"^ >=13.7 GeV^ 
and 56.5 GeV^. Comparison with predictions of different parton parametrizations: GRV LO 
(sohd hne), SaSlD (dotted line), LACl (dot-dashed line) (from ^). 



Table 3: 



<Q'> 
[GeV^] 


XBj 


(stat. + syst.) in % 


13.7 


0.003 - 0.023 
0.023 - 0.092 
0.092 - 0.213 
0.213-0.386 
0.386 - 0.786 


0.28 ±5.3 ±6.3 
0.34 ±7.5 ±5.3 
0.34 ± 10.8 ±8.0 
0.35 ± 13.6 ± 15.7 
0.43 ± 11.1 ± 10.5 



Comment: The new unfolding method was used in the analysis. It "leads to smaller 
statistical errors and a reduced model dependence compared to one- dimensional proce- 
dure". The GRV LO parametrization describes the F2 data well. At low x, "the SaSlD 
parametrization is found to be lower than data, whereas the LAC 1 parametrization is 
considerably higher. . . " . 

•DELPHI 96a [|5|] (LEP 1) 
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Data on F2 were taken in the period 1991-93 for between 4 and 30 GeV^ (< > 
= 12 GeV^) and for xsj down to 0.003. The so called p^'^'^^^^ was also measured for 
nonzero virtuality of the target photon (see sec. [4.2| for details). Estimated target 
photon virtuality was used in the unfolding of F2 0. The TWOGAM event gen- 
erator was used to simulate QPM events, another event generator was used to obtain 
the QCD correction (LL) to the pointlike contribution for light quarks in the FKP 
approach. The GVMD and the pointlike (FKP) contributions were studied with 
= 0.1 and 0.5 GeV. Results for F2 for the light quarks are presented in fig. |^ and in 
table |. 

Q 0.8 

a 

^ 0.7 

Q 
O 

a 

0.5 
0.4 
0.3 
0.2 
0.1 



"0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

X 

Figure 6: The unfolded F2 /a for the hght quarks together with OPAL 94 data. The curves 
show the sum of the GVMD model prediction multiplied by the threshold factor I-xbj SLnd 
the prediction of the FKP parametrization, with the parameter pj, equal 0.1 GeV (upper 
line) and 0.5 GeV (lower line). The bottom curves show separately the GVMD contribution 
with different target masses (from |3^j. 




Table 4: 



<Q'> 
[GeV^] 


XBj 


F^/a 
[stat. + syst.) 


12 


0.003 - 0.080 
0.080 - 0.213 
0.213 - 0.428 
0.428 - 0.847 


0.21 ±0.03 ±0.06 
0.41 ±0.04 ±0.05 
0.45 ±0.05 ±0.05 
0.45 ±0.11 ±0.10 



The averaged value of F2 /a over the xbj range between 0.3 and 0.8 was extracted: 



<Q'> 
[GeV^] 


< F^/a > 


12 


0.45 ±0.08 



For comparison with other measurements of the dependence of the F2 , see fig.0. 
It was found that although < >=0.13 GeV^, a fixed value of 0.04 GeV^ fits the data better. 
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Figure 7: F2 /a averaged over xbj between 0.3 and 0.8 as a function of Q^. The curves show 
the FKP parametrization predictions for different values of the parameter p% (from ,35]). 



The study of the F2 behaviour at < >=12 GeV^ in the low xbj domain leads 
to the following results (fig. || and table H): 



Table 5: 



<Q'> 
[GeV^] 


XBj 


[stat. + syst.) 


12 


0.003 - 0.046 
0.046 - 0.117 
0.117-0.350 


0.24 ±0.03 ±0.07 
0.41 ±0.05 ±0.08 
0.46 ±0.17 ±0.09 



^ 0.6 
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0.05 O.I 0.15 0.2 0.25 0.3 0.35 

X 

Figure 8: The unfolded F^/a for < >= 12 GeV^ as a function of xsj with the LO 
parametrizations LACl, GS, DO, and GRV (from ^Ej]). 

Comment: The effect of the nonzero was included in the unfolding o/F^. No rise 
of F2 at small xsj has been found. The GRV and GS leading order parametrizations 
of the quark density in the photon are in agreement with the data. 

•DELPHI 96b,conf [|§ (LEP 1) 

The measurement of the photon structure function F2' at < > = 13 and 106 GeV^ 
(data collected in the years 1991-95) for xsj down to 0.003 together with a study of 
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the hadronic final state (see next section) is reported. Three types of subprocesses 
have been considered: QPM (direct term with rriq ^ for all quarks), soft hadronic 
GVMD or VMD (the TPC/27-type parametrization) and the resolved target photon 
contribution (RPC) with GS2 parton parametrization. In addition the effects of the 
resolved probe photon were studied. The unfolding was done in the linear, and, for 
better sensitivity to Iow-xbj, logarithmic scale in xsj- 

Fig. 1^ shows the xsj dependence of F2 for two samples, compared with the 
predictions of various parton parametrizations. 




Figure 9: Unfolded /a as a function of xbj compared with QPM+GVMD+RPC (with 
parton parametrization GS2) (solid line) and predictions of the SaS4 (dashed line), GS2 
(dash-dotted) and GRV3 (dotted) parametrizations (from 

The value of F2 /a averaged in the range 0.3 < xbj < 0.8 was extracted for = 
13 and 106 GeV^ (see table |). 



Table 6: 



<Q'> 
[GeV^] 


< F^/a > 
(stat. + syst.) 


13 
106 


0.38 ±0.031 ±0.016 
0.576 ±0.081 ±0.076 



Comment: The importance of the final hadronic state topology was noticed and the 
resolved photon contribution was included. The study of a linear and logarithmic un- 
folding was performed. 
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•DELPHI 97a,conf (LEP 1, LEP 2) 

A study of F2 in the range between 3 and 150 GeV^, based on data from the 1994- 
95 runs for energies around the Z° mass and from 1996 for energies between 161 and 
172 GeV is presented. An analysis of the hadronic final state is performed using the 
TWOGAM Monte Carlo program, where QPM (with A^/ = 4), soft hadronic VMD and 
RPC (with the GS2 parton parametrization, both single and double resolved photon 
terms) parts are included (see also next section). 



The unfolded results for F2 as a function of xsj and are presented in figs. |TO]a 



and p!0|b, respectively. 
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^ LEPl 




L2 L4 L6 L8 
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Figure 10: a) F2/a versus xbj from the DELPHI experiment based on the LEP 1 
data for three values of < >: 6.5, 13, and 22 GeV^. The soUd hne corresponds to 
QPM+GVMD+RPC (with with parton parametrization GS2), dotted - QPM+GVMD and 
dashed - GVMD. b) F^/a averaged over XBj from 0.3 to 0.8 as a function of Q^. Results ob- 
tained from the LEPl and LEP2 data are shown together with results of a previous analysis 
of the LEPl data (from j. 



Comment: "There is some indication of deviations of TWOGAM predictions for the 
data in the high region. " 



►DELPHI 98,conf |^ (LEP 2) 



Data on F^ were taken in the period 1996-1998 for between 10 and 1000 GeV^ 
(for the energy 163-188 GeV). The three components were applied: the QPM for all 
quarks simulated by the TWOGAM generator, the soft hadronic VMD contribution 
and RPC for single and double resolved photon processes, and a good description was 
obtained (see also next section). The data are shown in table |^ and in fig. [Tl|. 
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Table 7: 







In 

{ ^fnf A- '^71 ) 


91 


U.U i — U. i 


u.oo m u.ui m u.uo 




U. i — u.o 


n zLi -1- n (T? 4- n 09 
u.^i Hz u.uo m u.uz 




3 — 08 


51 + 05 + 04 


49 


n ni n 1 

U.U i — U. ± 


n 41 -1- n ni -i- n n'? 
u.^i IE u.ui m u.uo 




u. i — u.o 


n 48 -1- n 09 -1- n 09 
u.^o Hz u.uz ZL u.uz 




0.3 - 0.8 


59 ± 03 ± 04 


99 


0.01-0.1 


0.45 ±0.06 ±0.02 




0.1-0.3 


0.52 ±0.05 ±0.03 




0.3-0.8 


0.73 ±0.05 ±0.03 


400 


0.01-0.1 


0.5 ±0.3 ±0.1 




0.1 - 0.3 


0.70 ±0.2 ±0.2 




0.3-0.8 


1.0 ±0.1 ±0.3 




Figure 11: The DELPHI structure function F2 /a as a function of xsj at < >= 22, 43, 
99 and 400 GeV^. Comparison with the parton parametrizations: GS (dash - dotted line), 
GRV (dashed hne) and SaS (dotted line); TWOGAM prediction is the solid line (from §^). 
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The dependence was studied in two xsj ranges: 



for 0.01 < XBj < 0.1 Fa^/a = (0.12 ± 0.02) + (0.073 ± 0.009) InQ^ 

while for 0.3 < xbj < 0.8 F:j/a = (0.09 ± 0.05) + (0.14 ± 0.02) InQ^. 

The data for F2 /a as a function of are presented in figs. [12], |13| in comparison with 
the Monte Carlo predictions and LO calculations. A suppression factor was introduced 
to describe the high data. 
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Figure 12: The DELPHI results for F2 / a as a function ofQ^ forO.Ol < xbj < 0.1 compared 
witli otiier data (from p^). 
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Figure 13: The DELPHI results for ja as a function of for 0.3 < xbj < 0.8 compared 
with other data (from J3^j. 
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Comment: "Prediction for the RFC component tends to have too high cross section 
for the high Q^" and the low XBj- "For the first time TWOGAM has been tuned to 
get a good agreement in the region of high " . 



•L3 98a gD|] (LEP 1) 

The F2 data for 1.2 < < 9.0 GeV^ collected in years 1991-95 are reported. Two 
bins: < Q'^ > = 1.9 GeV^ in the xbj range 0.002 < xsj < 0.1 and < Q'^ > = 5 GeV^ 
in the xbj range 0.005 < xbj < 0.2 were analysed. To improve the measurement of 
Wyis the kinematics of the tagged electron was included (the transverse momentum 
conservation). In the analysis three Monte Carlo generators: PHOJET 1.05c (p™" = 
2.5 GeV), HERWIG 5.9 and TWOGAM = 2.3 GeV) have been used. They led 

to different results for xbj, Pt distributions and the energy flow (for more details see 
next section). 

Two (!) data sets for F2 based on two Monte Carlo analyses: the PHOJET and 
TWOGAM, are given. They differ up to 14% at high xsj and up to 28 % at low xbj- 
The results are presented in fig. |T^ and in table || (the values of F2 /a are given at the 
centre of the xbj bin). 



Table 8: 



<Q'> 
[GeV^] 


XBj 


F^/a{setl) 
{stat. + syst.) 


F^/a{set2) 
{stat. + syst.) 


1.9 


0.002 -0.005 
0.005 -0.010 
0.010-0.020 
0.020 -0.030 
0.030 -0.050 
0.050 -0.100 


0.184 ±0.009 ±0.013 
0.179 ±0.007 ±0.009 
0.176 ±0.006 ±0.006 
0.191 ±0.008 ±0.004 
0.193 ±0.008 ±0.007 
0.185 ±0.007 ±0.015 


0.231 ±0.011 ±0.016 
0.199 ±0.008 ±0.010 
0.191 ±0.007 ±0.006 
0.193 ±0.008 ±0.004 
0.199 ±0.008 ±0.007 
0.206 ±0.008 ±0.017 


5 


0.005 -0.010 
0.010-0.020 
0.020 - 0.040 
0.040 - 0.060 
0.060 -0.100 
0.100 -0.200 


0.307 ±0.021 ±0.035 
0.282 ±0.014 ±0.027 
0.263 ±0.011 ±0.015 
0.278 ±0.013 ±0.007 
0.270 ±0.012 ±0.008 
0.252 ±0.011 ±0.029 


0.394 ±0.027 ±0.045 
0.318 ±0.016 ±0.031 
0.277 ±0.012 ±0.016 
0.279 ±0.013 ±0.007 
0.275 ±0.012 ±0.008 
0.287 ±0.013 ±0.032 
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Figure 14: Two sets of data for ja versus XBj for two values of < > : 1.9 and 5 
GeV^. Set 1 (2) corresponds to the unfolding based on PHOJET (TWOGAM). The solid 
line corresponds to GRV LO, dashed - SaSlD (for P'^=0, is denoted here as t^) and 
dotted - SaSlD (P'^=t^=0.075 GeV^j parton distributions. For the Q'^=5 GeV^ in addition 
the LACl and LAC2 predictions are shown (from RD/j. 



The Q"^ - dependence was also studied; results for 0.01 < xbj < 0.1 are presented 
in table P (see also fig. llSD, for the two sets of data. 



Table 9: 



<Q'> 


F^/a{setl) 


F^/a{set2) 


[GeV^] 


[stat. + syst.) 


(stat. + syst.) 


1.5 


0.173 ±0.004 ±0.009 


0.196 ±0.005 ±0.010 


2.4 


0.195 ±0.005 ±0.004 


0.208 ±0.005 ±0.004 


3.8 


0.245 ±0.006 ±0.007 


0.252 ±0.006 ±0.007 


6.6 


0.278 ±0.009 ±0.013 


0.292 ±0.009 ±0.014 



The following - dependence was found for setl: 

F2{Q'^)/a = (0.131 ± 0.012 ± 0.021) ± (0.079 ± 0.011 ± 0.009) ln(gV(GeV)^ 
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Figure 15: Two sets of data for F'^ja as a function of Q^. Set 1 (2) corresponds to 
the unfolding based on PHOJET (TWOGAM). The sohd Une corresponds to the fit to set 
1, dashed hne denotes the GRV LO parton parametrization prediction, dotted - SaSlD 
(P"^ =0.075 GeV^, P2 is denoted here as t^) and dot-dashed corresponds to LACl (LAC2 is 
similar to LACl ) (from 1^). 



The F2 data are not corrected for the virtuahty of the target photon since this 
correction is model dependent: 7.5% for GVMD and 10-20% for SaSlD. As estimated 



by PHOJET and TWOGAM (QPM in TWOGAM), < >= 0.075 (0.087) GeV^ [10|. 



Comment: "Because of the large discrepancies with the data, we do not use HER- 
WIG for the results." "At low values of xbj the data are above the prediction of the 
GRV LO and SaSlD models (...). The LAC model can reproduce the Xsj behaviour of 
F2 at < >=5 GeV^ but it predicts too fast a rise of F2 as a function o/lnQ^". 
Because of the significant differences in results for F2 obtained using different Monte 
Carlo for the first time two different data sets were 'derived' from the same sample of 
events. 



•L3 99a gT]] (LEP 2) 

Results of the measurements of F2 at the CM energy 183 GeV in the interval 9 < 
< 30 GeV^, XBj from 0.01 to 0.5 are presented (data collected in 1997). Three 
Monte Carlo generators were used: TWOGAM (p™"=3.5 GeV) and PHOJET 1.05c 
^^min_2_5 QgV) with JAMVG (the latter used for the charm quark contribution, as in 
PHOJET the charm quark is treated as massless). The pj^ cutoff in TWOGAM needed 
to describe data is 3.5 GeV, since the p™" = 2.3 GeV used in previous analysis (L3 
98a) "produces a too large cross section for hard processes in the high region P5|". 
The hadronic final state is reasonably described by the used Monte Carlo generators. 
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Comparison was made with the SaSlD, GRV LO, GRV HO and LACl predictions. 
The values of F^ja are not corrected for ^ 0. 

The results for F2 /a obtained using PHOJET are given in table and in fig. [16 . 



Table 10: 





0^ ■ 

XB] 


In 


\^pv'^^ 

\\jrc- V I 




\ovtlv. oyoL.) 


10.8 


0.01 - 0.1 


30 ± 02 ± 03 




0.1-0.2 


0.35 ±0.03 ±0.02 




0.2-0.3 


0.30 ±0.04 ±0.10 


15.3 


0.01-0.1 


0.37 ±0.02 ±0.03 




0.1-0.2 


0.42 ±0.04 ±0.01 




0.2-0.3 


0.42 ±0.05 ±0.05 




0.3-0.5 


0.35 ±0.05 ±0.08 


23.1 


0.01-0.1 


0.40 ±0.03 ±0.03 




0.1-0.2 


0.44 ±0.04 ±0.04 




0.2-0.3 


0.47 ±0.05 ±0.02 




0.3-0.5 


0.44 ±0.05 ±0.11 
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Figure 16: F2 /a as a function of xbj for three bins unfolded with PHOJET (plus 
JAMVG). The data are compared with predictions of the GRV LO, GRV HO, SaSlD and 
LACl parton parametrizations (from /^/j. 



The Q dependence of F2 was studied as well, results are shown in fig. [17 
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Figure 17: F2 /a as a function of for different xbj regions from LEPl and LEP2 data 
(the unfolding based on PHOJET with JAMVG). The sohd hne corresponds to the fit, the 
GRV LO parton paranietrization prediction is denoted as a dashed hne, SaSlD - dotted, and 
dot-dashed hne corresponds to the LACl prediction (from p^/). 

The - dependence was found to be, for 0.01< xbj < 0.1: 

^2'(Q^)/« = (0.13 ± 0.01 ± 0.02) + (0.08 ± 0.009 ± 0.009) ln(gVGeV)^ 

and for 0.1< Xbj < 0.2: 

^2'(Q^)/« = (0.04 ± 0.008 ± 0.008) + (0.13 ± 0.03 ± 0.03) In(gVGeV)^ 

Comment: High p^'"" cutoff (3.5 GeV) m TWOGAM is needed. 

"The values of F2 /a are not corrected for the fact that is not strictly equal to zero". 
"The agreement between data and the different Monte Carlo predictions is good except 
at pseudorapidity values, rj >3 andrj < —1", note however that these regions contribute 
negligibly, when F2 is extracted. 

"At low XBj, the rise of the F2 (with Q"^) is larger than predicted by the GRV 
and SaSlD models, thus requiring a modification of the gluon density in the photon 
structure function F2 . " 

•L3 2000 [g2|| (LEP 1) 

The measurement of the structure function for the real photon was performed based 
on 1991-1995 data for large between 40 to 500 GeV^ with < Q"^ >=120 GeV^. 
The virtual photon structure function was studied as well, see sec. p]2|. 
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The F2 data from single-tag events were obtained using the JAMVG generator 
modelhng the QPM with A^/=4, the PHOJET 1.05c with a cutoff p?J^'^=2.5 GeV and 
TWOGAM with p??*"=3.5 GeV generating three processes: QPM, VMD and QCD 
resolved photon contributions. The hadronic final state was investigated (see next 
section for details), and it is properly described by JAMVG (QPM) model for XBj >0.5 
and by PHOJET for smaller XBj- 

The F2 data are presented in table and fig. |T8|, where the results are compared 
with the predictions of QPM and of QCD calculations, the latter being needed to 
describe the xsj region below 0.5. 



Table 11: 



<Q2> 
[GeV^] 


XBj 


{stat. + syst.) 


120 


0.05-0.20 
0.20 - 0.40 
0.40 - 0.60 
0.60-0.80 
0.80-0.98 


0.66 ±0.08 ±0.06 
0.81 ±0.08 ±0.08 
0.76 ±0.12 ±0.07 
0.85 ±0.14 ±0.08 
0.91 ±0.19 ±0.09 
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Figure 18: F2 /a as a function of xsj for < >=120 GeV^ from single-tag events. The 
solid line corresponds to the QPM prediction, the GRVLO parton parametrization prediction 
is denoted as a dashed line, AGF - dotted line, and dot-dashed line corresponds to the LRSN 
prediction (from ^]). 



The averaged values of < ^2^/0 > as a function of < > are presented in table 



121 for two XBj ranges and in fig. |TS 
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Table 12: 











{stat. + syst.) 


{stat. + syst.) 




XBj = 0.05 - 0.98 


XBj = 0.3 - 0.8 


60 


0.73 ±0.11 ±0.07 


0.66 ±0.09 ±0.06 


90 


0.89 ±0.13 ±0.09 


0.79 ±0.14 ±0.08 


125 


0.85 ±0.11 ±0.09 


0.88 ±0.12 ±0.08 


225 


1.01 ±0.25 ±0.10 


1.18 ±0.22 ±0.11 



LL 

V 



2- 



= GeV^ 
• DATA 
— QPM 
-- LRSN 



0.05 <x< 0.98 








10 

(GeV^) 



Figure 19: The dependence of F2 /a averaged over xsj = 0.05 — 0.98 for single-tag data. 
Comparison with the QPM prediction (sohd line) and the LRSN NLO calculation (dashed 
line) (from /Bl/j. 



The virtuality of the photon is estimated to be P^=0.014 GeV^. 

Comment: "The structure function -^^(a;) of real photons shows an excess at low XBj 
over QPM and over several QCD calculations. " 



•OPAL 94 g (LEP 1) 

The measurement of F2 at < >=5.9 and 14.7 GeV^ was performed using data 
from the period 1990-92. The new Monte Carlo TWOGEN program was used to 
generate events according to chosen formula for F2{x, Q^, P^) or F2 [x, Q^). The VMD 
with two different distributions of qq pair (peripheral and QPM - type) and pointlike 
QCD contributions for light quarks (in the FKP approach) were studied. The cutoff 
parameter pj. separating these contributions was determined to be 0.27±0.10 GeV for 
both values. The Monte Carlo model with this parameter value was then used to 



unfold F2 . See also sec. O 



Results for unfolded F2 , with the charm contribution calculated using QPM, com- 
pared with the PLUTO and TPC/27 data are presented in fig. 2C. 
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Figure 20: The results for unfolded charm subtracted ja shown and compared with other 
experiments at similar < > values. The curves show contributions of (left) VMD (dots), 
QCD-based model (dashes) and their sum (line) for < >=5.9 GeV^, (right)< >=14.7 
GeV^ (from g^). 



Comment: In the analysis a simple formula F2 /a = 0.2(1 —x) was applied to describe 
the VMD contribution. The TWOGEN with a QPM formula for F^{x,Q'^, P'^) was 
used to obtain the correction for the ^ 0. 
"No increase of F2{x) is observed" at small xbj- 

•OPAL 97a ^ (LEP 1) 

The measurement of F2 was done for two samples, 6 < < 30 GeV^ and 60 < < 
400 GeV^, using all data at the Z peak (years 1990-95) with new type of unfolding 
as compared to OPAL 94. The detailed analysis of the hadronic final states was 
performed and sizeable discrepancies with the expectations were found especially at 
low XBj (see also next section). The influence of the choice of different Monte Carlo 
generators (HERWIG 5.8d, PYTHIA 5.718, also for comparisons and systematic checks 
the F2GEN) on the unfolded F2 was studied. The dependence of F2 on the P^ values 
was estimated based on the SaSlD parametrization. 

The results for unfolded F2 are presented in fig. pl|a,b,c and in table |13| (the value 
of F2 /a is given at the centre of the xbj bin). 
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Table 13: 



<Q2> 
[GeV^] 


XBj 


{stat. + syst.) 


7 5 


n nm — n oqi 

0.091 -0.283 
0.283 - 0.649 


28 + O^+^-y^ 
0.32 ± 0.02t[]:?^ 
0.38 ± OMt'^ofi 


14 7 


0.137-0.324 
0.324 -0.522 
0.522 -0.836 


38 ± 01+nVo 

0.41 ± omiiofs 
0.41 ± o.o3i[J:?? 
0.54 ± o.o5t[!i^ 


135.0 


0.100 -0.300 
0.300 -0.600 
0.600 -0.800 


0.65 ± omt'^ofQ 

0.73 ± 0.08tn^ 
0.72 ± 0.10tE];^7 
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Figure 21: The F2 /a data from the OPAL experiment for the number of flavours Nf- 
Curves in (a)-(d) show predictions of the GRV and SaSlD parametrizations (from J4^j. 



The measurement of F2 /a as a function of averaged over the XBj range, 0.1< 



XBj <0.6, gives results shown in fig. ETld and in table in. 



The slope d{F2 / a) / dlnQ^ is measured to be 0.13 



+0.06 
0.04- 



Comment: New measurements of F2{x,Q'^) are presented, "allowing for the first time 
for uncertainties in the description of the final state by different Monte Carlo models". 
The F2 is not correced for . Large discrepancies between the hadronic energy flow 
data and Monte Carlo simulations are observed at low for x^s < 0.1, when the 
results are presented versus pseudorapidity or azimuthal angle (see next section). 
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Table 14: 



<Q^> 
[GeV^] 


{stat. + syst.) 


7.5 
14.7 
135.0 


0.36 ± 0.021^:^^ 
0.41 ± 0.0lin4 

o.7i±o.o7i°:J^ 



•OPAL 97b [g (LEP 1) 

The measurement of F2 (years 1993-94) was done for 1.1 < < 2.5 GeV^ and 
2.5 < < 6.6 GeV^ (average (5^=1.86 and 3.76 GeV^) as a function of x^j, reach- 
ing the lowest, at that time, measured (center of bin on a logarithmic scale) value: 
=0.004. For a better sensitivity on the low XBj region the unfolding procedure was 
performed on a logarithmic scale. Final state topology was analysed as well (using the 
HERWIG 5.18d, PYTHIA 5.722, F2GEN - both with the pointlike and the peripheral 
qq distributions generators). The GRV LO and SaSlD parton parametrizations were 
used in the analysis. Charm contribution is treated differently in different MC pro- 
grams. Discrepancy between the data (hadron energy flow) and results from the Monte 
Carlo generators, as well as between different models was found for low x^s <0.05 (see 
sec. U). 

The obtained values of F2 / a are given in table |15| and shown in fig. ^ together 
with early measurements from PLUTO and TPC/27. 



Table 15: 



<Q2> 
[GeV^] 


XBj 


{stat. + syst.) 


1.86 


0.0025- 0.0063 
0.0063 - 0.020 
0.020 - 0.040 
0.040 - 0.100 


0.27 ±0.03^^^ 

0.22 ± omto^ol 
0.20 ± omto^l 

0.23 ± OMt^Q^i 


3.76 


0.0063 - 0.020 
0.020 - 0.040 
0.040 - 0.100 
0.100 - 0.200 


0.35 ± 0.031^:^^ 
0.29 ± 0.03t[l:[]^ 

0.32 ± omto^l 

0.32 ± 0.031H! 
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Figure 22: The OPAL F^J/a data (circles) as a function of xbj for < >= 1.86 GeV^ (a) 
and <Q^>= 3.76 GeV^ (b). Also data from PLUTO (crosses) for <Q^ >= 2.4, 4.3 GeV^ 
and TPC/2j (squares) for < >=1.31, 2.83 GeV^ are shown. The curves show predictions 
of the GRV HO (dots), GRV LO (line), and SaSlD (dashed) parametrizations. The range 
of the XBj-bins of the OPAL results are marked at the tops of figures (from ^^). 



• OPAL, < ) = 1.86 GeV^ 

* PLUTO, < > = 2.4 GeV^ 
□ TPC/2y, < > = 1.31 GeV^ 



GRV-HO 
GRV-LO 
SaS-lD (LO) 




1 1 1 1 1 

• OPAL, < ) = 3.76 GeV^ 


o PLUTO, < > = 4.3 GeV^ 










; □ TPC/2Y, < > = 2.83 GeV^ 


(b) 



Comment: No correction for y^O was made. GRV HO is consistent with the low-x 
OPAL results in the lower bin, hut at higher it underestimates the low x OPAL 
data. GRV LO and SaSlD describe the unfolded results worse than GRV HO. Shapes 
of measured F2 are flat within the errors, but a small rise in the low x region is not 
excluded. 



•OPAL 97c (LEP 2) 

New data on F2 from LEP 2 at the CM energies 161-172 GeV were collected in 
1996 in two samples Q'^=Q-2Q GeV^ (0.004 < xbj < 0.76) and Q2=20-100 GeV^ 
(0.012 < XBj < 0.94). Also the distribution of the final hadronic energy flow was 
studied using HERWIG 5.9, PYTHIA 5.718 and F2GEN, see next section. The LEPl 
OPAL 97a data were used for comparison. For unfolding the Monte Carlo based 
on HERWIG program with the GRV parton parametrization was used. Two types of 
binning were performed, with < > = 9, 14.5, 30 and 59 GeV^, and with < > 
= 11 (combined samples 9 and 14.5 GeV^) and 41 GeV^ (30 and 59 GeV^). 

The unfolded F2 /a as a function of x Bj and are presented in figs. ^ and ^ 



and in tables 16, 17 
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Table 16: 



<Q'> 


XBj 


( fttnf, -\- susf 1 

\ Cf V KAjV • 1 tj KJ fJ V t 1 


9 


0.02-0.1 
0.100 -0.250 
0.250 -0.600 


0.33 ± 0.03l^:^[! 
0.29 ± OMt^om 
0.39 ± O.OSlgili 


14 5 


n 09 — 1 
0.100 - 0.250 

0.250 -0.600 


f) 07 _|_ f) AQ+U.itj 

0.42 ± 0.05lEj:?l 
0.39 ± OMt^o.li 


30 


0.050 -0.100 
100 — 350 
0.350 - 0.600 
0.600 - 0.800 


0.32 ± 0.04+^:^^ 
52 ± 05+n?o 
0.41 ± 0.09t^il 
0.46 ± 0.15l|iil 


59 


0.050 - 0.100 
0.100 -0.350 
0.350 - 0.600 
0.600 - 0.800 


0.37±0.06t^:^? 
0.44 ± 0.07i[]:[J? 
0.48 ± 0.09+g J^ 
0.51 ± 0.14l°;^8 



Table 17: 



<Q'> 


< F^/a > 


[GeV^] 


{stat. + syst.) 


9 


0.36 ± 0.05i^:^^ 


14.5 


0.41 ± QMtlTi 


30 


0.48 ± 0.05lS:o? 


59 


0.16 ± o.o6l[!:SI 
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Figure 23: The XBj-dependence of F^ja for different values of < >: 9, 30, 14.5, 59, 
11, and 41 GeV^. For comparison the GRV LO (sohd hne) and SaSlD (dot-dashed hne) 
parametrizations are shown; dotted line represents the hadronic component and the dashed 
one the (augmented) asymptotic solution (ASYM) (from ^] ). 
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Figure 24: The Q^- dependence of F2 /a averaged over 0.1< xbj <0.6, for the energy 91 
GeV (open circles) and for the energies 161, 172 GeV (full circles). Predictions of the QCD 
calculation are shown by the lines: solid (GRV LO), dotted (SaSlD), and double-dotted 
(HO, based on the GRV HO parametrization for light quarks); dashed line corresponds to 
the (augmented) asymptotic solution (ASYM) (from p^/j. 



A fit to the new F2 data at the energies 161-172 GeV and the previous OPAL 
97a set at 91 GeV for from 7.5 to 135 GeV^, averaged over the Xsj range of 0.1-0.6 
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(see fig. 



has the form: 



F:j{Q^)/a = (0.16 ± 0.05^°:}^) + (0.10 ± 0.02+°:°^) ln(QVGeV2). 
A special study of the dependence of F2 is performed for different xsj ranges, 



see fig. 
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Figure 25: The dependence of ja in three XBj bins. The data are compared to HO 
predictions based on the GRV HO parametrization (from J^/}. 



Comment: No correction for ^0 was made. Discrepancies are observed in the 
hadronic energy flow between the data and the HERWIG and PYTHIA simulations, 
especially at x^s <0.1. Accuracy of the data does not allow to see the expected different 
slope of F2 versus InQ"^ for different xbj ranges. 

•OPAL 2000 [ID (LEP 1, LEP 2) 

The low XBj data of F2 obtained at the CM energies 91 GeV, 183 and 189 GeV (in 
1993-5 and 1997-8) were collected in the partly overlapping ranges: 1.5-30 GeV^ 
at LEPl and 7-30 GeV^ at LEP2. Three samples: LEPl SW {< > = 1.9 and 3.7 
GeV2), LEPl FD {< > = 8.9 and 17.5 GeV^) and LEP2 SW {< > = 10.7 and 
17.8 GeV^) were studied. In the analysis new Monte Carlo programs and improved 
unfolding methods were introduced. The final state was studied using PHOJET 1.05, 
HERWIG 5.9, "HERWIG 5.9 + power law pt (dyn)" and F2GEN, with the GRV LO 
parametrization as an input function, see next section for details. 

"HERWIG 5.9 + power law pt (dyn)" and PHOJET 1.05, giving better description 
of OPAL data, were applied to obtain the F2 structure function. To reduce dependence 
of results on the implemented MC model, two-dimensional unfolding was performed 
with Xcor used as a first variable and E^^/Etot as a second one. To investigate the 
\ow-x Bj range, the unfolding of the XBj variable was done on a logarithmic scale. 

The values of F2 /a were obtained as an average of results of the unfolding with 
"HERWIG 5.9 + power law pt (dyn)" and PHOJET 1.05c, see table {x%j corre- 
sponds to the log centre of the bin for XBj below the charm threshold. The F2 / a was 
unfolded in the quoted bins and corrected to the values). 
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The unfolded F2 /a as a function of xsj for (5^=1-9 and 3.7 GeV^ is slightly lower 
than the previous OPAL results obtained using HERWIG 5.8d (OPAL 97b), see 
fig. ^ In fig. ^ the comparison with other measurements and parton parametriza- 
tions predictions for range 1.9-5 GeV^ is shown. 

The resuhs for the unfolded F:^ at < > = 8.9, 10.7, 17.5 and 17.8 GeV^ (similar 
to the previous ones, OPAL 97b; comparison not shown) are compared in fig. p5|with 
results of other experiments and predictions of various parametrizations, in similar 
range. 



Table 18: 



<Q'> 
[GeV^] 




{log centre 
of bin) 


F,ya 
{stat. + syst.) 


1.9 


0.0006- 0.0028 
0.0028- 0.0143 
0.0143 - 0.0724 
0.0724- 0.3679 


0.0012 
0.0063 
0.0322 
0.1124 


0.269 ± 0.027^^:^^^ 
0.177 ±0.009lJ]:[JlI 
0.179 ± 0.007+Hnfi 

0.227 ± o.oiol[5:[jl^ 


3.7 


0.0015- 0.0067 
0.0067- 0.0302 
0.0302- 0.1353 
0.1353- 0.6065 


0.0032 
0.0143 
0.0639 
0.1986 


0.269 ± 0.033^^:^1^ 
0.232 ± 0.013l[5:[ji 
0.259 ± 0.010lo:o?3 

0.296 ± o.oi4i[5:[!i 


8.9 


0.0111 - 0.0498 
0.0498- 0.2231 
0.2231 - 0.8187 


0.0235 
0.1054 
0.3331 


0.221 ± 0.0171^:^^^ 

0.308 ± o.oi4i[5:[;i^ 

0.379 ± 0.022l[5:[!J^ 


10.7 


0.0009- 0.0050 
0.0050- 0.0273 
0.0273- 0.1496 
0.1496- 0.8187 


0.0021 
0.0117 
0.0639 
0.3143 


0.362 ± 0.0451^:^;^^ 
0.263 ± 0.015lJ]:[!! 

0.275 ± o.oiilJ]:[!! 

0.351 ± 0.012in?^ 


17.5 


0.0235- 0.0821 
0.0821 - 0.2865 
0.2865- 0.9048 


0.0439 
0.1534 
0.3945 


0.273 ± OmSt'^omg 
0.375 ± 0.023lH?|i 
0.501 ± 0.027tlf^l 


17.8 


0.0015- 0.0074 
0.0074- 0.0369 
0.0369- 0.1827 
0.1827- 0.9048 


0.0033 
0.0166 
0.0821 
0.3483 


0.428 ± 0.061^^:^?? 
0.295 ± 0.019l[5:[ji 
0.336 ± 0.013l[j:[j^^ 
0.430 ± 0.013ini 
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X 

Figure 26: The XBj-dependence of F2 /a measured by OPAL for the LEPl SW sample with 
< Q2 > equal to a) 1.9 GeV^ and b) 3.7 GeV^. For comparison the previous OPAL data 
(OPAL 97h) and the present sample, both obtained using HERWIG 5.8d, are presented. 
The GRV LO predictions are also shown (from ^/}. 



~T ' ' — 

• OPAL, <QV=1.9 GeV^ 
« L3, <Q^>=1.9 GeV' 

PLUTO, <QS=2.4 GeV^ 
o TPC-2T, <qV=2.S GeV^ 



GRV LO Q =1.9 GeV 

- - - GRV LO Q^=2.8 GeV' 
SASID 0^=2.8 GeV^ 

- - - QPM Q-=2.8 GeV^ 




OPAL, <0V=3.7 GeV^ 
L3, <QS=5.0 GeV^ 
PLUTO, <QS=4.3 GeV^ 
TPC-2Y, <QV=5.1 GeV^ 



-GRV LO Q =3.7 GeV 

- GRV LO Q^=S.l GeV^ 
SASID Q^=5.1 GeV^ 

- WHITl Q'=S.1 GeV^ 

- QPM Q^=5.1 GeV^ 



(b) 




Figure 27: The XBj-dependence of F2 /a measured by OPAL for the LEPl SW sample with 
< Q2 > equai to a) 1.9 and b) 3.7 GeV^. For comparison the (L3 98a), PLUTO 1^ and 
TPC/2^ ^ data are shown together with GRV LO, SaSlD and QPM parton parametriza- 
tions predictions (from ^]). 



41 



OPAL, <Q >=10.7 GeV 
OPAL, <0V=S.9 GeV^ 
PLUTO, <OS=9.2 GeV^ 
DELPHI, <OV=12.0 GeV^ 
L3, <Q^>=10.8 GeV' 



-GRV LOQ^=8.9 GeV' 

- GRV LO Q^=12.0 GeV- 
SASID Q^=12.0 GeV^ 

- WHITl Q'=12.0 GeV^ 

- QPM Q^=12.0 GeV^ 




OPAL, <qV=17.S GeV^ 
OPAL, <Q->=17.5 GeV^ 
TOPAZ, <Q->=16.0 GeV^ 
L3, <QS=23.1 GeV^ 
ALEPH, <QS=20.7 GeV^ 



-GRV LO Q^=16.0 GeV- 

- GRV LO Q^=23.1 GeV' 
SASID Q^=23.1 GeV^ 

- WHITl Q'=23.1 GeV^ 

- QPM Q-=23.1 GeV- 




Figure 28: The XBj-dependence of F2 /a measured by OPAL for < > equal to a) 8.9 
(LEPl FD) and 10.7 GeV^ (LEP2 SW), and b) 17.5 (LEPl FD) and 17.8 (LEP2 SW) GeV^. 
For comparison the ALEPH 99a, DELPHI 96a, L3 99a, PLUTO §/ and TOPAZ 94 
data are shown with predictions of GRV LO, SaSlD, WHITl and QPM (from /g|/j. 



Comment: No correction for ^0 was made. The F2 measurement at "the lowest 
attainable xbj values", uses new MC models and two-dimensional unfolding, on a 
logarithmic scale. 

"The GRV LO and SaSlD parametrizations are generally consistent with the OPAL 
data in all the accessible Xbj and regions, with the exception of the measurement 
at the lowest scale, < > = 1.9 GeV^ , where GRV is too low." 
The precision of the measurement is insufficient to determine whether or not there is 
a rise in F2 at low x Bj ■ 

•AMY 1.5 95 (TRISTAN) 

The upgraded AMY 1.5 detector was used to perform a high measurement of the 
photon structure function F2 (with averaged = 73 and 390 GeV^). For previous 
measurement done with AMY 1.0 see 0. In the Monte Carlo analysis the FKP 
approach was used to describe light quarks (with pj. ^ 0.5 GeV), the QPM for c and 
b-quarks, and VMD contribution for soft hadronic processes were introduced. Results 



for unfolded F2 are presented in table |l^ and in fig. where predictions of various 



parton parametrizations as well as individual components are shown. 
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Table 19: 



< 02 > 




J. 2 / " 






{stat. + syst.) 


73 


0.25 


0.65 ±0.08 ±0.06 




0.50 


0.60 ±0.16 ±0.03 




0.75 


0.65 ±0.11 ±0.08 


390 


0.31 


0.94 ±0.23 ±0.10 




0.69 


0.82 ±0.16 ±0.11 




Figure 29: The F^jo- data from the AMY collaboration. Comparison with the following 
contributions to the photon structure: a) and b) FKP (u,d,s) with various pj^, QPM (c, 
b) with various charm masses, VMD. Comparison with parton parametrizations: c) and 
d) WHITl, LACl, DC; e) and f) GRV, AFG. Upper a),c),e) and lower b),d),f) figures 
correspond to averaged Q^=73 CeV^ and Q^=390 GeV^, respectively (from p^/). 



The values of F2 averaged over 0.3< xbj <0.8 are shown in table |20 



Table 20: 



<Q'> 
[GeV^] 


< F^'/a > 


73 
390 


0.63 ±0.07 
0.85 ±0.18 



Comment: The observed XBj-behaviour of F2 is consistent with the GRV parametriza- 
tion and with the FKP one. The fitted parameter is equal to 0.51±0.39 GeV if only 
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the AMY data are included, and 0.45±0.07 GeV if all available, at that time, data 
are taken. 



•AMY 1.5 97 O (TRISTAN) 

The measurement of F2 at < >=6.8 GeV^ was performed. In the Monte Carlo 
analysis the FKP approach was used to describe hght quarks (with ^ 0.5 GeV), the 
QPM for c-quark, and VMD contribution for soft hadronic processes were introduced. 
The results for F2 versus x^j, together with earlier data, are presented in fig. |3D|. A 
comparison of the data with various parametrizations is given in fig. pTl. 



K PLUTO Qi' = 4.3 GeV^ 
o PLUTO - 9.3 GeV" 
o TPC/27 q'^ = 5.1 GeV^ 



□ TOPAZ q'^ = 5.1 GeV^ 
y- OPAL q'^ = 5.9 GeV^ 
• AMY = 6.8 GeV^ 




0.0 

0.2 0.4 0.6 0.6 1 

X 

Figure 30: The F2 /a data from the AMY coUaboration in comparison with the results of 
other experiments at ~ 4-9 GeV^ (from /p^/). 
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Figure 31: The F^ja. as a function of xbj- Left: Comparison with the SaSlD, SaS2D 
and GRV LO parton parametrization predictions. Also shown are the c-quark contribu- 
tions from the SaS and GRV LO parametrizations. Right: Comparison with the WHITl, 
WHIT4, WHIT6 and LACl predictions. The c-quark contribution is shown for the WHITl 
parametrization (from |5^/j. 



Comment: "The xbj behaviour of the measured F2 is consistent with the QCD-based 
predictions such as SaS, GRV LO and WHIT models, but inconsistent with the LACl 
prediction for x sj values around 0. 01. " 

•TOPAZ 94 g (TRISTAN) 

The photon structure function F2 has been measured for averaged values from 5.1 
to 80, and also 338 GeV^. The Monte Carlo with FKP contribution included for light 
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quarks (with = 0.1 - 1 GeV), the QPM one for charm quark and VMD contribution 
for a soft particle production. The one and two-jet events ("hadronic" events with a 
stuck quark and the pointhke ones based on 7*7 qq, respectively) in this sample 
are studied, for the first time in DlSg-y process (see also sec. and p^ ). 



Table 21: 



<Q'> 
[GeV^] 


g2 
[GeV^] 


< XBj > 


XBj 


F,ya 
{stat. + syst.) 


5.1 


3-10 


0.043 
0.138 


0.010-0.076 
0.076 - 0.20 


0.33 ±0.02 ±0.05 
0.29 ±0.03 ±0.03 


16 


10 - 30 


0.085 
0.24 

0.555 


0.02 -0.15 
0.15 -0.33 
0.33 -0.78 


0.60 ±0.08 ±0.06 
0.56 ±0.09 ±0.04 
0.46 ±0.15 ±0.06 


80 


45 - 130 


0.19 

0.455 

0.785 


0.06 -0.32 
0.32 -0.59 
0.59 -0.98 


0.68 ±0.26 ±0.05 
0.83 ±0.22 ±0.05 
0.53 ±0.21 ±0.05 
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Figure 32: The F2 /a versus xsj for < >=5.1, 16, 80 GeV^ and a comparison with the 
prediction of FKP(u,d,s), QPM (c), VMD contributions with various pj^ parameters (upper 
set) and with parton parametrizations: DG, LACl, GRV (lower set) (from Jsll/j. 



The results for < > = 5.1, 16 and 80 GeV^ are presented in table pil and 
in fig. |3^ where the individual components as well as predictions of various parton 
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parametrizations are shown. To study the dependence of the averaged values 
of F^ja. in the x^j-range from 0.3 to 0.8 were extracted at < Q'^ >= 16, 80 and 338 
GeV^ (see table |2^; the numbers in parentheses are results for the light quarks alone). 



Table 22: 


<g2> 




[GeV^] 




16 


0.47 ±0.08 (0.38 ±0.08) 


80 


0.70 ±0.15 (0.49 ±0.15) 


338 


1.07 ±0.37 (0.72 ±0.37) 



Comment: The final hadronic state described by QPM(c) + VMD+FKP(u,d,s) with 
Pj'=0. 1,0.5 and 1 GeV, and also by the GRV, DC and LACl parton parametrizations, 
was studied. 



For an overall comparison the collective figures of F2 /a versus xbj (fig- 0) and 
F2 /a versus (fig. 0), containing also earlier data not discussed here, are presented 
based on |5^. In these figures the comparison with theoretical predictions for xbj 



dependence of F2 based on the SaSlD (LO) and GRV NLO parton parametrizations 
(in fig. 1^ in addition the prediction of the asymptotic solution) are shown. For 
comparison with theoretical predictions of the evolution for different xbj ranges 
see also figs. 0, |10|d, |1|,|11, ^ ^ 0, gl|d, H |5[ 

The effective parton density, as measured at HERA collider by the HI group ||5^ 
in the jet production from resolved photons, is compared with the F2 data in fig. ^ 
as well. 
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Figure 33: The photon structure function F2 /a as a function of x bj in bins of compared 
to the GRV NLO (sohd hne) and SaSlD (LO) (dashed hne) parametrizations of parton 
distributions in the photon (from J5^/ ). 
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Status: 3/2000 (home.cern.ch/s/soldner/www/) 
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Figure 34: The photon structure function F2 /a (data from various experiments) as a 
function of (from j53/j. 
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Figure 35: The scahng violation of the effective parton densities measured at HERA com- 
pared with the dependence of F2 (averaged over the range 0.3< xbj <0.8) (from j5j/j. 
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2.3 Hadronic final states in the DlSe^ experiments 

Although a detailed analysis of hadronic final states, including those involving large 
Pt particles/jets, has not been the main aim of DlSg^ experiments, it has proved 
crucial in extracting the unfolded F2 as mentioned in the previous section (sec. p.2.2| ). 
The obtained experimental results, and problems that appear while describing the 
hadronic energy fiow, transverse energy out of event plane, pseudorapidity and other 
distributions within the existing Monte Carlo models, deserve close attention and a 
separate treatment - this section is devoted to this subject. 

The F2 corresponds to the cross section for hadron production in 7*7 collisions, 
eqs. (^, H), and obviously the bulk of the contribution is due to the soft processes, i.e. 
production of particles with a relatively small transverse momentum p^- One should be 
aware that in the hadronic final state in DlSe-y experiments two types of large scale may 
appear: andp^ [^. The bulk of the data, as it was already mentioned, corresponds 
to events with not very large pr] if for these events the relation ^ p^ ^ holds, 
the interpretation in terms of the photon interaction between one direct (7*) and one 
resolved real (7) photon may be introduced. Then it is not clear what scale should be 
used in the parton density for a real photon, f^^x, Q"^): = p^, Q^, or 7 Moreover 
the processes corresponding to -C p'^ should be treated as being resolved from the 
point of view of both the real and the virtual photon. 

Some of the problems that will appear here are common with those presented in sec. 



for the jet production in 7*7 collision and in sec. 2.4 where dedicated measurements 



of the large pt jet production in real 77 and 7p processes are presented. 

2.3.1 Modelling of the hadronic final state in 7*7 collisions 

In the DlSe-y (single-tag) events the whole hadronic final state is not seen in the 
detector. A correct unfolding method is necessary to reconstruct the true kinematical 
variables (for definition see beginning of sec. p.2.2|) . At this stage models to generate 



separate samples for the different components of the photon are needed, as we have 
already mentioned in sec. |2.2.2| . The modelling is performed using different MC 



programs and it is still a large source of systematic errors in F2. 

We start with a short description of the FKP approach used in the early analyses 
which bases on the QPM, (soft) VDM and RPC contributions. The QPM contribu- 
tion is generated using the e~^e~ — * e^e~qq matrix element for massive quarks, for 
an arbitrary value of Pt of produced quarks. The final state transverse momentum 
distribution is here da/dp"^ ~ p^'^. The QPM was applied mainly to the c and b quark 
production. The (soft) VMD (called sometimes hadronic) contribution, in which the 
initial real photon is treated as the vector meson(s), say p, with the typical soft inter- 
action with a photon probe, leads to the production of the final quarks (then hadrons) 
with the da/dp\ ~ e~^^T distribution (where 6 ~ 3 - 5 GeV~^). This contribution 
should be negligible for the production of particles with large px where, on the other 



hand, other contributions should be included. In the FKP approach |^ the cutoff pa- 
rameter p9p is built in to separate the mentioned above VMD component of the initial 
photon target and the resolved photon (QCD) contribution, which develops from the 
pointlike coupling of the 7 to the qq pair. This parameter, p^., corresponds therefore 
to the intrinsic transverse momentum of the partons in the photon. The additional 
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cutoff p™" is introduced in the analysis, to define hard (2 2) subprocesses where the 
pQCD can be used. In the description of the particle production with p-r > p^*" the 
resolved photon contribution is introduced (see also sec. 2.4). In the FKP approach, 
where the (soft) VMD contribution was introduced without the dependence, both 
cutoff parameters are of the same order. 

In the DELPHI analysis, the TWOGAM generator ("fixed" in 1993) was used, with 
the QPM, VMD for soft production and RPC contributions included within the FKP 
approach |^, This generator treats exactly the kinematics of scattered electron 
and positron, uses exact expression for the two photon luminosity function, and also 
in the (single and double) resolved process the kinematics of the partonic system is 
exact for any photon virtuality. The TWOGAM generator was used also in the last 
L3 measurements. 

In many OPAL analysis the generator F2GEN (developed from TWOGEN used 
in older analyses) was used, which can generate events according to "any selected 
formula for 0]a. In this generator one can easily put a specific assumption about 
the final state. For example different angular distributions may be assumed in the 
7*7 CM system for the produced pair of quarks: "pointlike", i.e. as for lepton pair 
from two real photons, "peripheral" with an exponential distribution of the transverse 
momentum, or the "perimiss" combination . 

In recent analyses general purpose Monte Carlo generators are being used in anal- 
ysis of hadronic final state. HERWIG and PYTHIA (the implementation of the SaS 
approach) apply chosen parton parametrizations. In the PHOJET generator, where 
hard interaction is described using QCD, the description of soft interaction bases on 
the Regge theory. The very recent ALEPH, L3, OPAL and TRISTAN F2 data have 
been obtained using the modified HERWIG (PYTHIA) program, where the additional 
power distribution of the transverse momentum of partons in the photon, of the form 
dPt/iPt +Pto)^ is included, see Q and ALEPH 99a, ALEPH 99b,conf or OPAL 
2000. 

As far as the global hadronic variables like Wyis, Q^, Mjj (invariant mass of two 
jets) and other distributions are concerned, there is a fairly good description of the 
data by existing Monte Carlo generators. The problems arose for the transverse energy 
out of the tag plane (the plane defined by the initial and tagged electrons) or for the 
energy flow per event as a function of pseudorapidity, r] = — ln(tan(6'/2)) 0. The 
discrepancies were very pronounced especially for small xbj- 

The first observation of the disagreement between the data and Monte Carlo models 



was made by the OPAL collaboration |^. "The serious discrepancies between the 



data and any of the available Monte Carlo models are seen both within the central 
region of the detector {\ri\ < 2.3), where the energy flow is well measured, and in the 
forward region, where the energy can only be sampled." It is clear that the unfolding 
of F2 "will have large errors as long as the energy flow from different models remains 
in clear disagreement with the energy flow in the data, in particular in the region of 
Xyis <0.1 and <30 GeV^" (from [0). A similar effect has been seen by now by 
other groups, see below. 

The improved description of the mentioned data is obtained in the "HERWIG 

^^The polar angle is measured from the z — axis, and the tagged electron(positron) is by definition 
at the negative rapidity. 



50 



(PYTHIA) + power lawpt" generators. The PHOJET 1.05 and TWOGAM generators 
can also satisfactorily describe the data (DELPHI, L3). Note however that the values 
of pj?'" fitted in these analyses are high (2.5 - 3.5 GeV). Another approach which 
usually improves the agreement with the data is the inclusion of multiple interaction. 
The effect is similar in all these approaches and consists in increasing pj- in processes 
involving photons. 

The recent progress in describing the hadronic final state in single-tag events at 
LEP by using the two-dimensional unfolding led to a reduction of the Monte Carlo 
modelling error (OPAL 2000). As follows from first combined results of the newly 
formed LEP-wide Two Photon Physics group, which used the "HERWIG 5.9 + 
power law pt" and PHOJET 1.05, a good agreement is found between experiments in 
the central regions. "Unresolved problems remain in the small angle regions where 
the results of the different experiments showed greater discrepancies than could be 



accounted for by statistical fluctuations..." [0. It was also found that to describe 
low data an improved, modified model "HERWIG 5.9 -|- power law pt (dyn)", with 
a dynamical limit p'^j^ax — + Pt ^o be introduced. 

For the general discussion and data - Monte Carlo comparison in the context of 
modelling the hadronic final state see |^, see also 0, |n| and [jl2|, where the jet 
production in the photon-induced processes is studied. 

2.3.2 Measurements of the hadronic final state 

The analyses of the hadronic final state accompanying the measurements of dis- 
cussed in sec. p.2.2| , as well as the results from independent analyses of large px hadron 
production in single tagged events are presented below. Variables used in the unfold- 
ing procedure are defined and discussed in sec. ^^2^21 . Note that in some analyses the 



jet topology of the final state was studied (see e.g. TOPAZ 94, OPAL 97d,conf). 

For the two photons involved, we introduce the following notation for their squared 
(positive) virtualities: (= Q^) and P2 (= -P^), with the relation to variables used 
in the discussed DlSe^ scattering in parentheses, i.e. with the inequality > P2 
understood. 
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DATA 



•ALEPH 97a,conf |]5§ (LEP 1) 

Single tagged events were collected in years 1992-1994 at the average Q^=14.2 GeV^. 
A dedicated study of the hadronic final state assuming the QPM (for u, d, s and c 
quarks) + VMD (for the target photon) and using the HERWIG 5.9 generator, with 
the GRV parametrizations, was performed. Also the modified HERWIG generator, 
"HERWIG 5.9 + power law pt" , with ptQ = 0.66 GeV was used. A cone jet-finding 
algorithm was used in the analysis. 



Results for energy fiow are presented in figs. ^T\. The discrepancy is observed 
in the energy fiow in the positive pseudorapidity region. The approach "HERWIG 5.9 
-|- power law pt" leads everywhere to a better description of the data. 
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Figure 36: The energy (out of the plane of the tag and the beam) distribution. Comparison 
with different Monte Carlo models: HERWIG 5.9 (shaded), "HERWIG 5.9 + power law pt" 
(solid line) and QPM + VMD (dashed line) (from §^). 



-\- ALEPH data 
HERWIG 5.9 

HERWIG 5.9 + power law pt^ 

---- QPM + VMD 




Figure 37: The energy Bow versus pseudorapidity of the final hadrons. Comparison with 
different Monte Carlo models as in fig. ^ (from ^dj]). 



In order to pin down the source of the observed discrepancy the analysis of the num- 
ber of cone jets was performed in the hadronic final state (see fig.p^ for results). The 
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standard HERWIG model, having additional production mechanisms with a resolved 
photon, should in principle give larger number of two-jet events than the QPM+VMD 
model. On the contrary, it leads to a lower number. 

To check the presence of the single and double resolved photon events, the 
and Xtag distributions were measured (not shown). Note, that double resolved photon 
events contain the resolved virtual photon contribution (see also sec. ^). Also the 
energy not assigned to jets in two-jet events, -v^^ag studied (not shown). The 

"HERWIG + power law pt^ approach was used successfully to describe these data. 
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Figure 38: The number of events classified in four bins referring to various numbers of jets 
in an event. Comparison with different Monte Carlo models as in fig. (from 
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Comment: "The ALEPH data confirms recent problems reported by OPAL in modelling 
the hadronic final state. " 

" 'HERWIG + power law pt ' is better at modelling the region of large Et^out " o-nd the 
peak in the energy flow at positive rapidity. 

The "double resolved processes are not present to an large extent in the data at this 



•ALEPH 99a J3§ (LEP 1) 

Analysis of the hadronic final state during the F2 measurement is reported. Monte 
Carlo models: QPM (for all quarks) + VMD (for a target photon) and "HERWIG 5.9 
+ power law pt \ with GRV LO parametrization, were studied. Results are presented 
in figs. p9| and HO. 
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Figure 39: The data in which the scattered electron is detected (tagged) in luminosity 
calorimeter LCAL, compared with the QPM + VMD model (solid histogram) and the "HER- 
WIG 5.9 + power law pt" model (dashed histogram) (from j3^/j. 
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Figure 40: The data in which the scattered electron is detected (tagged) in electromag- 
netic calorimeter ECAL, compared with the QPM + VMD model (solid histogram ) and the 
"HERWIG 5.9 + power law pt" model (dashed histogram) (from j^). 
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Comment: The HERWIG program with modified pt distribution and the QPM + VMD 
model give "an acceptable description of the data". 



•ALEPH 99b,conf fSJ (LEP 2) 

A new unfolding method was used in the measurement of F2 at ^/s=l'&?) GeV based 
on 1997 data (see also previous section). Two samples of single-tag events with the 
tag defined as a cluster in one of the small angle luminosity calorimeters SICAL or 
LCAL, for Q2 ranges: 7-24 GeV^ and 17-200 GeV^, were studied. 

The new unfolding procedure, based on the principle of maximum entropy, was 
used to correct the distributions for finite resolution and acceptance. The unfolding 
was done for XBj and in additional for a second kinematical variable, i^iy, defined as 
the total energy of particles produced at an angle with respect to the beam direction 
smaller than 17°. In fig. ^ results of unfolding of the x variable for different values 
of En obtained using "HERWIG 5.9 -|- power law pt' (SaSlD) are presented (LCAL 
data). One can see that the small Eu gives better result. 




Figure 41: The relationship between reconstructed and true values for xsj for four bins of 
En (from /H/j. 



Although the two-variable unfolding reduces the dependence on the Monte Carlo, 
some model dependence still remains. This point was studied by comparing several 
kinematic distributions obtained using "HERWIG 5.9 + power law pt" with the GRV 
LO and SaSlD parton parametrizations and in addition using the PHOJET model. 
The SICAL-tag data (not shown) are described reasonably by the used models, while 
for LCAL data (fig. H^) "PHOJET is in poor agreement with the data". 
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Figure 42: Comparison of data and simulations for the LCAL tagged data. Predictions of 
the "HERWIG 5.9 + power law pt" model with two parton parametrizations: GRV (solid 
line), SaSlD (dotted line), and of PHOJET (dashed line) are shown (from 



Comment: This unfolding method "leads to smaller statistical errors and a reduced 
model dependence compared to one- dimensional procedure". 

•DELPHI 95 (LEP 1) 

First evidence of the hard scattering process in the single tagged events in the data 
from the 1991-1992 run is reported. The values of Et of observed jets were larger than 
1.5 GeV, while the magnitude of the mass of virtual photon was equal to < >~ 
0.06 GeV^. This corresponds to the standard resolved (almost real) photon process, 
from the point of view of both photons. 

Comment: This is not a typical DIS experiment since the photon probe is almost 
real. The analysis is not typical, either, for the standard large px jet study. "The 
data are consistent with the predictions for quark and gluon density functions in the 
GS parametrization. The sum of the contributions: VMD + QPM (Quark Parton 
Model) + [QCD - RFC] is needed in order to describe the data; the DO and LACS 
parametrizations do not adequately describe the data. " 

•DELPHI 96b,conf [|§ (LEP 1) 

The results for the averaged Q^=13 and 106 GeV^ were considered and compared with 
the TWOGAM generator, for results see sec. ^ The QPM, GVMD contributions 
and RFC with the GS2 parton parametrization were included. 
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Comment: It was shown that the resolved photon contributions are needed to obtain 
the description of the data. 



•DELPHI 97a,conf (LEPl, LEP2) 

Analysis of the hadronic final state in the F2 measurement (see sec. |2.2.2| ) was 
performed for energies around the mass, and for 161-172 GeV (1996 run). The 
< >= 6.5, 13, 22 and 17, 34, 63 GeV^ were studied, respectively. The hadronic 
final state topology with events containing jets was studied using the TWOGAM gen- 



erator. For various distributions obtained for LEP 2 data, see fig. In fig. ^ the 
event energy fiow as a function of the pseudorapidity is shown for LEPl and LEP2 
data. 



Comment: "All variables were found to be in good agreement with TWOGAM predic- 
tion. " 

Note that the px range of jets may be of the order of < > for small samples. 

LEPl - <Q^>=34 {GeW/cf 




Etag/Ebeam 



•S 











6 - 












4 - 






















,n 


E- 




P; (GeV ) 




I^'iGeV ) 



Figure 43: Comparison of data and the Monte Carlo prediction for < > = 34 GeV^ 
(LEP2) for the variables : (shown first left, then right in each row) Energy of tagged particle; 
The invariant mass; The number of jets; The jet angle with respect to the tagged particle; 
Pt for jets in the same hemisphere as the tagged particle; px for jets in the opposite hemi- 
sphere. Curves show the Monte Carlo predictions: GVMD-hQPM+RPC (GS2) (solid line), 
GVMD+QPM (dots), GVMD (dashes) (from §^). 
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Figure 44: The event energy How as a function of the pseudorapidity for < > = 13 GeV^ 
(LEPl) and < Q"^ > = 34 GeV^ (LEP2). The notations are as in fig. |3| (from 



•DELPHI 98,conf §^ ||] (LEP 2) 

The hadronic final state analysis accompanied the F2 data taking in the period 1996-98 
for Q2 between 10 and 1000 GeV^ (for the e+e" energy of 163-188 GeV), see sec. |]2]2. 
The three-component description of DIS events was applied: the TWOGAM generator 
was used for the QPM (all quarks), the soft hadronic GVMD and RFC (with the GS2 
parton densities in single and double resolved photon processes) contributions. The 
results are presented in figs. E5| and Good description of the data was obtained. 
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Figure 45: The cross section for averaged Q^=42 GeV^ as a function of a) The number of 
jets, b) the jet angle with respect to the tagged particle and of the transverse momentum of 
the jet for c) tagged and d) untagged events. Curves show the Monte Carlo predictions with 
QPM + GVMD + RFC (solid line), QPM + GVMD (dotted line) and GVMD (dashed line) 
(from 
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Figure 46: The cross section for averaged Q^=400 GeV^ as a function of a) number of jets, 
b) the jet angle with respect to the tagged particle, and of the transverse momentum of the 
jet for c) tagged and d) untagged events. The notations are as in fig. |4^ (from /|3^j. 
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•L3 98a g (LEP 1) 

The study of the hadronic final state was performed while the F2 data were collected 
in years 1991-95 for 1.2< <9.0 GeV^ (see sec. ^^T^). Two bins: < >=i.9 GeV^ 
in the xbj range 0.002< xbj < 0.1 and < >=5 GeV^ in the xbj range 0.005< xbj < 
0.2 were analysed. A new method for reconstructing kinematical variables, based on 
imposing a transverse momentum conservation, was used in the analysis. 

The invariant mass of the hadronic system obtained using PHOJET 1.05 - W^^, 



as well as Wyis and Wrec (from the new method) are compared in fig.^ In further 
analysis of xbj, Pt distributions and of the energy fiow three Monte Carlo generators : 
PHOJET 1.05, HERWIG 5.9 and TWOGAM have been used. The results are shown 



in figs. ^ and [49 



Two data sets for F2 based on the PHOJET and TWOGAM analyses are presented 
(they differ up to 28 %), see previous section. 

Comment: "A significant improvement is seen in the Wrec variable which uses the 
constraint of transverse momentum conservation. 

HERWIG disagrees with data in the small xbj region." The distribution of the trans- 
verse momentum of charged particles and energy flow were reproduced by PHOJET 
and TWOGAM, the spectrum of HERWIG was found to be too soft. 
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Figure 47: Monte Carlo (PHOJET) comparison of the generated two-photon mass W-y-y 
distribution with the distributions W^s and Wrec (from pOJ). 
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Figure 48: The px distribution for charged hadrons: data and Monte Carlo predictions 
(from g^). 
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Figure 49: Energy flow for two Q bins and different x regions as a function of the pseudo- 
rapidity, with the polar angle defined with respect to the tagged electron. Comparison with 
three Monte Carlo generators is shown (from j^]). 



•L3 2000 |42| (LEP 1) 

The analysis of the hadronic final state was performed in the measurement of the 
structure function for the real photon and virtual photon at ^/s= 89-92 GeV (see 
also sees. p.2.2| and |33). The large events , between 40 to 500 GeV^ with < 
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>=120 GeV^, were studied. The mean virtuality of the target photon is estimated 
to be < >=0.014 GeV^ for single-tag events and < >=3.7 GeV^ for double-tag 
events. "A new feature in the present work is a fit to two-photon kinematics imposed 
on each event." 

Single-tag events were analysed using the JAMVG generator modelling the QPM 
contribution with Nf=4, the PHOJET 1.05c with a cutoff p™"=2.5 GeV (and charm 
generated using JAMVG), and TWOGAM with |}™"=3.5 GeV generating three pro- 
cesses: QPM (all quarks), the soft hadronic VMD and QCD resolved photon contribu- 
tion (RPC). Kinematical variables were reconstructed in a fitting procedure applying 
constrains from four-momentum conservation. The correlation between generated and 



fitted values for the invariant mass and xbj is shown in fig. pO 



The Monte Carlo predictions and the data are shown in figs. ^ and |52| for single- 
tag events. The JAMVG (QPM) is in agreement with the data for xbj above 0.5. 
PHOJET agrees with the data at low xsj- Similar distributions were studied for 
double tagged events, see fig. p3, where some disagreement is seen. 
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Figure 50: a) Correlation between values of W generated and obtained in the unfolding 
after the fit, b) the same for xbj- Based on the JAMVG Monte Carlo (from J^/). 
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Figure 51: Distributions for the single-tag data. Comparison with PHOJET, TWOGAM, 
JAMVG predictions and with the background estimation is shown. The Monte Carlo distri- 
butions are normalized to the same number of events as the data (from p^/j. 
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Figure 52: Results from the single-tag events, a) xju b) yjn c) and d) Wfu distributions. 
Comparison with predictions of the PHOJET, TWOGAM and JAMVG models and with an 
estimation of the background is shown (from J4^j. 
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Figure 53: Results from double-tag events, 
distribution (from 



(GeV^ 

a), b) and c) as in fig. |5^, and d) the 



Comment: "...it is not clear if (this) successful analysis (i.e. unfolding) can be applied 
also to and Xsj regions not dominated by the direct process" (from j^). "PHOJET 
and TWOGAM expectations exceed the data at low xbj values" for double-tag events. 

•OPAL 94 [g (LEP 1) 

Data for F2 were collected in 1990-92 for averaged Q^=5.9 and 14.7 GeV^. Early 
analysis of the final hadronic state (the TWOGEN generator with the contribution 
based on QPM, the soft hadronic VMD process and on the FKP parametrization for 
the QCD contribution) was performed. The parameter p^ = 0.27 ± 0.10 GeV was 
determined. 



•OPAL 97a (LEP 1) 

The hadronic final state was analysed in details in the measurement of F2 with one 
photon highly virtual (Q^ between 6-30 GeV^ - low sample, and between 60-400 
GeV^ - high sample), the other being almost real. For generating the hadronic 
final state the Monte Carlo programs HERWIG 5.8d, PYTHIA, and for comparison 
F2GEN were used with the GRV and SaSlD parton parametrizations. 

The energy (i?t,out) transverse to the plane defined by the beam axis and the tag 
direction, and other quantities for the low sample are presented in fig. The 
discrepancy found for the Eyis and Et^out distributions in this sample (figs. |5^ b and 
d, respectively) is absent in the high events (not presented). 

The hadronic energy Et,out distributions in the three xbj bins are shown in fig. |55. 
The failure of the models in the low region is most visible at low xbj- 

To establish a source of the discrepancy the energy fiow per event in the low 
sample was also studied as a function of the pseudorapidity (see fig. EB) . The distribu- 
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tion of pseudorapidity r] for the low sample, corrected for the experimental effects, 
is shown in fig. |5^. 
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Figure 54: Comparison of the data event quantities in the low-Q^ sample with HERWIG, 
PYTHIA and F2GEN (PL) Monte Carlo simulations, a) the distribution of the visible 
invariant mass; b) the total visible energy of the event; c) the transverse momentum of the 
event in the tag plane; d) the energy out of the tag plane (from B^j. 
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Figure 56: The hadronic energy Bow per event as a function of the pseudorapidity rj for the 
data and various Monte Carlo simulations, in various ranges of x^is and Otag for the low 
sample (from 
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Figure 57: The measured energy Row in the low-Q'^ sample corrected for the detector 
inefhciencies, as a function of the pseudorapidity rj, compared to the values generated by 
various Monte Carlo models (from p^/j. 



Comment: None of the generators represents the final state accurately; Eyis,Et^out - 
distributions as well as the hadronic flow per event show a clear discrepancy. 
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"The failure of the models in the low region is most marked at low xbj- " 
The differences between the Monte Carlo models and the data in the low region are 
more pronounced for the energy flow per event as a function of the pseudorapidity and 
the azimuthal angle. 

"Particular attention will need to be given to the angular distribution of partons in 7*7 
system". 

The relation between p\ and in the low sample may indicate a need to take 
into account the structure of the virtual photon. 

"The differences between the data and the models contribute significantly to the sys- 
tematic errors on F2 ". 



•OPAL 97b [gD (LEP 1) 

The measurement of F2 and the modelhng of the 7*7 fragmentation into hadrons at 
low region (1.1 to 6.6 GeV^) and very small xbj - bins from 0.0025 to 0.2 is re- 
ported. The hadronic energy flow as a function of rj for three xbj regions is plotted in 

fig. m 

Differences between the data and Monte Carlo models (HERWIG, PYTHIA, F2GEN) 
in the energy flow distributions versus the pseudorapidity (fig. ^8]) and in the summed 
energy transverse to the tag plane (not shown) are found for < 0.05. 
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Figure 58: The measured hadronic energy Bow as a function of the pseudorapidity. Com- 
parison of data with various Monte Carlo predictions using the CRV LO and SaSlD parton 
parametrizations for three bins in Xms ctnd two bins in 6 tag (from p^/j. 



Comment: Discrepancy between data and Monte Carlo models in small-x region is 
found. 
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•OPAL 97c [|^ (LEP 2) 

The data on at the CM energies 161-172 GeV were collected in 1996 in two samples: 
Q2=6-20 GeV2 (0.02 < xbj < 0.6) and Q2=20-100 GeV^ (0.05 < xbj < 0.8). The 
analysis is similar to that of OPAL 97a. The final hadronic energy fiow was studied 
using the HERWIG, PYTHIA and F2GEN models. The LEP 1 data OPAL 97a, 
97b were used for comparison. For unfolding the HERWIG program with the GRV 
parton parametrization was used. Two types of binning were performed, with < > 
= 9, 14.5, 30 and 59 GeV^, and with < > = 11 and 41 GeV^. 

Comment: "A slightly better agreement between the predictions for the hadronic energy 
flow of the various models and the data in the region 0.1< x^s <0.6 is found for the 
data taken at ^/s^=161-172 GeV than for the data collected at ^/s^=91 GeV. At 
Xvis < 0.1 significant differences persist, as the data prefer a more pointlike hadronic 
energy flow than assumed in either HERWIG or PYTHIA. " 



•OPAL 97d,conf [|5|] (LEP 1) 

An analysis of the hadronic final state was done, in which the discrepancies between 
the data and predictions, reported in OPAL 97a |^ , were examined from the point 
of view of the number of produced jets. The data for ^ 6 — 30 GeV^ taken in 
years 1994-95 were compared to the results of the HERWIG, PYTHIA and F2GEN 
generators. 



The numbers of events for different numbers of jets, divided by the sum of all 
events, are presented in fig. BU. 
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Figure 59: The fraction of events with 0-3 jets. The points show the data with the back- 
ground subtracted, with statistical errors; histograms were obtained with the HERWIG, 
PYTHIA and F2GEN generators (from 



The results of further studies of the energy fiow versus the pseudorapidity rj for 



events with different number of jets are shown in fig. 60 
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1 jet events 2 jet events 

Figure 60: The hadronic energy Eow per event in bins of the pseudorapidity rj. Plots a), b), 
c) and d) show the average energy Eow per event for summed events, and for events with 0, 
1 and 2 jets, respectively. The samples are represented as in fig. (from ^^/}. 



The number of events versus the energy transverse to the tag plane, obtained by 
the different generators and observed in the experiment, are compared in the fig. ^ 
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Figure 61: The number of a) one-jet and b) two-jet events as a function of the jet transverse 
energy, Ex- The HERWIG, PYTHIA and F2GEN events are normalized to the a) one-jet 
and b) two-jet data events (from J5^j. 
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Comment: "There is a marked difference between the data and the Monte Carlo samples 
in the number of events with 2 jets". "While the F2GEN sample, generated using the 
pointlike approximation, overestimates the 2 jet rate by a factor of 2.4, HERWIG and 
PYTHIA are too low by factors of 2.4 and 3.6 respectively." 

Disagreement between the data and each of the generators for the hadronic energy flow 
versus the rapidity is seen for all types of events: 0, 1 and 2-jet. 
"All of the Monte Carlo samples model well the Et^out distributions for events with 1 
jet, but PYTHIA underestimates the Et^out of the events with 2 jets. " 



•OPAL 2000 g (LEP 1, LEP 2) 

Hadronic final state analysis accompanying measurement of F2 was performed at the 
CM energies 91, 183 and 189 GeV (data from 1993-5 and 1997-8). Two detectors were 
used: the small-angle detector SW allowing to access the lowest xbj region, and the 
forward detector (FD) providing the LEPl data in the same range as the LEP2 
data. Three samples were studied: LEPl SW with < >= 1.9, 3.7 GeV^ and 
0.0006 < XBj < 0.6065, LEPl FD with < >= 8.9, 17.5 GeV^ and 0.0111 < xbj < 
0.9048 and LEP2 SW with < >= 10.7, 17.8 GeV^ and 0.0009 < xbj < 0.9048. 
In the analyses new Monte Carlo programs, and improved unfolding method - the 
two- variable unfolding, were used. The unfolding of xbj was done on a logarithmic 
scale in order to study the low-x^j region in more detail. 




Figure 62: The comparison of the LEPl SW data with MC predictions ("HERWIG 5.9 + 
power law pt (dyn)" and PHOIET 1.05) for a) Etag/Eb, where Etag is the energy of tagged 
electron and E^ is the energy of the beam, b) electron scattering angle 9tag, cj and d) 
Ea/Eh, where Ea is the highest energy cluster opposite to the tag (from p^/j. 

The final hadronic state was studied using HERWIG 5.9, "HERWIG 5.9 + power 
law Pt (dyn)", F2GEN (all with GRV LO parametrization as an input function) and 
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PHOJET 1.05 programs. 

In fig. ^ the LEPl SW data are compared to the "HERWIG 5.9 + power law pt 
(dyn)" and PHOJET 1.05 distributions for variables related to the scattered electron. 
Both the generators underestimate the data (similar comparisons for LEPl FD and 
LEP2 SW samples give better agreement - not shown). 

For the variables related to the hadronic final state there are "large discrepancies" 
in comparison of data and Monte Carlo, what can be seen in fig. ^ where various 
measured distributions for the LEPl SW sample are shown together with "HERWIG 
5.9 + power law pt (dyn)" and PHOJET 1.05 simulations. 




Figure 63: The comparison of the data (LEPl SW) with the MC predictions ("HERWIG 
5.9 + power law pt (dyn)" and PHOJET 1.05) for a) invariant mass Wyis, b) Inxms, c) 
number of tracks in the event, and d) E^^ - the energy out of the tag plane. The dominant 
background contributions are also shown (from J5^/). 

The hadronic energy flow as a function of pseudorapidity in various MC models 
for LEPl and LEP2 data is presented in fig. The HERWIG 5.9 generally under- 
estimates data, while F2GEN overestimates data in the rapidity range from -1 to 2.5. 
The PHOJET 1.05 gives the best description of the data for LEPl SW and LEPl FD 
samples. For LEP2 SW sample "HERWIG 5.9 + power law pt (dyn)" seems to be in 
better agreement with data than PHOJET 1.05. 

The correlation between the Monte Carlo generated invariant mass and its visible, 
reconstructed and corrected values was studied (not shown). 
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11 



Figure 64: The energy How per event as a function of the pseudorapidity is shown for the a) 
LEPl SW, b) LEPl FD and c) LEP2 SW samples in comparison with predictions ofF2GEN, 
HERWIG 5.9, "HERWIG 5.9 + power law pt (dyn)" and PHOJET 1.05 (from g^). 



In fig. ^ the usefulness of two choices of the second, beside xsj, unfolding variable 
is studied for the LEPl SW sample. Both studied variables Efor/Etot and / Etot 
are sensitive to the angular distribution of the final hadrons. Similar results for LEPl 
FD and LEP2 SW were obtained as well (not shown). 

Test of the unfolding methods was performed, see figs. ^ and where F2 ob- 
tained using different MC models and unfolding variables is shown for the LEPl SW 
and LEP2 SW samples, respectively. Similar test was performed for LEPl FD data 
(not shown). These results allowed to choose the most reliable unfolding variables, 
which are Xcor and Efor/Etot- 
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Figure 65: The distribution of the observed energy in the forward direction divided by 
the total observed energy, Ejor/Etot (s;C,e), and the transverse energy out of the tag plane 
divided by the total observed energy, E^V-'^/Etot (b,d,f) for three bins ofxms- The comparison 
of the data with the predictions of F2GEN, HERWIG 5.9, "HERWIG 5.9 + power law pt 
(dyn)" and PHOJET 1.05 is shown (from ^). 
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Figure 66: The unfolding of xsj studied for the LEPl SW sample. Results based on one- 
variable unfolding with a) x^s, b) Xcor and on two-dimensional unfolding with c) Xrec dind 
E^^/Etot, d) Xcor and E^^/Etot, e) Xrec and Efor/Etot, f) Xcor and Efor/Etot- The unfolding 
was done using F2GEN, HERWIG 5.9, "HERWIG 5.9 + power law pt" (dyn) and PHOJET 
1.05 (from J4|/j. 
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Figure 67: As in fig. 66 for tfie LEP2 SW sample (from p|/j 



Comment: "HERWIG 5.9 + power law pt (dyn)" and PHOJET 1.05 describe the 
OPAL data better than the F2GEN and HERWIG models. 

To reduce the MG '^modelling error, two-dimensional unfolding has been introduced, 
using E^*' / Etot as a second unfolding variable." 

"Monte Garlo modelling of the final state is still a significant source of systematic 
error, but it no longer dominates all other sources. " 



•TOPAZ 94 g (TRISTAN) 

The jet production (one and two jets) has been studied for the first time in deep 
inelastic cj scattering (3.0 GeV^ < <30 GeV^), see also sec. |2.2.2| . Events with 



the transverse momentum of jets between 2 and 8 GeV and for the \r]\ < 0.7 were 
studied using the jet cone algorithm with R=l (see next section for the definition of 
the jets). The Monte Carlo simulation was based on the VMD formula for the hadronic 
part. For the perturbative (pointlike) part a QCD model was used with FKP (uds) 
and QPM (c) contributions, with p^ =0.5 GeV. 

The pointlike and hadron-like configurations resulting in the different final-state 
jet topologies (two jets and single jet, respectively) were studied. The pointlike part 
was simulated with the same angular distribution as for 77 — > qq. The hadron-like 
part involves the struck quark along the 7* direction and the remnant one along the 
7 direction. 

Comment: The production of events with two high-px jets is consistent with the point- 
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like perturbative contribution. 

"...an excess over the pointlike component is observed in the one-jet sample, which 
is direct experimental evidence for the existence of the hadron-like component in deep 
inelastic 67 scattering" . 
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2.4 Large pr processes in 77 and "fp collisions. Resolved pho- 
tons 



The structure function F2 for the real photon considered in sec. |2.2| is sensitive mainly 
to the combination of quark densities, and, moreover, due to its strong dependence 
on the charge, to the up-type quark distributions. Therefore it is of great importance 
to have an additional and quite different from the DlSg-y-experiments source of in- 
formation. Hard photon(s)-induced processes provide such a source, with the cross 
section given by a combination of different parton densities in the photon(s) convo- 
luted with the cross sections for basic partonic subprocesses (see below). The role of 
these resolved photon processes, e.g. in determining the gluon density in the photon, 
is unique indeed. They are also crucial in probing the heavy quark content of the 
photon, in particular the charm quark density; this topic we discuss separately in sec. 

As we have already mentioned, the large px jet production in 77 collision leads to 
the complementary information on the photon structure to that coming from the F2 . 
Moreover, the study of the hadronic final state, with the production of, predominantly, 
small pt particles (jets), is a standard part of the structure function analyses performed 
in DISe7 measurements at LEP and TRISTAN (see discussion in sees. p.2| and p.3|) . 

The production cross section of the large px jets, which are the hadronic repre- 
sentations of the hard partons produced in the basic subprocess, does depend on the 
applied jet definition. Although this introduces an additional uncertainty in the de- 
scription of the data, the relatively high rate makes large px jet production, apart from 
F2 , the basic source of information on the partonic structure of 7. The main goal in 
such measurements is to extract the individual parton density in the real photon]^. 

In this section we focus on the dedicated large pr single and double jetf^ production 
in the resolved 77 and 7p collisions, where the real photon(s) may interact through 
its partons[5 The bulk of data is coming from the HERA collider, where the pho- 
toproduction of jets with higher transverse momentum than at LEP and TRISTAN 
colliders was analysed. The data on the three or more jets, the photon remnant or the 
jet shape are also available. 

Beside the jet production also the inclusive one-particle production is sensitive to 
the partonic content of the photon, but it depends in addition on the fragmentation 
functions, or, as in the prompt photon production, on an additional assumption (the 
isolation criteria). Nevertheless they play a crucial role in establishing the structure of 
the photon, as the newest charged particle production measurements at HERA, where 
the gluon content has been derived. In addition we collect also the results for the pho- 
toproduction of the prompt photon, which may consist in future an additional source 
of information on the gluon density in the photon. The prompt photon production 
was studied also in 77 collision at TRISTAN (see TOPAZ 98,conf). 

In the previous section we have discussed the single tagged events with virtuality of 



^^At present most of the data deal with the almost real photon; for the jet production in processes 
with virtual photon, see sec. ||. 

^'^We count only jets from hard process, not remnant jets. 

^^The first evidence for the resolved photon contributions and a need for the gluonic content of 
photon in 77 collision can be found in AMY 92. See HI 92 and ZEUS 92 for a first evidence for 
hard processes in the collision. 
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photons radiated by the tagged electrons larger or much larger than 1 GeV^ (with the 
exception of one experiment, DELPHI 95). In this section the initial photons may 
be considered real, as the median virtuality of the photon(s) is smaller than Ag^^, 
while the transverse momenta of the observed jets are chosen to be much larger than 

^QCD- 



2.4.1 Theoretical description 

Below we discuss shortly how the single and double jet production in the photon- 
induced processes is described in QCD. The specific case of the the charged hadron 
photopro duct ion at HERA will be described shortly in HI 99a, and the prompt 
photon photoproduction will be discussed in the separate subsection |2.4.5| . 

The production of large pr jets (or individual particles) may proceed directly (the 
direct or the (Q)PM contribution), or as a result of partonic interaction between 
components of the photon(s) leading to resolved photon contribution. In some analyses 
the soft VMD-type contributions are included separately. 

The resolved photon processes with large pt jets can be characterized by the scale 
of hardness Q^, at which the photon is resolved. It is usually provided by the pt of the 
final jets, i.e. ~ p^. (We introduce here the notation in order to distinguish 
it from the DIS scale Q^, which is equal to the virtuality of the photon probe.) To 
apply the perturbative QCD to such resolved almost real photons, i.e. with ~ 0, 
the scale should be much larger than A.q^j^. Events where virtualities fulfill the 

relation Aq^^ ^ ^ are discussed in the next section, where the concept and 
data on the structure of virtual photon will be introduced. 




Figure 68: The inclusive large px jet production in the resolved 'jp collision (left) and in 
the (double) resolved 77 collision (right). 

Assuming the factorization between the hard subprocess cross section a and the 
parton densities, the generic LO cross section for the two-jet production in jp collision 
(fig. left) is given by 

^^,P^,et, jet, = Y^j j dx^dxpf,/^{x„ Q')fjUxp, Q^)a'^-^''' (20) 

id 

where /j/^(z^, Q^) describes the LL probability of finding a parton of a type "i" in the 
photon. The corresponding function fj/p{xp,Q'^) is related to the proton. The x^(xp) 

^^Thc collection of useful formulae for the cross section a for 2^2 body partonic processes can 
be found in |]60| . 



78 



variable is by definition the part of tlie four-momentum of tlie plioton (proton) carried 
by its parton. The direct interaction of the initial photon corresponds to = 1 and 
f^l^ = 6{1 — x^). The scale of hardness plays here a role of the factorization scale 
between the hard and soft processes in the 7p collision. Note that the precise definition 
of the factorization scale can only be given in the QCD calculation going beyond the 
LO approach. 

For the two-jet production in 77 collision, 77 — > jeti jet^ X (fig. right) a 
similar LO formula holds: 

If in the final state a specific hadron is measured, the fragmentation function will 



appear in addition in formulae (^) and (|2lD , together with the corresponding z vari- 
able, describing the part of four-momenta of final parton taken by the final hadron. 

In the NLO approach, needed to describe the present data, the NLL parton densi- 
ties should be used and the a should contain 0:5 corrections. Note that the separation 
between the direct and resolved photon contribution holds only in the LO approach. 
The theoretical predictions based on the LO - [|17| or NLO (HO) [|17| - [M 



QCD calculations are available for the inclusive one and two-jet photoproduction in 
e^e~ and ep experiments. 

Here we would like to point out that the way of counting the order of the QCD 
calculation in the photon-induced processes is still a subject of discussion (see e.g. 
p7| , ||192|| c and earlier papers ||189|| ef). The origin of the problem is the presence of 



the definite PM prediction, in form of alog((5^) contribution, leading to the inhomo- 
geneous term in the evolution equations for the quark densities in the photon. One 
approach bases on the treatment of the quark density in the photon as being of order 
q'^ ~ a /as, while in the second one g''' ~ a. This difference leads to different sets 
of diagrams which formally should be included in the NLO calculations, see e.g. the 
prompt photon production at HERA, where two types of NLO predictions are com- 
pared to data (ZEUS 2000a). We will not discuss here this issue any further, and in 
the presentation of the results we will use the terminology describing the accuracy of 
the analysis, LO or NLO, as given in the original papers. 

The multiple interaction between partons can also be taken into account in the 
analysis of the jet production in 77 and 7p processes, in analogy to pp processes. 
Since the direct events should be free from such multiple interactions, the 77 and jp 
collisions may offer a unique laboratory to study this problem. The multiple interaction 
seems to be needed for a good description of the data, see below. 



2.4.2 Measurements of the resolved photon processes 

At present, the 77 collisions arising in e~^e~ colliders in no-tag0 or antitagged condi- 
tions correspond to the 77 events where both real photons can be described by the 
Weizsacker - Williams energy spectra. In the OPAL experiment the maximal squared 
target photon masses were equal Pi2= 0-8 GeV^, while in the ALEPH experiment - 
around 7 GeV^, and in TOPAZ - 2.6 GeV^. Still, the median virtuality in 77 collision 

^°The name "untagged" is also used. 
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is small enough {e.g. 10~^ GeV^ in TOPAZ). In the jp scattering at ep collider HERA 
a similar Weizsacker - Williams spectrum describes the flux of photons coming from 
the electron (in year 1996/7 the positron, in 1998 the electron again^)- here < 4 
GeV^ (but with the median 0.001 GeV^), recently also < 1 GeV^ or below 0.01-0.02 
GeV^ if the dedicated detector is used. 

Note that, contrary to the DlSe^ experiments, in the resolved photon processes 
there is no one basic observable, analogous to F^. Usually the pr or Et and pseudo- 
rapidity f] distributions for final particle or jet are compared with the QCD predic- 
tions using a specific parton parametrization. The agreement with the data justifies 
a posteriori the correctness of the applied parton parametrization and of the QCD 
calculation. Only in few cases the parton distribution is extracted from data (and 
only for 7p events, see below). 

For comparison with QCD not only the above distributions are important but also 
the study of the structure of the jets (transverse energy fiow around the jet axis, i.e. 
jet profile ) and the structure of the underlying event (jet pedestal) 0, corresponding 
to the hadronic activity outside the jets where effects due to the beam remnants , 
from spectator partons, may be seen. Also an additional interaction between these 
remnants, leads to the higher energy level of the underlying event. 

For a more detailed study of properties of the resolved photon processes the sepa- 
ration between the direct photon events, where the photon participates directly in the 
hard subprocess, {x^ ~ 1), and the resolved photon events, is needed. 

In the two-jet events, Xp and x^ (for HERA measurements) and x^'^ (for LEP 
and TRISTAN data) distributions, sensitive to parton densities, were studied, see e.g. 
ZEUS 95b, 98a and OPAL 97e. The angular distribution {e.g. in the dijet CM 
system) da /d cos 6*, on the other hand, is expected not to be sensitive to these ingre- 
dients. So this measurement may help to verify the expectation of different angular 
distributions due to different subprocesses which correspond to the direct and resolved 
photon contributions, see ZEUS 96 for the first results, also OPAL 99a. The first 
extraction of the gluon density (HI 95a) and of the effective parton distribution (HI 
98) was reported by the HI Collaboration. The data on prompt photon production 
have appeared recently both from the TRISTAN and HERA colliders. 

The main Monte Carlo generators used by the experimental groups for these kinds 
of processes are the general purpose generators PYTHIA, HERWIG and PHOJET 
(with or without the multiple parton interactions). The DELPHI collaboration uses 
the TWOGAM generator to simulate the QPM, VMD and resolved photon contribu- 
tions. All of these MC programs use the cutoffs, separating the soft and hard 
processes in the simulation, being fixed by the requirement to reproduce the total cross 
section. (In case of the multiple interaction p™* parameter is used instead of p^*" ). 
For both e^e~ and ep colliders the additional power law pt spectrum has been included 
recently in the HERWIG and PYTHIA programs. They were used in e.g. ALEPH 
97a,conf, 99a, 99b,conf; OPAL 2000; TOPAZ 98,conf; ZEUS 95c (where the 
spectrum was first introduced). 

The single and double jet samples were studied, reaching the highest px of jets in 

^^Note that starting from 1998 new energies are accessible: 27.5 GeV for e" and 920 GeV for p, 
see 2.3.4. 

^^z.e. the transverse energy determined outside of the jet cone 
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77 processes up to 17 GeV. Jets at the HERA collider are observed up to 75 GeV 
(ZEUS 98c,conf and 99b,conf). 

The jet distributions are in general in agreement with the expectations for both 
single and double jet production, both at e~^e~ and e^p colliders. In the pseudorapidity 
dependence the discrepancy is observed for the jet rates. The transverse energy flow 
around the jet axis, especially for small and small is not properly described 
by existing Monte Carlo generators. Also for the two-jet events it has been found 
that distribution at the small lacks proper description. The agreement with the 
data is improved when the additional power law pt spectrum is taken into account in 
the analysis, as in F2 and hadronic final state analysis. A need for extra pt in the 
simulation of the whole event can also be taken as a hint that the multiple scattering 
may be important in the photon-induced processes. The multiple interaction included 
in the experimental analysis is modelled usually as in the pp experiments (see e.g. 
HI 96a). Global study of the photon-induced jet production with the aim of tuning 
the general purpose Monte Carlo models was performed recently in []T2[. "The main 
ingredients investigated in the tuning are the description of the "underlying event" and 
the choice of photon structure" . The conclusion of this analysis is that "the agreement 
without the multiple interaction is very poor" . 

The strong dependence on the proper choice of the jet definition and the parameters 
like the cone size R (see below), has also been observed in the jet production at HERA. 
There is a possible relation of this effect to the problem of describing the underlying 
events observed both in e~^e~ and in ep collisions, that was mentioned above. 

The following two jet definitions are used in the analyses reported below. The 
first is the jet cone algorithm with the fixed value of the cone variable R, defined as 

R = \J (50)^ + (5//)^, with 50 and 5rj describing the differences between the cone (jet) 
axis and the particle direction in the azimuthal angle and pseudorapidity, respectively. 
It is used in PUCELL and EUCELL algorithms. The second one corresponds to the 
kr - cluster algorithm on which the KTCLUS approach is built up. Beside R, also 
Rsep is included in some analyses using the cone algorithm. It corresponds to the 
additional separation between partons, and the Rsep = 2 i? means no restriction. (The 
discussion of jet definitions in the NLO calculation can be found e.g. in [rS7].) 



The difference of the two jets' transverse energy AEt = Eti-Et2 and analogous 
difference for the pseudorapidities, Ar] = jr^i — are introduced in two-jet events. 
(For the analysis of the energy of the underlying event the variable 67] = rjceii — rj is 
used.) In the analysis of two-jet events there appear also the average pseudorapidity 
fj = {rji + 772)72 (denoted also by 77jet(s)) and the average transverse energy of jets, Et 

{Ejet(s))- 

It is worth noticing the difference between the variables describing the produced 
parton and the corresponding quantities for the jet, representing the considered parton 
on the hadronic level. The unfolding procedure is needed in order to obtain the 
"true" partonic variables, e.g. the transverse momentum of the parton (pt) from the 
transverse momentum or energy of the jet {pt or -E't)0- 



23 



We keep in subscripts small letters for quantities related to the partonic level 
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2.4.3 Jet production in 77 collisions 



Here we present the results for the jet production in 77 collisions, where one or two 
resolved photons may interact. 

As in the previous section we will denote the virtualities of two involved photons 
by and Pg^. The antitagged events at LEP correspond to < 4 GeV^, with 
median approximately 10~^ GeV^. For the antitagging conditions in OPAL 97e the 
maximum value P^ = 0.8 GeV^ was used. At TRISTAN the P^ value in antitagging 
case lies between 10~^ and 2.6 GeV^ with the mean value ~10~^ GeV^. 

In e^e~ machines the parton momentum fractions, variables x^'^ = x^, are deter- 
mined from the two final jets with the highest Et-, according to the formulae: 

< = y""^'',^l'\ - (22) 

^hadronsK^ ^ Pz) 

Variables referring to the dijet centre of mass system are (in OPAL 99b) denoted by 
star, e.g. rj*. 

Below we present the LEP and TRISTAN results. 



DATA 

•ALEPH 93 [H] (LEP 1) 

The presence of hard scattering processes is demonstrated in the sample from 1990- 
1991 data. The antitagging condition corresponds here to < P^^2 > — 0-23 GeV^. The 
VMD, QPM and QCD contributions, the latter with p^*" cutoff, were studied using 
Vermasseren program for QPM, and TOPAZ MC generator for QCD part. The p^*" 
was fitted to be 2.5±0.1 ± 0.3. 



Comment: "The data are well described by a combination of three models, namely 
VMD, QPM and QCD". 

The fitted parameter is "higher than the values previously obtained at TRISTAN 
but consistent with recent HERA photoproduction results" . 



•DELPHI 94 ^] (LEP 1) 

Data were collected in the period 1990-92 with < Pf >~ 0.12 GeV^. Events are 
groupped into a minimum bias sample (sample I) and a high pt jet events sample 
(sample II). Jets, defined according to the Lund cluster algorithm, with px greater 
than 1.75 GeV, and a polar angle between 40° and 140°, were observed. The three 
component model (VMD for a fully nonperturbative contribution, QPM, and QCD 
with the p^™) was used in the analysis. 

Results for the jet pt distributions are shown in fig. To model the 
(resolved photon) contribution the FKP approach 



QCD 



see sees. 



2.2.2, 2.3) was used 



with DG and DO parton parametrizations. For a comparison also the LACl and GS 
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densities were used. 
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Figure 69: Distribution of jet transverse momentum compared with the Monte Carlo pre- 
dictions (TWOGAM generator). Dark area - QPM only; dots - QPM + VMD; Ml line - 
QPM + VMD + DG (p??*" = 1.45 GeV/c); dashed line - QPM + VMD + DO (p^*" = 1.22 
GeV/c) (from 1^). 




Comment: " A three-component model which includes QCD hard scattering subpro- 
cesses is needed to describe the data. " 

Remnant jets, generated along the direction of the initial quasi-real photon, are stud- 
ied within the FKP approach. 

"Each parametrization is associated to a special value of p™*", constrained by the 
description of the visible total cross section". 

•DELPHI 95 (LEP 1) 

The measurement of hard processes in 67 colhsion was performed based on 1991-1992 
data, as described in the previous section. The average virtuahty of the photon from 
the tagged electron, registered by VSAT, was equal to < Q'^{= Pf) >~ 0.06 GeV^ 
(single-tag events). The values of Et of observed jets were larger than 1.5 GeV. 

Comment: "This study confirms recent results from no-tag experiments, indicating 
that the photon has a significant partonic content. " 



•DELPHI 97b,conf [H (LEP 1, LEP 2) 

Based on the LEPl and LEP2 data (1995 and 1996 runs), with CM energy 172 GeV, 
various distributions for the hadronic final state in anti-tag events were measured. The 
TWOGAM generator was used in the analysis with soft VMD, QPM and QCD-RPC 
(with p^*"- cutoff) contributions. Different parton parametrizations for the photon 
were used assuming the "leading order QCD factorization", which lead to different 
values of p??'": p7"(G'52) = 1.88±0.02 GeV and p?^*"(GW) = 1.58±0.018 GeV. 
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Results are presented in fig. |7D|. 
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Figure 70: Distributions of: a), d) invariant mass, b), e) pT, c) number of reconstructed 
jets, f) pseudorapidity. Solid and dotted lines - predictions based on VMD + QPM + RFC 
with the GS2 and GRV parton densities, respectively. The dashed line - RFC contribution 
with the GS2 parton densities; the hatched histogram - the VMD part, the double hatched 
- the QPM part. Upper (lower) figures correspond to ^/s=91 (172) GeV (from 



Comment: "Monte Carlo prediction (based on TWOGAM generator) gives perfect agree- 
ment at CM energy = 91 CeV and slightly exceeds data at 112 GeV " . 



•OPAL 97e [|5p (LEP 1.5) 

The inclusive one and two-jet cross sections have been measured at y^s^ = 130 and 
136 GeV (based on the 1995 run). The antitagging condition was applied to the ini- 
tial photons corresponding to the maximum photon virtuality P^„^ ^ 0.8 GeV^. The 
transverse energy of jets E't is taken to be larger than 3 GeV and the pseudorapidity 
lies within \'!]\ < 1. Generators PYTHIA 5.721 (with the SaSlD parton parametriza- 
tion) and PHOJET 1.05 (with the GRV LO parton parametrization) were used. The 
cone jet-finding algorithm with R=l was used (for the first time in photon - photon 
collisions at LEP). 

For the two-photon events, the distribution and the transverse energy flow 
around the jet, studied separately for the direct (mm{x1^, x~} > 0.8) and the resolved 
{mm{x^ ,x~} < 0.8) photon contributions, were obtained (see fig. |7l] and fig. |7^, 
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respectively). 




Figure 71: The number of two-jet events as a function of compared to the PHOJET 
(sohd hne) and PYTHIA (dashed hne) simulations. The hatched histogram is the direct 
contribution to PYTHIA events (from j6^/ j. 
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Figure 72: The transverse energy Eow around the jets in two-jet events, b), d) Direct events 
(min{x^ , x~} > 0.8). a), c) Resolved events (min{x^ , x~} < 0.8). The data are compared 
with the PHOJET (solid line) and PYTHIA (dashed line) simulations. The rapidity 5rj 
distribution is shown in a) and b), while the 6(f) angle distribution in c) and d) (from j§^). 



The one-jet and two-jet cross sections were measured up to Et= 16 GeV. 



Results are presented in fig. 73 
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Figure 73: The inclusive one-jet (left) and two-jet (right) cross section as a function of Et 



for the jets with \r]\ < 1 compared to the NLO calculations ^182]. The solid line is the sum 
of direct, single resolved and double resolved contributions that are also shown separately 
(from MW). 



The r] distributions for one and two-jet events were studied as well (see fig. |7^ . 
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Figure 74: The inclusive one-jet (left) and two-jet (right) cross section as a function 
of \r]\ for the jets with Et > 3 GeV compared to the LO QCD calculations of the 
PYTHIA (with the SaSlD, LACl and GRV parton distributions) and PHOJET (with 
the SaSlD parton distributions) generators (from 



Comment: "The data on rj distributions agree well with NLO QCD calculations based 
on GRV parametrization. The GRV LO and SaSlD parametrizations describe the data 
equally well, the LAC 1, however, gives twice the observed value" . 
The large difference between Monte Carlo generators with a fixed parton parametriza- 
tion, and for a particular Monte Carlo generator with different parton parametrizations 
is observed. 



•OPAL 99a [||] (LEP 2) 

The dijet production in two-photon collisions at e'^e~ energies 161 and 172 GeV in 
the "antitagging condition" was measured. The transverse energy of the jets is taken 
to be Et > 3 GeV and the pseudorapidity lies within \r]\ < 2. The Monte Carlo 
generators PYTHIA 5.722 and PHOJET 1.05c were used with the SaSlD and GRV 
LO parton parametrizations, respectively (with p^ = 1.4 GeV, as a default value). 
Both generators use the multiple interaction to model the underlying events (a default 
p^*=1.4 GeV). The cone jet-finding algorithm with R=l was used. 

The contributions of direct and resolved photon subprocesses (x^ < 0.8) were stud- 
ied for various distributions. 

The distributions were measured for different Et bins, from 3 to 20 GeV, 
showing the large contribution of the resolved photon processes (fig. [75]) . 
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Figure 75: The uncorrected x-y distributions in bins of Et, compared to PYTHIA and 
PHOJET predictions (from Mj). 



The transverse energy flow around the jet axis is presented in fig. ^ for the double 
resolved and direct photon processes. Larger hadronic activity for the resolved photon 
sample is observed as expected. 
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Figure 76: The transverse energy Bow around the jet axis as a function of 5r] in two-jet 
events; a) double resolved events, (b) direct events. Comparison with PYTHIA (solid line) 
and PHOJET (dashed line) (from ^). 



The angular dependence plotted separately for the three basic resolved photon 
subprocesses and for the direct events is presented in comparison with matrix element 
predictions [^| in fig. 77. A comparison with the NLO calculation with the GRV HO 
parton parametrization was also made (not shown). 
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Figure 77: The angular distribution of events with Eti = Et2 and colhnear in (p, separated 
into the direct and the double resolved contributions, together with the LO QCD expectations 
based on matrix elements pD/ (from |6^j. 



The cross sections for the direct, single resolved and double resolved two-jet events 
versus Et, compared to the NLO calculation with the GRV HO parton parametriza- 
tion, are shown in fig. [78| . 
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Figure 78: The inclusive two-jet cross section as a function of Ex for jets with \ri\ < 2, 
compared to the NLO calculations , 18^ . Solid line is the sum of the direct, single resolved 
and double resolved cross sections that are also shown separately (from /|6^/). 



The infiuence of the "underlying event" has been studied. The transverse energy 
flow outside the jets was measured as a function of x^, in order to study the possible 
contribution from multiple interactions. The Monte Carlo simulations (PYTHIA, 
PHOJET) were performed with various parametrizations and values of the - cutoff 
parameter for multiple interactions (fig. [791). 



The two-jet pseudorapidity distribution is presented in fig. |80| together with the pre- 
dictions from the PYTHIA and PHOJET generators, with various parton parametriza- 
tions. Here the asymmetric cuts were applied, to avoid infrared instability in the NLO 
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calculations. The results are also shown for two samples, corresponding to the direct 
and double resolved photon processes. 
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Figure 79: The transverse energy Bow outside the jets for \r]*\ <1 as a function of for 
various parton parametrizations, MC models and values of (from pW). 
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Figure 80: The inclusive two-jet cross sections as a function of |r/| for jets with E^^^ > 5 
GeV and E^^'^ > 3 GeV: a) for all events, b) for direct events, c) for double resolved events 
(from j6^/j. 



To reduce the influence from multiple interactions and hadronization effects the 
two-jet cross section was measured also for E^*'^ > 5 GeV (not shown). 

Comment: "The Et dependent two-jet cross section is in good agreement with NLO 
QCD calculation. " 

"Using PYTHIA and PHOJET the GRV LO parametrization is able to describe the 
two-jet cross section; prediction based on the SaSlD is too low, and that for LACl - 
significantly too high. " 



89 



"Within the errors of the measurement we are unable to differentiate between models 
with and without multiple interactions. " 



•AMY 1.0 92 (TRISTAN) 

The measurement of the high px hadron production in the quasi-real 77 colhsion, for 
the e"'"e~ colhsion energy between 55 and 61.4 GeV, was performed. In the observed 
three and four-jet events there are one or two spectator jets coming from the resolved 
photon(s). The QPM+VMD contribution, and the MJET one, which corresponds to 
the RPC contributions, were included in the analysis. The DG parton parametriza- 
tion was used with the p^^^ equal to 1.6 GeV. Various distributions were studied, the 



results for the px distribution are shown in fig. 

The sensitivity of the results to the gluonic content of the photon was studied. 




P^r^' (GeV/c) 

Figure 81: The pt distribution compared with the predictions of QPM (dotted histogram), 
QPM + VMD (dashed histogram) and QPM + VMD + MJET (soUd histogram) with p??™ 
= 1.6 GeV (from §^). 



Comment: The "first clear" evidence is found for the hard scattering of hadronic con- 
stituents of photons in the photon - photon collisions at TRISTAN, 
"...without the gluonic component it is impossible to reproduce the data". 

•AMY 1.5 94 (TRISTAN) 

AMY 1.5 (the upgraded AMY detector at the energy 60 GeV) measured the inclu- 
sive single and double jet cross section ^ in quasi-real 77 collisions. Jets with px 
above 2.5 GeV up to 8 GeV and \t]\ < 1 were studied. The QPM + GVMD (for 
the diffractive production of the hadrons) and the improved MJET generator for the 
resolved photon processes (using DG, LAC 1,2,3 and GRV parton parametrizations, 
and various p^*" parameters) were applied in the analysis. The cone algorithm with 
R=l was used. 

In fig. ^ the px distributions for one and two jets, for |?7| < 1.0, are presented. 
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Figure 82: The inclusive jet cross section as a function of px integrated over \r]\ < 1.0 for 
one-jet (left) and two-jet events (right). The curves represent the sum of the QPM (direct) 
and MJET (resolved) cross sections using the parton paranietrizations: LACl (full line), 
GRV LO (double dot - dashed), DG (dotted line), LACS (dashed line) and LACl without 
the gluon component (dot - dashed line). The short - dashed curve corresponds to the QPM 
cross section (from 



The r] dependence for the jet with px larger than 2.5 is presented in fig. q3 
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Figure 83: The inclusive jet cross section as a function of rj integrated over \pt\ > 2.5 GeV. 



The lines convention is the same as in fig. p2 (from [68]). 



Comment: "The data are in good agreement with LO QCD calculations based on either 
LAC 1, DG or GRV parametrizations of the parton densities in the photon. The 
calculation based on LAG 3 disagrees with the data and we do not see the deviation 
that is observed in the rj distribution by the HI experiment" (here denoted as HI 93j. 

•TOPAZ 93 [H] (TRISTAN) 

The inclusive jet cross section was measured at the collision energy y^s^=58 GeV 
and virtuality of photon below 2.6 GeV^. Jets with px between 2.5 GeV and 7.5 
GeV and \ri\ < 0.7 were analysed. The direct and resolved photon contributions were 
included, as well as the VMD one. The resolved contributions were calculated with 
various parton parametrizations (DO, DG and LACl, 2, 3), using the p™*" parameter. 
The BASES/SPRING generator was used. To describe jets, the cone algorithm with 
R = 1 was apphed. 
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The energy flow within the jet (per jet) as a function of 6r] and 6(f) was also studied 
and found in agreement with the Monte Carlo simulation (not shown). 

In figs. |8^ and |85| the pt distribution of the jet and the energy flow per event as a 



function of the polar angle of the jet (9) are presented, respectively. Fig. ^ shows the 
inclusive one-jet and two-jet cross sections. 



t3 




Figure 84: The jet pt distribution for a no-tag sample. Predictions are shown for the sum 
of direct + resolved photon (with the DG parton parametrization) + VMD processes for 
different p^*" values. Contributions of the direct (dashed line) and the VMD (dot - dashed 
line) processes are shown. Parton parametrization according to DG (from 
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Figure 85: The energy flow versus the polar angle 9 for the jet sample. The data are 
compared with the Monte Carlo predictions for the sum of the direct and the resolved 
processes (solid histogram) and the direct process only (dashed histogram) (from ^/ }. 
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Figure 86: The inclusive a) jet and b) two-jet cross sections as a function of px for the 
pseudorapidity \ rj \< 0.7. TOPAZ Collaboration results are compared with the predictions 
of different parton parametrizations: LACS, DO + VMD, DG, DG (without gluon), LACl 
and LAC2. Also a direct contribution is shown (from J6^/ ). 



Comment: The polar angle distribution gives .."direct evidence of the presence of the 
resolved processes. The data exclude the parametrizations of LAC 3 and DO + VMD, 
which predict a very large gluon content even at large x^." 

• TOPAZ 95,conf ^ (TRISTAN) 

The updated results of the jet study in (quasi) real 77 collisions are given, for a 
larger data sample than in [^. Anti-tag condition corresponds to the maximum 
photon virtuality of P^ax = 2.6 GeV^. Production of jets with pt between 1 (2) 
and 8 GeV, and |?7| < 0.7 was measured. The direct (QPM) + VMD + resolved 
(RFC) contributions were included in the analysis of the data. The parton density 
parametrization used is that of DG, and LACl, with optimalized cutoff parameters 

FT ■ 

The Pt distribution for the inclusive jet production has been studied down to 1 
GeV, and the data show clearly the evidence for the resolved processes (see fig. |87|) . 
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Figure 87: The inclusive jet distribution as a function of pT for \ri\ < 0.7. The solid line is 
a prediction of the direct + resolved (with DG parton parametrization) + VMD processes 
^pmin _ -^ j GeV). The contributions of the direct and of direct + VMD processes are shown 
as dashed and dot-dashed histograms, respectively (from ^/). 



The evidence for the remnant-jet activity in small-angle region was observed. 

Pt jets in the 



Fig. |8§ shows the particle energy flow for the sample with two high 
central region, where the energy deposit from direct (or VMD) processes is rare. 
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Figure 88: The energy How distribution for the central two-jet sample with |r/| < 0.7 and 
Pt^"^ GeV. Comparison with prediction of Monte Garlo models as in fig. (from [TC]). 



Comment: The existence of the resolved processes in 77 collisions is confirmed by the 
jet rate and the direct detection of the energy clusters coming from the remnant jet. 
The rate of the remnant jet is also consistent with the Monte Carlo expectations. 
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• TOPAZ 98,conf ^ (TRISTAN) 

A study of dijet production in (quasi) real 77 collisions (anti-tag conditions: from 
10"^ to 2.6 GeV^ with the mean value ~ 10"^ GeV^) at ^j;; = 58 GeV based on the 
full data taken in the years 1990-95 was performed. (Also a prompt photon production 
cross section was measured.) Limits: > 3 GeV and > 2 GeV for two jets with 
the highest transverse energy were applied. The Monte Carlo programs PYTHIA (also 
with additional power law pt distribution) and the PHOJET with GRV LO and SaSlD 
parton distributions were used. The jets were reconstructed using a cone jet-finding 
algorithm with fixed value of cone radius R = 1.0. 

The study of the hadronic final state (the transverse energy distribution, the jet 
energy profile and the azimuthal angle dependence of the number of dijet events) was 
carried out by separating the direct and resolved sample by a reconstruction method 
developed at HERA jp experiments (see e.g. sec. |2.2| and ^.4.4|) . The number of dijets 
as a function of is shown in fig. ^ The predictions of the PYTHIA generator with 
the GRV LO structure function agree with the data, but that with SaSlD are too 
low by a factor 2 for x™*" (defined as mm{x^ ,x~)) ~ 0.5. This observation appears 
for both PYTHIA and PHOJET programs, and also for different values of parameter 
p™", varying from 1.6 to 2 GeV. 



a 
2 



© 

s 

s 

z 



400 



350 



300 



250 



200 



150 



100 



50 



0,1 0,2 0.5 0.4 0.5 0,6 0,7 0.8 0.9 



X 



2 



.Si 
s 

© 

s 

S 

z 



400 



350 



300 - 



250 



200 



150 - 



100 



50 



Data 

PYTfflA GRV (P^""=1.6GeV) 
PYTHIA GRV (P^''"=2.0GeV) 
PYTfflA SaSlD (P^™=1.6GeV) 
PHOJET SaSlD (P7"=2.0GeV) 




0,2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Figure 89: The x!^*" distribution of the dijet sample. Left: Comparison with the PYTHIA 
predictions with the GRV LO photon structure function; the shaded histogram - contribution 
from the direct process, the dashed one - from the resolved process, the dotted one - sum 
of both. Right: The same distribution compared with the various Monte Carlo predictions 
(with p% = 1.6 GeV or 2 GeV) (from /[^j. 



The measured Et distributions for the resolved - enriched sample (x™*" < 0.8) and 
for the direct - enriched sample (x™" > 0.8) are shown in fig. 00. 



95 




3 X5 4 4.5 5 5.5 6 6.5 7 7.5 S 

[GeV] 



3 3.5 4 4.5 5 5,5 6 6,5 7 7,5 8 

[GeV] 



Figure 90: The Et distributions for (a) the direct - enriched sample and (b) the resolved - 
enriched sample. The histograms are the PYTHIA predictions with the GRV LO structure 
function (from /[Zl]/). 



The transverse energy flow as a function of the pseudorapidity 61] was measured 
with respect to the jet direction for different dijet samples, see fig. 0. The direct - 
enriched sample is nearly symmetric and there is almost no activity outside the jet, 
while the resolved - enriched sample shows considerable activity outside the jet cone 
for \6ri\ > 1. The energy flow outside the jets is well modelled by the PYTHIA gener- 
ator. 
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Figure 91: Transverse energy flow relative to the direction of each jet in dijet events as a 
function of 6ri for (a ) direct - enriched and (b ) resolved - enriched sample. The histogram is 
the result of the PYTHIA simulation (from 



The transverse energy flow versus the azimuthal angle around the jet direction is 



shown in fig. 92 for the two above samples together with the Monte Carlo predictions 
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with and without the multiple interactions (MI). The activity observed for the resolved 
- enriched sample is much larger than for the direct - enriched; at \5(f)\ ~ 1.5 the 
difference is by a factor 8. 
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Figure 92: The transverse energy How relative to each jet direction as a function of the 
azimuthal angle 5(f) for (a) direct - enriched and (b) resolved - enriched sample. The PYTHIA 
(with and without multiple interaction, solid and dashed line, respectively) and PHOJET 



(dotted line) predictions are shown (from [71]). 



The dependence of the azimuthal angle difference between two jets with the highest 
Ej- was studied. Here the effect of the primordial pt distribution of partons - gaussian 
and the power-law type - was analysed using the PYTHIA program (with the GRV 
parton parametrization). These data (not shown) are reproduced very well when the 
initial state radiation (IR) is included and the power-law pt distribution is assumed. 

Comment: Too low prediction for distribution from the PYTHIA and PHOJET 
programs if based on the SaSlD parton parametrizations. 

"The PYTHIA predictions on the direct and resolved processes reproduce the trans- 
verse energy flow reasonably well. While the PHOJET prediction for the resolved sam- 
ple shows too small activity in the middle region between two jets, the direct sample is 
reproduced well. " 

The transverse energy flow around the jet shows that... "In 77 collisions at TRIS- 
TAN, the PYTHIA predictions with and without multiple interactions are almost iden- 
tical, indicating that there are no effects of the multiple interaction at the TRISTAN 
energy region" , although they have been found at HERA in 7p reaction. 

The data on the azimuthal angle between jets are very well described when the power 
law Pt distribution in PYTHIA is assumed. 
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2.4.4 Jet production in 7j9 collisions 



Here we present the data for the jet production in jp scattering taken at HERA collider 
at y/s^ ~ 300 GeV 0, where only one photon may be resolvedQ. The photoproduction 
events correspond here to the limit of virtuality of the initial photon: ~ 4 GeV^ 
(with median 0.001 GeV^), at present also ~ 1 GeV^ or (with a special condition) ~ 
0.01-0.02 GeV^ (see below). We will concentrate below on the single and double jet 
production and in addition we will include in the presentation some results for the 
charged particle production in the context of extraction of the gluonic content of the 
photon. 

The partonic variables, related both to the initial and final states in the hard 
partonic subprocesses, are reconstructed from the corresponding quantities for the 
final state hadronic jets. The relation between variables corresponding to these two 
levels depends on the order of the perturbative QCD calculation, as was mentioned 
before. The parton momentum fraction x^*^* (called also x™* or x"*^) is in the LO case 
equal to x^. In practice it is reconstructed using two jets with the highest transverse 
energy Et (or E^^^) in the event, using the following relation: 

" ' ^^^^ 

where also jet pseudorapidities, r] = — ln(tan(6'/2)), and the energy of the photon E^ 
(= yEe) enter. (The positive pseudorapidity corresponds to the proton direction.) The 
scaled energy y of the initial photon (eq. |^) is measured from the transverse energy 
E!^ and pseudorapidity of hadrons rj^ according to the formula 



E^Te"^ , (24) 



2^e h 

where the sum is over the produced hadrons. 

The rapidity in the jp CM system is usually denoted by 1]* (typically rj — rj* ~ 2, 
where rj is the value in the laboratory frame). In the angular distribution analysis 
(ZEUS 96, 98d,conf) the polar angle in the dijet CM system, 9*, is introduced. 



DATA 

•HI 92 [|72l (HERA) 

The evidence for the hard photoproduction of a jet and single particle with Et > 10 
GeV is reported. The "soft" contribution and the hard processes were simulated, the 
latter ones using the PYTHIA program with the MT Bl parton parametrization for 
the proton and DG for the photon, and p™" = 2.5 GeV. The jet cone algorithm with 
R < 1 was used in the analysis. 

^''The energy of the electron at the beginning of running of the HERA colhder was equal to 26.7 
GeV; starting from 1994 the positron replaced the electron and its energy increased to 27.5 GeV, 
with the energy of the proton 820 GeV. Starting from 1998 the electron beam with energy 27.5 GeV 
is used again while the proton energy is higher than before, 920 GeV. 

^^The 77 events leading to the large pT jets are rare at HERA. 
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•HI 93 J73l (HERA) 

First measurement of the inclusive jet cross section at the ep colhder HERA (based 
on the 1992 data) is reported for the photon virtuahty smaller than 0.01 GeV^. 
Events with the scaled energy y of the initial photon between 0.25 and 0.7 were col- 
lected. The photoproduction of jet was studied for Et from 7 to 20 GeV and the 
pseudorapidity interval —1 < t] < 1.5. The PYTHIA 5.6 program (with the GRV LO 
parton parametrization for the proton and GRV LO for the photon) with the jet cone 
algorithm assuming R=l was used in the analysis. 

The transverse energy flow around the jet axis was studied and the discrepancy 
was found in a form of too large transverse energy flow on the forward side of the jet 
(see fig. ^3]). Among others, multiple parton interactions were mentioned as a possible 
explanation of this effect. 




A?) i<p 

Figure 93: The transverse energy Bow as a function of 5r] (integrated over \5<f)\ < 1.0) (a, 
c) and as a function of Sep (integrated over \6r]\ < 1.0) (b, d). Figs, a) and b) correspond to 
-1.0 < 7] < 0.5, c) and d) to 0.5 < rj < 1.5 (from f^). 



The inclusive jet cross sections da/dE^ versus Ej- and da/dt] as a function of rj 
integrated over the corresponding range of the rj and Et, respectively, were measured 
and compared with the LO calculation using the following parton parametrizations: 
LAC2, LACS and GRV LO for the photon, and GRV LO for the proton (see fig. |93). 
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Figure 94: The inclusive jet Et spectrum integrated over the rj interval -1.0 < rj < 1.5 (a) 
and inclusive r/ spectrum (h) for jets with Et > 7 GeV. The LO QCD predictions from the 
PYTHIA generator using parametrizations LACS (dashed line), LAC2 (dash-dotted line), 
GRV LO (full line) and GRV LO (without gluons, dotted line) (from j^). 



The shape of the da/dEx is well described in the range of rj between -1 and 1.5. It 
is a problem to describe the r/ distribution of jets. 

Comment: "In the (pseudorapidity) range 0.5 to 1.5 the data show larger average val- 
ues of the transverse energy flow outside the jet cone on the forward side of the jet 
than predicted by the Monte Carlo. " 

LACS gives cross section higher by factor 3 than data for the Et distribution of the 
jets. 

"None of the models describe well the measured rj dependence (for jets)." 

This discrepancy is absent in HI 96a analysis, where multiple interactions are added 

in the Monte Carlo programs. 

•HI 95a []7|] (HERA) 

The photoproduction of two-jet events (the 1993 data) was studied for the photon 
virtuality smaller than 0.01 GeV^. The scaled energy y of the initial photon was 
between 0.25 and 0.7. The jet Et range from 7 to 20 GeV and the pseudorapidity 
interval < ?7 < 2.5 with |A?7| < 1.2 (between the most energetic jets) (see also HI 
96a) were investigated. 

The cone jet-finding algorithm with R=l (and 0.7 for cross checks) was used. The 
PYTHIA 5.6 generator with the GRV LO parton parametrizations for the proton and 
the photon, with pjp^^=2 GeV was used in the analysis. 

For the first time the inclusive (LO) cross sections were derived for the parton level, 
and the gluon density in the photon was measured. 

The transverse energy flow around the jet direction per event versus the rapidity 
distance from the jet direction was studied for 7 < Et < 8 GeV and < rj < 1 and 
found to be asymmetric and different for samples with > 0.4 and < 0.4. The 
data were compared with the PYTHIA (with MI) predictions (not shown). 

The transverse energy flow versus the azimuthal angle around the jet direction and 
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the transverse energy of the underlying events outside the jets (here named Ef^ ^'^ " ) 
was studied for 7 < Et < 8 GeV in a cone size R = 1. Results (see fig. |9^) compared to 
the prediction of the PYTHIA program (with and without MI) lead to the conclusion: 
"the multiple interaction gives an improved description" . 
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Figure 95: (a, c) The transverse energy How versus the azimuthal angle with respect to the 
jet direction in two rapidity bins, (b, d) Distributions of the transverse energy measured 
outside of jets. Histograms show the PYTHIA simulations with (full line) and without 
(dashed line) the multiple interactions (from J^j. 



To achieve the goal which was here the extraction of the gluon distribution, the 
single parton cross section da/dpt integrated over the parton rapidity range, as well as 
the single parton cross section da/drj were studied and compared with the prediction 
based on the LO parametrizations for the photon: GRV, LAC 1 and LAC 3 (see 

fig. m). 
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Figure 96: The single parton cross sections: (top) da/dpt integrated over the pseudorapidity 
range < rj < 2.5, (bottom) da/drj forpt > 7 GeV. The sohd hne - the LO QCD calculation 
with the GRV LO parametrization for partons in the proton and the photon. The dashed 
(dotted) line - the same for the LACl (LACS) parametrization for the photon (GRV LO for 
the proton) (from |^/). 





• HI 









To extract information on the subprocesses, the full two-jet kinematics was used 
and the distributions of Ar], AEt, and of and Xp were studied. In fig. ^we present 
the distribution of x^. 
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Figure 97: The distribution of Xy in the photon. The solid line - the contribution from the 
quark-resolved photon processes; the dashed line - the direct photon contribution, from the 
PYTHIA Monte Carlo (from ^). 



The (LO) gluon distribution in the photon was derived, at the average (factoriza- 
tion and renormalization) scale < >=< pf >=75 GeV^ for 0.04< x^ <1, see table 
23| and fig. ^ for results. 



102 



Table 23: 





< > 


T 1 T 1 / rv 






(stat. + syst.) 


75 


0.059 


1.92 ± 0.87 ± 1.68 




0.14 


1.19 ±0.34 ±0.59 




0.33 


0.26 ±0.24 ±0.33 




0.59 


-0.12 ±0.15 ±0.33 




0.93 


-0.08 ±0.61 ±0.30 




Figure 98: The gluon distribution extracted from the resolved photon processes at 

comparison the GRV LO (full line), LACl (dashed) 



and LACS (dotted) gluon parametrizations are shown (from 



Comment: "The multiple interaction option gives an improved description" of the jet 
profiles and pedestal distributions ( still some "deviations from data at large jet rapidi- 
ties 2 < rj <2.5" remain). 

In the extracting the gluon density cp was taken as determined by the two-photon ex- 
periments at LEP and TRISTAN, in the form given by the GRV LO parametrization. 
"A high gluon density at large parton momenta as suggested by the LAC 3 parametriza- 
tion is clearly excluded. The strong rise of the LACl parametrization below x-y < 0.08 
is not supported. " 



•HI 96a []7|] (HERA) 

The jet production with > 7 GeV (and -1 < r/ < 2.5) was measured in ep colhsions 
(data from 1994) with the scaled photon energy 0.25 < y < 0.7 and below 0.01 
GeV^. The properties of the hadronic final state and the distribution of the transverse 
energy are studied in detail. The PYTHIA 5.7 (with p™" = 2 GeV), HERWIG 5.8 
(p^™ = 2 GeV) and PHOJET 1.0 (p^*" = 3 GeV) generators with the GRV LO 
parton parametrizations for the proton and the photon) were used in the analysis. 
The PYTHIA program were used also with MI, then with p™* = 1.2 GeV. The cone 
algorithm with R=l was used. 

The integration over the 7p CM system pseudorapidity -2.5 <ri* < 1 leads to the 
total transverse event energy distribution shown in fig. p9| a (here 0.3 < y < 0.7). The 
average transverse energy flow versus rj* for the total Et range between 25 and 30 
GeV was also measured (see fig. |99|b). The shape of both distributions may indicate 
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the need of the muhiple interactions. 
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Figure 99: a) The differential transverse energy cross section integrated over the pseudora- 
pidity (-2.5 < f]* < 1). Tfie fiistograms are tlie results of the simulations with interactions 
of the beam remnants (full line - PHOJET, dashed - PYTHIA) and without them (dotted 
line - PYTHIA). b) The corrected transverse energy How versus r]* (rj* > corresponds to 
the proton direction). The pseudorapidity range and histograms as in a) (from 



To get an insight into the details of the considered events the transverse energy 
flow outside of the two jets with the highest Et was studied as a function of for 
the |77*| < 1 and Arj < 1.2. Results for the transverse energy density can be found in 

fig- unq. 
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Figure 100: The corrected transverse energy density in the region 
as a function of x^. The histograms are as in fig. |9| (from j^/j. 



< 1 outside the jets, 



The distribution of the transverse energy around the jet axis was also measured as 
a function of 5(j) (not shown). The jet width obtained in this analysis is similar to the 
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corresponding quantity in the pp collision. 

Results obtained for the jet cross section (here 0.25< y <0.7) are presented in 
fig. |101| a, where da/dEx for the jet production in two t] regions is shown, and in 
fig. p,01| b, where the distributions da/drj for the events with the transverse jet energy 
Et > 7, 11, 15 GeV are presented. The comparison with the PHOJET and PYTHIA 
simulations, with and without multiple interactions, was done for both kinds of distri- 
butions. Note that the rapidity distribution is sensitive rather to the photon structure 
functions, while the Et cross section to the matrix elements for the hard processes. 
Note also that "for Et bigger than 7 GeV previous measurements (HI 93) suffered 
from a defect and are superceded by this new measurement" . 
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Figure 101: a) The cross section da/dEx for the jet production for Et > 7 GeV, in two r] 
ranges: -1 < r] < 2 and -1 < rj < 1. The curves show the Monte Carlo simulations with 
interactions of the beam remnants (full line - PHOJET, dashed - PYTHIA) and without 
them (dotted - PYTHIA). b) Cross section da/dr] versus ij for different thresholds in Et: 7, 
11 and 15 GeV. The curves are as in a); additional dash-dotted curve - PYTHIA with the 
LACl parametrization (from |7^/ ). 



Comment: In addition to the primary hard scattering process, the interaction between 
the two beam remnants is included in the analysis (e.g. using PYTHIA with pip^ > 1.2 
(for GRV LO) or 2.0 GeV (LACl)). It gives '^adequate description of data" for the 
transverse energy versus pseudorapidity and the average energy flow obtained in this 
analysis. 
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"For the first time the underlying event energy has been measured in jet events using 
direct and resolved photon probes. " 

"The strong influence of the underlying event energy on the measured cross section is 
demonstrated. " 

The multiple interaction seems to improve also da/dEx and da/drj distributions 
for jets; within this approach the low and positive rj range is still not properly 
described by the two considered LO parton parametrizations: LACl and GRV within 
the PYTHIA program. PHOJET describes these data. 



•HI 98 J76[ (HERA) 

The two-jet events (the 1994 data) corresponding to to P^ lower than 4 GeV^ and the 
scaled energy y between 0.2 and 0.83 were measured. All jets have E^ > 7.5 GeV, with 
Et ranges above 10 GeV and Ai?^ < 0.25.^^, and pseudorapidity < < 2, with 
At] < 1. The cone algorithm with R=0.7 was used. The PYTHIA 5.7 generator with 
the multiple interaction using the GRV LO parton parametrizations for the proton 
and the photon was used in the analysis (p™=1.2 GeV). The PHOJET 1.06 was also 
used with the same parton parametrizations (with the pr cutoff 2.5 GeV). 

The double differential dijet cross sections, da / dx^/ d\og{E'^/ E^), as a function 



of E'^ for few ranges of are compared with the NLO QCD calculation [ |183| ] and 
PYTHIA (GRV) simulation, see fig. |10^ . 




Figure 102: The double differential cross section as a function of the square of the averaged 
jet transverse energy for the different x.y bins, compared to the PYTHIA (GRV) and the 
NLO calculation (GRV and GS) (from J^j. 
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Figure 103: The effective parton density in tlie plioton (from J5^j. 

The method based on the Single Effective Subprocess Approximation was used 
to extract, for the first time, the effective (LO) parton density of the photon: 

a-%{q^ + ^G^), (25) 

with q'' = Y^{q'' + (p). 

This was done for 63 GeV^ < pi <1000 GeV^, for 0.2 < < 0.4 and 0.4 < < 
0.7. The dependence of the effective parton density on the (= pi) scale is shown 
in fig. |103| and in table ^ (where the statistical and systematical errors are added in 
quadrature). 



Table 24: 











0.3 


112 


3.11 ±0.88 




224 


4.10 ± 1.17 




447 


3.91 ± 1.34 


0.55 


112 


2.25 ±0.55 




224 


3.36 ±0.82 




447 


3.33 ±0.82 




891 


5.18 ± 1.38 



Comment: "Satisfactory overall description (of the double differential cross section for 
jet) except for x^ > 0.6". 

"The effective parton distribution grows with the scale pf, although the increase appears 
slightly steeper than expected from the GRV LO parametrization. " 

•HI 99a QT^ (HERA) 

The new method of extracting the gluon density in the (real) photon from the processes 
with charged particles is introduced. Events collected in 1994 with 0.3 < y < 0.7, 
< 0.01 GeV^, \ri\ < 1 and pt > 2,3 GeV were used. The reconstruction of the 
x-y variables from the highest pt charged tracks follows closely the analysis HI 95a. 



107 



The PYTHIA 5.7 (with and without MI) using the GRV LO parton parametrization 
for the proton, and for the photon : GRV LO, SaSlD, LACl, with the corresponding 
p^*" cutoff 1.2 GeV, 1.0 GeV and 2 GeV. The pt and 1] distributions were measured. 



The result on the LO gluon density at =< pf >=38 GeV^ is presented in fig. ITO^ 



and compared with the jet data based on the 1993 runs at < p^ >=75 GeV^ (see also 

fig- m- 



O HI 1993 (dijet events) 

• HI 1994 (high p^. track evenls) 




logio(xp 

Figure 104: The LO gluon density data extracted from charged particles events at 
< _p2 Qgy2 f-,^ together with predictions of the GRV LO, LACl LO and SaSlD 

LO parametrizations. Also the HI 95a / [7^ data extracted from dijet events at < >=75 
GeV^ are shown (o) (from |7^/). 



•HI 2000a [|78[ (HERA) 

The more precise dijet data taken in 1996 are analysed in regions 0.5 < y < 0.7, 
< 0.01 GeV^ and Et > 4: GeV or Et > 6 GeV. The rapidities lie in the range 
—0.5 < fj < 2.5 and Ar] < 1. In the cross section for Et > 4 GeV the cut on the 
invariant mass of the jets was introduced: Mjj > 12 GeV. The two Monte Carlo 
models were used : PHOJET and "PYTHIA 5.7 + power law pt" (with pt^ = 1.55 
GeV). Multiple interaction was included in both generators. For PYTHIA p™" = 
pjp' = 1.2 GeV, while for PHOJET pr cutoff equals 2.5 GeV. The GRV LO parton 
parametrizations in the photon and the proton were used. The CDFCONE algorithm 
with R=0.7 was assumed for jet reconstruction. 

The measured cross section as a function of for E^ > 4 GeV is presented in 
fig. 105. Similar distribution for Et > 6 GeV (obtained after correcting for pedestal 
energy due to the underlying event) is shown in fig. |106| , where the contribution from 
direct, and resolved - quark and gluon initiated - processes are indicated separately. 
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Figure 105: The distribution of dijet events with Et > 4 GeV (from pB^) 
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Figure 106: The distribution of dijet events with Et > 6 GeV (from |7^/j. 



Using the Single Effective Subprocess Approximation [54| the LO effective parton 
densities of the photon were determined as a function of for < > (=< pf >)=74 
GeV^. The result is presented in fig. |107| (left) and in table |2^. The gluon density 
is determined from the effective parton density by subtracting the quark densities (as 
given by GRV LO in agreement with F2 data). Results of this analysis are presented 
in fig. |107| (right) and in table The obtained gluon density is in agreement with 
that from the single particle data presented in HI 99a. 



109 




Figure 107: The distributions at < >=< pi >=74 GeV^ of the effective parton 
density (left) and the gluon density (right). The comparison with GRV LO (left) and with 
GRV LO, GRS LO, SaSlD, LACl predictions and the HI 99a data (right) is shown (from 



Table 25: 



<Q'> 


< a;^ > 


-^7/7,6// (•^7)/'^ 




[GeV^] 




total{stat.) 


total{stat.) 


74 


0.053 


10.1 ±4.9(3.3) 


4.0 ±2.1(1.4) 




0.094 


6.6 ±2.0(1.1) 


2.4 ±0.9(0.5) 




0.17 


3.5 ±1.2(0.3) 


0.99 ±0.55(0.12) 




0.30 


2.2 ±0.5(0.1) 


0.34 ±0.30(0.03) 




0.50 


2.3 ±0.8(0.2) 


0.33 ±0.40(0.08) 




0.79 


2.1 ±0.4(0.1) 


0.02 ±0.33(0.03) 



Comment: The data reach parton fractional energies down to x^ = 0.05. Leading 
order QCD gives a good description of the Ex > Q GeV data (after subtraction of the 
underlying event energy) which makes possible a determination of the effective parton 
density in the photon. "This quantity is dominated by the gluon density for x^ < 0.2 
which is found to rise strongly towards small x^." 
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•ZEUS 92 []7|] (HERA) 

The evidence for the hard scattering (jet production and resolved photon processes) in 
the photoproduction has been found. The cross section for jets with Ej- > 10 GeV at 
HERA was obtained. The PYTHIA and HERWIG generators, with GS and MT Bl 
parton parametrization for the photon and the proton, respectively, were used (with 
p^*"=1.5 GeV). In the analysis the cone algorithm with R = 1 was applied. 

Comment: "Evidence for the photon remnant jets was found." 



•ZEUS 94 g (HERA) 

The measurement was based on the 1992 data for the single and double jet photo- 
production, with y between 0.2 and 0.7. For tagged events was below 0.02 GeV^, 
otherwise below 4 GeV^. The analysis of the direct and resolved photon processes was 
made using the HERWIG 5.7 and PYTHIA 5.6 generators. In generation of events the 
parton parametrization GRV for the photon and MRSDO for the proton (in addition 
also DG, LAC and MRSD_) were used, with pip'"" =2.5 GeV. The jet-finding cone 
algorithm with R=l was applied. 
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Figure 108: Inclusive jet distributions for (a) transverse energy of jets, (b) pseudorapidity 
of jets. Comparison with the HERWIG prediction for the direct and resolved contributions 
is shown (from JSCj/ ). 



The results for the single jets' concerning the E^ distribution up to E^ = 18 GeV, 



integrated over rapidity t] below 1.6, are presented in fig. |108| a. Fig. |108| b shows the 
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da/drj data where the disagreement with the Monte Carlo prediction occurs for the 
positive 77. 

Dijet production has been studied by selecting events with two or more jets with 
Et>5 GeV, for r] smaller than 1.6. (fig. [T09|) . 

The and Xproton distributions were studied as well for events with lAr^l <1.5, 
I A(/)| >120° and the invariant mass of two jets Mjj larger than 16 GeV (not shown). 
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Figure 109: Kinematic distributions for events with two or more jets: (a) the jet pair 
invariant mass, (b) the transverse energy of jets, (c) the pseudorapidity, (d) cos 6* of jet 
angles in jet-jet CM with respect to the proton momentum for events with Mij > 16 GeV. 
The comparison with the Monte Carlo simulations is shown (from p^j. 



•ZEUS 95a J8T| (HERA) 

The 1993 data for the production of at least one jet with Et > 6 GeV (up to 41 GeV ) 
are presented. Events correspond to: below 4 GeV^, y between 0.2 and 0.85 and the 
jet pseudorapidity between -1 and 2. The PYTHIA 5.6 and HERWIG 5.7 generators 
with the GRV parametrization for the photon and MRSD_ for the proton (and in 
addition LACl and MRSDq) were used. The diffraction contribution was modelled by 
the POMPYT program (p?^^"=3 GeV). 
The cone algorithm with R=l was assumed for jets. 

The transverse energy flow around jet axis was studied. Results in three r] ranges 
are presented in fig. |110| , where "there is some discrepancy for the forward-going jets 
in the 6t] > 1" . 

The Et distributions integrated over two different pseudorapidity ranges and da/drj 
distribution integrated above three Et thresholds: 8, 11 and 17 GeV were measured. 



The results are presented in figs. Ill and 112, respectively, and compared to the 
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PYTHIA prediction (for LACl, GRV HO, ACFGP HO, GS parton parametrizations), 
with p^*" =5 GeV. The discrepancies between the measurement and the LO QCD 
prediction of PYTHIA are restricted to low Ex-, very forward jet. 
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Figure 110: Transverse energy profiles as functions of Srj integrated over |5</)| < 7r/2 (top 
row) and 5(f) (bottom row) (see the text). Results from the PYTHIA simulation (with both 
resolved and direct processes) are shown (from 
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Figure 111: The Et distributions for jets (see text). Results of the PYTHIA simulations 
with the LACl, ACFGP HO, GS HO and GRV HO parametrizations of the parton densities 
in the photon are shown (from 
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Figure 112: The jet pseudorapidity distributions: (left) integrated above three energy 
thresholds, Et > 8, 11 and 17 GeV (see text); the comparison with the PYTHIA sim- 
ulations using the LACl and GRV HO parametrizations for the photon is shown; (right) the 
Et > 8 sample compared in addition with the ACFGP HO and GS HO parametrizations 
for the photon (from Jgl[/j. 



Comment: "In the jet profiles, there is a significant excess of the transverse energy 
density in the data with respect to the Monte Carlo expectations for jets in the region 
1 < ?7 < 2. This excess is located outside of the jet in the forward direction, i.e. 
Ari > 1." 

"Except for the region of very forward, low Et jets, these measurements are fully 
consistent with LO QCD in this new kinematic regime of the structure of the photon" 
. The result (for da/drj with Et > 8 GeV^ and —l<ri<l) "does not support 
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the discrepancy of dcr/dr] with respect to LO QCD calculations observed by the HI 
Collaboration /Hi 93/". 



•ZEUS 95b J82[ (HERA) 

The photoproduction of dijets, with at least two jets of Et larger than 6 GeV , is 
considered in the 1993 data. Events corresponding to the scaled energy y between 0.2 
and 0.8 and lower than 4 GeV^, with median ~ 10^^ GeV^ (for |A?7| < 0.5) were 
grouped in the resolved and direct processes samples. The cone algorithm with R=l 
was used within the HERWIG 5.7 and PYTHIA 5.6 generators with the GRV LO 
parametrization for the photon and the MRSD_ for the proton (with p™"=2.5 GeV). 

The distribution was studied, for more recent results see ZEUS 96, 98a. The 
cut on the x^, equal to 0.75, was introduced later to enhance the resolved or the 
direct photon contributions, and a few distributions were studied separately for these 
samples. 

In fig. |113| the transverse energy flow around the jet axis versus 5ri is shown, for 
the first time separately for the resolved photon and direct photon contributions (with 



the failure to describe low data) 
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Figure 113: The transverse energy How around the jet axis versus Stj: (left) x-y > 0.75 
(direct events); (right) Xy < 0.75 (resolved events). The solid (dashed) line - the PYTHIA 
(HERWIG) simulation (from ^). 



The da/drj was also measured for the direct and resolved photon events (not shown, 
see below for comments and new data in ZEUS 96). 

Comment: "Both simulations fail to describe the transverse energy flow in the for- 
ward region (see also ZEUS 95a and HI 93/" 

"The LO QCD predictions (with DC, GRV and GS2 parton parametrizations) lie 
below the dijet cross section da/dfj data by factor 1.5-2". 

The importance of the NLO calculation is stressed. 

•ZEUS 95c (HERA) 

Photoproduction of three jets is studied for 7p centre of mass energies in the range 
130-270 GeV (1993 data). Events with two high-pr jets (px > 6 GeV, r] < 1.6) were 
selected. A third cluster in the approximate direction of the electron beam is isolated 
using a kx clustering algorithm and identified as the photon remnant. Its properties 



115 



{e.g. transverse and longitudinal energy flows around the axis of the cluster) were 
studied. The Monte Carlo PYTHIA both with default and with harder pt spectrum 
was apphed (the power law with pto = 0.66 GeV). For comparison also HERWIG was 
used. The GRV LO (also LACl) and MRSD_ parton parametrizations were used for 
the photon and proton, respectively. The photon remnant jet was isolated for the first 
time. 



Comment: Properties of the photon remnant jet isolated for the first time "are con- 
sistent with those commonly attributed to jets, and in particular with those found for 
the two jets in these events. The mean value of the photon remnant pt with respect to 
the beam axis is measured to be 2.1 ± 0.2 GeV, which demonstrates substantial mean 
transverse momenta for the photon remnant. " 

"The leading order QCD Monte Carlo simulation, with default parameters, does not re- 
produce the pseudorapidity distribution or the transverse momentum distribution (with 
respect to the incident photon) of the photon remnant. The mean value of pt for the 
photon remnant, 2.1 ±0.2 GeV, is substantially larger than the Monte Garlo expecta- 
tion. Better agreement can be obtained by increasing the average intrinsic transverse 
momenta of the partons in the photon to about 1.7 GeV." 

•ZEUS 96 g (HERA) 

Analysis of the 1994 data for dijets (in events with two or more jets) for Et above 
6 GeV and with the jet pair invariant mass above 23 GeV was performed. Events 
correspond to the range of y between 0.25 and 0.8 and below 4 GeV^. The cone 
algorithm with R=l was used. The PYTHIA 5.7 and HERWIG 5.8 (with the multiple 
interaction) generators (with the MRSA parton parametrization for the proton and 
the GRV LO for the photon) were applied with the pt cutoff 2.5 GeV. 

To obtain the scattering angle cos 6* distribution, sensitive to the parton dynamics 
and not parton densities as in analysis above, the cut was introduced not on At] (as 
in previous analysis) but on f/. 

The results for the uncorrected x^, Xp distributions and for the transverse energy 
flow as a function of St] are shown in fig. |114| . Both the HERWIG (with MI) and 



PYTHIA (without MI) describe the jet profiles data. Due to the cut on f) the absolute 
value of rj is restricted to be below 1.8. Note that the applied cut on the invariant 
mass suppresses events with low x-y. 
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Figure 114: The uncorrected distributions in Xp (a), x-y (b) and of the transverse energy 
How (c,d). The resolved (x^ < 0.75) and the direct (x-y > 0.75) events as a function of Srj 
are presented separately (c and d, respectively). The Monte Carlo (HERWIG) results are 
also shown (from j§^). 

The important results concerning the angular distributions due to various partonic 
subprocesses were obtained for the first time in the large px resolved photon processes. 
In fig. |115| the angular distributions da /d cos 9* for the resolved and direct processes 
together with the LO and NLO calculation based on the CTEQ3M parametrization 
for the proton and the GRV for the photon are presented. The comparison was also 
made with the HERWIG and PYTHIA simulations (not shown). 

Comment: The transverse energy flow is described properly, the "requirements of high 
mass and small boost remove the disagreement in the forward flow between data and 
the simulations which has been reported elsewhere in hard photoproduction at HERA. " 

The dijet angular dependence is well described by the LO and NLO QCD calcula- 
tions, and also by the HERWIG and PYTHIA models. 



117 





30 


"in" 








c 




3 


25 


>. 
k_ 




re 




i_ 


20 


15 




S- 


15 




'© 




tf) 


10 


o 


o 




■D 




e 


5 


■o 










ZEUS 1994 



• 


RESOLVED 




DIRECT 




_ NLO QCD h 




LO QCD ^ 

/ 


L-t — 


, , ,'l , , , 1 , , , 



0.2 



0.4 

Icose'i 



0.6 



Figure 115: The dijet angular distributions for the resolved and the direct contributions 
(see text) (from 



•ZEUS 98a |85| (HERA) 

This is the extension of previous analyses (ZEUS 95b, 96) based also on the 1994 
data on the dijet production with lower than 4 GeV and y between 0.2 and 0.8. 
Results for the production of at least two jets with the pseudorapidity between -1.375 
and 1.875 and for E^^^ = 6, 8, 11 and 15 GeV are presented, assuming |A?7| < 0.5. 
In the data analysis different jet-finding algorithms were applied: the cone algorithms 
EUCELL and PUCELL (both with R=l) and the kr - cluster algorithm KTCLUS. 
The Monte Carlo simulation of HERWIG 5.8 and PYTHIA 5.7 with or without the 
multiple interaction was performed using MRSA and GRV LO parton parametrizations 
for the proton and the photon, respectively. The cutoff for no MI option was p™"=2.5 
GeV, while in the MI approach: =2.5 GeV (HERWIG) and 1.4 GeV (PYTHIA). 
The data were compared to the predictions from a NLO QCD calculation, with the 
additional parameter describing the separation of jets Rsep = R or 2R. 

The resolved cross section was measured in the range 0.3 < a;.^ < 0.75 and the 
direct one - for > 0.75. Analysis of the distribution was performed; jet profiles 
in form of the transverse energy flow around the jet axis and the da/dfj for the various 
jet definitions and transverse energy thresholds were studied. 

The corrected x.y distribution obtained from the two-jet events indicates a need for 
the resolved photon contribution, see also ZEUS 95b and 96. The present analysis 
is based on the KTCLUS algorithm, for results see fig. |116| . The small x-y region is 
not properly described by the Monte Carlo simulations both with and without the 
multiple interaction. 

The transverse energy flow obtained using the KTCLUS algorithm is presented for 
different Ej- in fig. |117| , with a similar discrepancy in comparison to MC (HERWIG, 
with and without MI) seen in the forward direction, as in previous measurements. 
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Figure 116: The corrected x°l'^ distribution. Solid line - HERWIG with the multiple inter- 
action, dashed line - PYTHIA with the multiple interaction, dotted line - HERWIG without 
the multiple interaction (from J§^/). 
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Figure 117: The uncorrected transverse energy Bow around the jet. Solid line - HERWIG 
with the multiple interaction, dashed line - HERWIG without the multiple interaction (from 



The dijet cross section da/dfj obtained under condition \ Ari\ < 0.5, for Ex > 6 GeV 
was measured and the effect of the MI was studied (not shown) . This cross section was 
also studied for different Et thresholds and with different jet definitions, separately 
for > above and below 0.75. The results are plotted in fig. |118| (left) together with 
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the predictions of the NLO QCD approach (the CTEQ3M parton parametrization for 
the proton and GRV HO for the photon; Rsep=2 or 1) [|183|| . 

The same cross section, now with the KTCLUS jet definition, is plotted in fig. |118| 
(right). The predictions of the NLO QCD approach ||183|| with the different parton 
parametrizations for the photon (GS and GRV HO) are compared with the data. 
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Figure 118: Left: Dijet cross sections obtained using different jet algorithms: EUCEL, 
PUCEL and KTCLUS. The curves show the results of the NLO calculation / |J83|/ using GS 
parton parametrization with Rgep = 1 (solid line) and Rsep = 2 (dashed). Right: Dijet 
cross sections obtained using the KTCLUS jet algorithm. The curves show the results of 

1, using different NLO parton 



the NLO calculations K&K UMl and H&O OM] with R 



^sep 



parametrizations: GRV (KSzK - solid, H&O - dotted line) and GS (K&K - dotted line) (from 



"The difference between the parton distributions is largest in the direct photoproduc- 
tion region. This is due to differences between the quark distributions in the photon 
for > 0.8". 

Comment: There is a lack of a good description of the distribution for low x^. 
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The discrepancy in the rj distributions for resolved photon contributions was found for 
events with E™" > 6 GeV. 



•ZEUS 98b J8§ (HERA) 

The inclusive single jet cross section for the below 4 GeV^ with y between 0.2 and 
0.8 and for three energy regions for W between 134 and 277 GeV was studied. The 
measurement of the transverse energy of the jets bigger than 14 GeV and the pseu- 
dorapidity range from -1 to 2 was performed. The Monte Carlo generators HERWIG 
5.8 and PYTHIA 5.7 (default) {pip'=l GeV, and p??^ for two partons equal to 8 GeV), 
including the multiparton interaction, were introduced (using jet cone algorithm with 
R=l, 0.7 and 0.5) in the analysis. Events were generated with the MRSA parton 
parametrization describing the proton structure and the GRV HO and LACl parton 
parametrizations for the photon. 

Transverse energy profiles for jets with Et > 14 GeV are given for R=1.0 and 0.7 
using the PYTHIA (with and without MI), see fig. |119| . The discrepancies are observed 
for the jets with the lowest Et and for rj > 1 for R=l, they are reduced when R = 0.7 
is used. The multiple interaction improves the description of the data for forward jets 
with low Et, but is worse in describing the smaller t] region for R=l, and for whole t] 
region when R=0.7 is used. For jets with Et > 21 GeV no significant discrepancies for 
jet profiles are observed between data and the Monte Carlo simulations (not shown). 
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Figure 119: The uncorrected jet profiles as a function of Sr] (integrated over \5(j)\ < ir/2), 
for Et > 14 GeV, in three r] regions, compared with PYTHIA (with and without MI) for 



R=l and 0.7 (from 



The distributions of jets with E^^^=14, 17, 21 and 25 GeV as a function of t] 



(— 1 < < 2) were measured, see for results fig. 120. They are not properly described 
in the forward low Et region by the NLO calculation |p.83|| if R=l is used (upper 
figures), for R=0.7 the agreement is obtained (lower figures). In the analysis Rsep = R 
or 2R was apphed. 
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Figure 120: The differential cross section da/dr] integrated over Et for four tliresholds: 
Et > 14, 1 7, 21 and 25 GeV. Curves based on tlie NLO calculations using the GRV 
HO (Rsep = R), GS HO (Rsep = R SLnd 2R) and AGF (Rsep = R) parton parametrizations 
for the photon and the CTEQ4M for the proton with R=l (upper figures), R=0.7 (lower 
figures) are shown (from JS^j. 



The r/ distributions in three regions of energy W are not in agreement with the 
NLO calculation for the R=l, as can be seen in fig. |121| . The data for R = 0.7 (not 



shown) are in agreement with a QCD calculation, the same as in the previous analysis, 
ZEUS 98a. 
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Figure 121: The differential cross section da/drj for three regions of the energy W. The 
PYTHIA results with the MRSA for the proton and with different parton distributions in 
the photon (with Rsep = R or 2R) are shown (from i8t]). 



The cone radius dependence of the cross section was studied as well to elucidate 
effects of a possible underlying event, see fig. |122. 
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Figure 122: The cross section dependence on the R parameter for Et > 21 GeV and -0.5 
< rj < 2 obtained using various factorization scale Q (= //) in the LO and NLO calculations 
(with the CTEQ4 and GS96 parton parametrization for the proton and photon, respectively) 
m (from iH/j. 
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Comment: "...the uncertainties on the jet measurements due to possible underlying 
event contributions become reduced at high (Ej- > 21 GeV) or when using a reduced 
cone radius (R = 0.7)." 

"The measured cross sections for jets with R=0.7 are well described by the [NLO 
QCD] calculations in the entire measured range of rj and Ej-." 



•ZEUS 98c,conf §^ (HERA) 

The Et measurement for jet photoproduction (P^ < 4 GeV^) at 134 <W < 277 GeV 
is presented. Data collected in 1995-1997 lie in the large Et range: 17-74 GeV and 
-0.75 < 7] < 2.5. Two methods of jet identification were used - with the iterative cone 
algorithm (R=0.7 and 1) and with the fcr cluster algorithm. Comparison was made 
with the NLO calculation (with Rsep=R), based on the CTEQ4M and the GRV HO, 
GS 96 and AGF parton parametrization for the proton and photon, respectively. 



The Et distributions are presented in fig. |123| for various jet definitions: cone 
algorithm with R=l (left), cone algorithm with R=0.7 (center) and fcy cluster method 
(right). Fig. p.24| shows the comparison of different ratios of da/dET- The cone 
algorithm with R=l and the fcr cluster approach give similar results, larger than 
results obtained with R=0.7. 
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Figure 123: The Et distributions obtained using cone algorithm (left with R=l and center 
with R=0.7) and the Ut cluster method (right), compared to the NLO calculation /[i83|/ (from 
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Figure 124: The ratios of the Et distributions obtained with the jet definitions, as used in 
fig. p3 (from JSf j. 
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The comparison with the theoretical predictions using the CTEQ4M parton parametriza- 
tion for the proton and the GRV HO, GS96 and AGF densities for the photon, with 

Rft 



''Sep 



R, and different jet definitions is presented in fig. |125 
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Figure 125: The fractional differences between the measured da/dEx for jets and the NLO 
calculations (using the GRV HO parametrizations - black points ) with the cone algorithm 
(left with R=l, center with R=0.7) and the cluster method (right). The comparison with 
the NLO calculations with different parton parametrization for the photon: GS96 (dashed 
lines) and AGF (dash-dotted lines) is also shown (from IpDJ). 



Comment: "The NLO calculations give a reasonable description of the measured dif- 
ferential cross section in magnitude and shape. " 

•ZEUS 99a JS§ (HERA) 

The 1995 measurement of the photoproduction of dijets for < 1 GeV ^, and y 
between 0.2 and 0.85, corresponding to 134 < W < 277 GeV was performed, also 
a narrower y range 0.50 to 0.85 was studied. The aim was to constrain the parton 
densities in photon at high x^. Events were separated in to the direct sample [x^ > 
0.75) and the resolved photon sample (x^ < 0.75). The asymmetric cut is applied on 
the two highest transverse energy jets. The cross section is symmetrized in the pseu- 
dorapidities of these jets. Events correspond to the highest jets with the threshold 
Erieading = 14 GeV and for the second jet the threshold is given by ETsecond = 11 
GeV . The rapidity range —1 < rj < 2 is considered. Two Monte Carlo generators 
were used: HERWIG 5.9 and PYTHIA 5.7 with the fc^ clustering method for the jet 
identification. 

The uncorrected x^ distribution is shown in fig. |126 , together with Monte Carlo 



predictions for a direct and the full (direct+resolved) contributions. 
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Figure 126: The uncorrected distribution for the direct (defined as x-y > 0.75) and the full 
(direct+resolved) contributions obtained in the PYTHIA 5.7 and HERWIG 5.9 generators 
for Et > 14 GeV (from m. 



The energy flow around the jet was studied for different regions of and for 



different ET,ieading thresholds, see fig. 127 
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Figure 127: The transverse energy Bow around the jet axis (integrated over 
different ranges of transverse energy and x-y (from 



< 1), in 



Below the results for the Et distribution in different regions of 771,2 are presented 
(fig. |12(^ ), with the direct contributions displayed separately. Comparison with the 
NLO calculation (with CTEQ4M and the GRV HO, GS96 HO, AGF HO parton 
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parametrizations for the proton and photon, respectively) is shown. 
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Figure 128: Left: Dijet cross section da /dETieading with one jet within the rapidity range 
1 < r]i < 2 and the other jet in three rapidity ranges, described in the figure. Right: Dijet 
cross section da/dEx with one jet within rapidity range < rji < 1 and the other jet in 
three rapidity ranges, described in the figure. Comparison with the NLO calculation (with 
GRV HO) JJg3[ / - thick line (thin line shows the direct contribution) (from 



For the leading jet, with ETieading > 14 GeV, the rapidity distribution da/drj2 was 
measured as well. It was done for a full y region measured in this experiment and, 
to enhance the sensitivity to the choice of parton density in the photon, for a large 
y region (i.e. 0.5 < y < 0.85 corresponding to 212 < W < 277 GeV^). Three rji jet 
rapidity ranges were studied, see fig. p.29| for results. Here the comparison was made 
with the NLO calculation [|18D with the GRV HO, AGF HO and GS96 HO parton 
parametrization for the photon (the CTEQ4M parametrization for the proton). Also 
the comparison between the different NLO calculations ( ||176| , |188| , |185| , |186| , |187|| ) is 
made using the particular (GRV-HO) parton parametrization. 
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Figure 129: The rapidity distribution da/dr]2 for 0.2 < y < 0.85 (upper set of figures) and 
0.5 < y < 0.85 (lower set of figures). The extra cut for the leading jet, Exieading > ^4 GeV, 
was applied. Three r]i jet rapidity ranges were considered. The comparison with the NLO 
calculation with different parton parametrization for the photon is made in a), b) and c). In 
d) comparison of different NLO calculations with the particular parton parametrization for 
the photon (GRV-HO) is shown, see text (from Jg^/ j. 



Comment: The underlying events "play no role in the present kinematic regime. " 
For the full and for the higher energy range NLO QCD calculations using the GRV 



128 



HO, AGF HO and GS96 HO parametrizations for the parton densities in the photon 
can describe properly the distribution in ETiaeding- However "at central and forward 
pseudorapidities , both for the full and for the high range, the data lie above the 
NLO calculations" . It "suggests that in this kinematic region the parton densities in 
the photon are too small in the available parametrizations" . "This region has not been 
studied in F2 measurements. " 



•ZEUS 99b,conf (HERA) 

This is the continuation of the previous analysis of the dijet events for Et up to 75 
GeV. The data were collected in 1996-1997 with the kinematical cuts as in ZEUS 
99a. The fcr clustering algorithm was used for jet identification. The thresholds on 
ETieading froui 14 to 55 GeV was introduced. 

The x-^ distributions for different transverse energy thresholds are presented in 



fig. |130| ; they are well described by Monte Carlo generators (also the transverse energy 
flows, not shown). 
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; 10000 




Figure 130: The x-y distribution (with the statistical errors only) for different thresholds on 
the highest transverse energy of the leading jet. The second one has ETsecond > 11 GeV. 
The solid histograms represent predictions of HERWIG 5.9, the shaded ones represent the 
direct contributions (from JS^j. 



Both the transverse energy distribution for Exieading and the rapidity distribution 
were studied. Below, the results for the Exieading distribution in different regions of 171,2 



are presented (fig. |131|) . The direct (defined as x-y > 0.75) contributions are displayed 
separately. The results for the rapidity distributions for various thresholds on ETieading 
from 14 to 29 GeV are shown in fig. |131| , and in fig. |132| for a narrower high y region. 
The data are compared with the NLO calculation using CTEQ4M (proton) and AGF 
HO (photon) parton parametrizations, and some discrepancies with data are found. 



129 




Figure 131: Left: Dijet cross section da /dExieading with one jet within the rapidity range 
1 < r]i < 2 and the other jet in three rapidity ranges, described in the figure. Right: Dijet 
cross section da /dETieading with one jet within rapidity range < rji < 1 and the other 
jet in three rapidity ranges, described in the figure. Comparison with the NLO calculation 
|J8t^/ with the AGF HO parametrization for the photon and the CTEQ4M parametrization 
for the proton - solid line (dashed line shows the direct contribution) (from JS^j. 



Comment: These results confirm those in ZEUS 99a. 

"...the measured cross sections remain higher than the NLO Q CD predictions when the 
transverse energy threshold on the leading jet is raised up to 29 GeV" suggesting "that, 
in the kinematic region of the measurement presented here [high and high Et], the 
available parametrizations of the parton densities in the photon are too small. " 
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Figure 132: The rapidity distribution da/d'r]2 for 0.2 < y < 0.85. Five different cuts for the 
leading jet ETieading were applied from 14 to 29 GeV. Three rji jet rapidity ranges were con- 
sidered in a), b) and c), with the AGF HO parametrization for the photon and the CTEQ4M 
parametrization for the proton - solid line (dashed line shows the direct contribution). The 
particular parton parametrization for the photon (AGF HO) was used (from 
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Figure 133: The rapidity distribution dcr/d7]2 for 0.5 < y < 0.85, see Gg. \13^ (from 



•ZEUS 98d,conf g (HERA) 

The high mass {Mjj > 47 GeV, up to 140 GeV) dijet cross section was measured 
for 134 < < 277 GeV (the 1995-1997 sample). Events with smaller than 4 
GeV^ were collected and compared with the same NLO calculation ||183|| as in ZEUS 
99b,conf. In the present analysis the cone algorithm with R=l was used as well as 
algorithm. 

Results for the cos^* distribution obtained within these two-jet algorithms are 
presented in fig. p. 34] . In the fig. |135| the invariant mass distribution (integrated over 
the I cos^*| < 0.8) is shown. 
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Figure 134: The cos 9* distribution for events with Mjj above 47 GeV obtained using two- 
jet algorithms: the cone algorithm (left), the clustering (right). The comparison with 
the NLO calculation using the GRV HO parton densities (with Rsep = 1.4i? for the cone 
algorithm) and GS96 (Rsep = R for the cone algorithm) is made (from |9^/ }. 
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Figure 135: The Mjj distribution for events with \cosO*\ < 0.8 obtained using two-jet 
algorithms: the cone algorithm (left), the kx clustering (right). The comparison with the 
NLO calculation as in fig. \134 (from j^). 



Comment: "NLO QCD calculations account reasonably well for the shape and magni- 
tude of the measured da /d cos 6* and da/dMjj". 
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2.4.5 Prompt photon production in the resolved 7p colHsion 



The prompt photons, i.e. large-pr photons (predominantly direct), produced in the jp 
collision give an additional information as compared to jet(s) production discussed 
in the previous sections. Although the corresponding rate is much smaller (~ a^), 
the advantage is the possible direct access to the hard subprocess. However, there 
are also indirect contributions due to the quark and gluon fragmentation into photon 



189 1 . In practice the isolated photon production is considered. It corresponds to the 



restriction of the hadronic production around the photon direction. Usually the cone 
algorithm, with R defined as for the jets, is used with a limit on the allowed hadronic 
energy {e.g. 

EE^ < O.IE^) within the cone |[L9CI|| . 



The photoproduction of a hard direct photon may occur due to the direct interac- 
tion of the initial 7, or due to the resolved initial 7. The sensitivity to the particular 
subprocesses and to the form of parton density in the photon differs in such processes 
from that in the jet production, discussed in sees. |2.4.3| and |2.4.4 The domination 
of one contribution due to the gluonic content of the initial photon leading to the 
production of the forward, very energetic photons was pointed out in ||189|| g. This can 
be particularly useful in extraction of the gluonic content of the real (and also slightly 
virtual photon) ||191|| . Also there is a chance to test a special photon-gluon fusion 



process, namely the 'jg 'jg box diagram contribution. 

All these make prompt photon processes a source of valuable independent infor- 
mation on the photon-hadron interaction at high energies. New data which appeared 
recently (ZEUS 2000a) can be compared to the two versions of the NLO calculations 
|192|| . The agreement of the MC models and the NLO calculation with data is rea- 



sonably good for rapidities of the final photon in the range 0.1 < r] < 0.9, while for 
—0.7 < 7] < 0.1 the excess of the data over the predictions is seen. 

DATA 

•HI 97a,conf 0] (HERA) 



The measurement of the high Ej- photons {Ej^ > 5 GeV and —1.2 < t]^ < 1.6) 
in quasi-real jp collisions was performed, based on 1996 data. The isolation of the 
photon within a cone with R=0.8 around the photon direction was imposed. The total 
hadronic energy in this cone was restricted to be not higher than 10% of the photon 
energy. The balancing hadronic jet with transverse energy Et was also measured. 
The PYTHIA generator with the GRV LO parton densities in the proton and in the 
photon was used. The results for the transverse energy and for the pseudorapidity of 



the final photon distributions are presented in fig. |13(:. 



^^called also the deep inelastic Compton (DIG) process, as it corresponds to the jq jq subprocess 
in the Born approximation 
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Figure 136: Transverse energy Ej, (a) and the pseudorapidity rp (b) distribution of the 
produced photon (from /p3/j. 



The number of events as a function of the azimuthal angle between the final photon 
and the jet (A0(^jet)) and the transverse energy of the accompanying jet are presented 
in fig. |137| . 
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Figure 137: a) The angle between the final photon and the jet, A(/)(^ jg^-j, distribution, b) 
the transverse energy of the accompanying jet distribution (from ^/ j. 




The measurement of the photon and the jet allows to reconstruct the and Xp 
variables. Their distributions are shown in fig. p.38| . 
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Figure 138: The (a) and Xp (b) distributions compared with the MC predictions (from 



Comment: "It is shown that the prompt photon process at HERA is dominated by the 
direct '-yp interaction as expected from QCD. " 



•ZEUS 97 g (HERA) 

The photoproduction of the isolated prompt photon with high Ej^, together with a 
balancing jet (described by E^ and rj) has been observed for the first time at HERA. 
The data (collected in 1995) correspond to < i GeV^, 5 < < 10 GeV and 
Et>5 GeV, and in addition: -0.7 < 77"^ < 0.8 and -1.5 < 77 < 1.8. 

The isolation cone was imposed around the photon candidate; within the cone with 
R=l the hadronic energy was assumed to be smaller than 10% of the photon energy. 
The Monte Carlo PYTHIA program with the MRS A (proton) and GRV LO (photon) 
parton densities was applied. Results are presented in fig. [L39| for A0(7, jet), and 
AET{'y,jet) distributions. In fig. |140| the and Xp distributions are shown. 
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Figure 139: The A<j){'y,jet) (a) and AEri'yJet) (b) distributions (from Mj). 
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Figure 140: The and Xp distributions are shown for the isolated 7 plus jet production 
(from 



Comment: For the first time the isolated high - Et photons accompanied by balancing 
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jets were observed at HERA. The distribution is in good agreement with the LO 
QCD prediction of the PYTHIA generator. 



•ZEUS 2000a Q (HERA) 

The first inclusive measurement of pliotoproduction of isolated prompt photon ob- 
tained from run 1996-1997 (P^ < IGeV^ and 0.2 < y < 0.9) is reported. Events corre- 
spond to the following ranges of kinematical variables: Ej^ > 5 GeV and —0.7 < t]"' < 
0.9 for the photon. The PYTHIA 5.7 and HERWIG 5.9 generators, with p^" = 2.5 
GeV and no implementation of multi-parton interactions, were used. The proton and 
photon parton parametrizations, MRS A and GRV LO, respectively, were used in the 
analysis. 

The isolation cut for the photon was imposed: in the cone R=l the total hadronic 
energy was restricted to be below 0.1 Ej.. The Ej. distribution is presented in Fig. 
|141| together with the MC predictions and the NLO calculations by Gordon and by 
Krawczyk & Zembrzuski ||192|| . In Fig. |142| the data for the Vj"^ distribution together 



with the predictions of the MC models and of the NLO calculations ||192|| are shown. 
In order to find an origin of the discrepancy between data and both MC and NLO 
calculations the rapidity distributions were studied in three y bins (not shown). The 
disagreement is strongest for the smallest y bin {y between 0.2 and 0.32) 
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Figure 141: The Ej, distributions for -0.7< irf <0.9 and 0.2< y <0.9, is presented together 
with the MC predictions and the NLO calculations by Gordon and by Krawczyk&Zenibrzuski 
|J9^/ with the GRV HO parton parametrizations (from |9^/). 
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Figure 142: The -rf distribution for 5< <10 GeV and 0.2< y <0.9, together a) with 
the predictions of the PYTHIA and HERWIG models, b) with the NLO calculations by 
Gordon (LG) and by Krawczyk&Zembrzuski (K&Z) g9|/ with the GRV HO and GS parton 
parametrizations for the photon, (from 



Comment: "The [MC] models are able to describe the data well for the forward (proton 
direction) photon pseudorapidities, but are low in the rear direction. The disagreement 
is strongest in the W interval 134-110 GeV. This result, together with the disagreements 
with the NLO predictions seen also in recent dijet results at HERA /ZEUS 98sl], would 
appear to indicate a need to review the present theoretical modelling of the parton 
structure of the photon. " 
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3 Partonic content of the virtual photon 



The notion of partonic content of the virtual photon has appeared in high energy 
interactions soon after the related concept for the real one (@a, [^). From the very 
beginning it has a status of a unique test of QCD, as without additional experimental 
or model assumption the definite, singularity free predictions can be derived for both 
the X and dependence with the LO and NLO accuracy]^. By virtue of the analytical 
continuation (dispersion relation) one expects that a limit of zero virtuality will lead 
to the partonic content of the real photon. 

Not only the DlSe^* can be performed for the virtual photon, where the correspond- 
ing structure functions of the 7* are measured, but also the resolved virtual photon 
processes sensitive to the individual parton densities in 7* can be studied. 

The deep inelastic scattering on the virtual photon with hadronic final state was 
studied experimentally already in 1984 (PLUTO 84). Recently new data on DISe7* 
have appeared from the LEP accelerator (L3 2000). Also, the measurement of the 
leptonic final state in the deep inelastic scattering on the virtual 7* has been performed 
at LEP (OPAL 99b). It shows the importance of the interference terms which makes 
the extraction of the corresponding structure functions ^^^^^ unfeasible in some 
kinematical regions (see sec. ^). 

The bulk of the resolved virtual photon data, with up to 80 GeV^ is just arriving 



from the HERA collider (single and double jet production), see sec. |3.3| . The early 
analyses treated these events only as a special class of DISep events - therefore the 
standard (for such analyses) types of Monte Carlo generators were used. In the recent 
analyses the generators incorporating resolved virtual photon interaction are used. The 
recent measurement of jet production in the virtual photon-proton scattering allowed 
to extract the effective parton density in 'y* (HI 2000b). Still some discrepancies are 
observed, e.g. in ZEUS 2000b, where the dijet production in the 7*p collision for 
below 4.5 GeV^ cannot be properly described by Monte Carlo models nor by the 
NLO QCD. 

A dedicated study of the forward production of particles and jets, performed at 
HERA with the aim to establish a dynamics of the parton radiation at small xbj 
uses a large sample of the Monte Carlo generators - those typical for DlSgp events 
and those for processes with a resolved virtual photon. Here the possible signal from 
the BFKL evolution for the partons from a proton is confronted with the importance 
of the contribution due to the resolved 7* interaction. We discuss this topic in the 
separate section (sec. |3l^ ). 



3.1 The virtual photon- induced processes 

The description of hadronic high energy processes with the initial virtual photon de- 
pends on its characteristic virtuality: on how much it differs from zero and how big 
it is as compared to other relevant scales, e.g. Aqcd ot pt- From these relations 
one can derive the role of the virtual photon in the considered process. However, the 
interpretation is highly nontrivial, as it may depend on the reference frame, and even 
on the chosen gauge. We will not discuss these theoretical aspects any further. Below 



^^However, the importance of nonperturbative components of the virtual photon was pointed out, 



see below (sec. 3.2.1). 
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we try to underline only the roles in which the photon appears in the experimental 
analyses. 

We start with the short discussion of the flux of virtual photons. 



3.1.1 The flux of virtual photons 

As for today, we always deal with virtual photons arising from the leptonic beam(s)0. 
Assuming the factorization holds between the leptonic source and the hadronic hard 
subprocess, one can introduce the flux of the transverselv polarized virtual photons 
(e.g. ia, H) 

fr/eiy^P ) = ^[ ]■ (26) 

The flux of longitudinal virtual photons , 

is expected to be small for the present experimental setup and it is usually neglected 
in the analyses [ 198 , 27 1. 



Both the flux of virtual photons emitted by the electron (eqs. |2g, ^), and the cross 
sections involving initial virtual photons {e.g. a'^*^ and a^*^) depend on convention 
used in the definitions^. This ambiguity is however absent in the ee or ep cross 
sections, where a convolution of the virtual photon flux and the virtual photon cross 
section appears {e.g. a^^ ~ /e®0"'^'^). Note that the concept of the resolved electron 
may happen to be very useful here, see sec. E. 



The virtual photons described by eqs. (|26| , |27|) undergo further interaction. They 



may interact directly, or via their partonic agents, provided that in a process that 
follows a probe exists, hard enough to resolve the virtual photon. 



3.1.2 The virtual photon as a (non-elementary) probe 

The description of the inclusive processes with 7* as a probe of another object, e.g. 
of a real photon in DISe-y, is relatively simple if global quantities, like total 7*7 cross 
sections or F2{xBj,Q'^) for a real photon, are considered. Here the virtuality of the 
photon probe (we denote it consequently through the survey by Q^) provides the 
hardest scale in the process, and a direct interaction of 7* with partonic content of 
the target dominates. However, when one tries to reconstruct the different classes of 
hadronic final states in the DIS-type of measurements, among them the rare events 
with Px ^ Q^, one should include the interaction due to the partonic content of the 
probe - 7* (see also discussion in sec. p.3| ). 

Similarly, in the case of ep collision the highly virtual photon, being a probe for 
the proton target (DISep events with standard variables and xbj), can be resolved 
itself by the large pt jet (particle) production^ (see below). 

^^Protons, due to their larger mass, are much weaker source of the virtual photons, 
^^as they contain the convention - dependent flirx of the initial virtual photons. 
^"This problem ("Who is probing who?") was underlined e.g. in ]9^ . 
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3.1.3 The virtual photon as a target 



Sec. ^ covers the dedicated measurements of the "structure" of the virtual photon. 
Such a photon, with virtuahty P^, plays here a role of a target - with the structure 
to be resolved. The structure functions for such object are measured in the DlSg^* 
experiments at e~^e~ colliders in the double-tag events, with ^ P^. Jet production 
resolving the initial virtual photon, i.e. with p"^ ^ P^, corresponds in e~^e~ colliders 
to the single (7*7) and double-tag (7*7*) events, and in the ep case - to the single 
tagged events 7*p. In order to pin down the partonic content of the virtual photon 
also large pt particles (among them a real photon ||191|| ) can be used. 

In practice an admixture of the resolved virtual photon contribution is always 
present in measurements on the "real" photon discussed in the previous section. Data 
are usually corrected for this effect (see sec. H). 

As it was mentioned above, there is a possibility to introduce here the structure 
function of the electron (see sec. P), given as a convolution of the flux for virtual 
photons (eqs. |2^, |2^) and the structure function of 7*. 



3.2 DISe7* experiments 



Measurements of the virtual photon structure functions in the deep inelastic scattering 
at e'^e^ colliders are performed using double-tag events. These events are selected in 
the kinematic region where one of the virtual photons (the probe) has, on the average, 
a large virtuality and the other, the target, a small one, P^ <C Q^. 

For general double-tag events the hadronic state X is produced in the lepton beams 
collision, e{pi)e{p2) e{p[)e{p'2)X , via the 7(g)7(p) collision with definite helicity 
states of the photons. The corresponding cross section for the unpolarized lepton 
beams is given by (see 0a,b,c) 
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where T denotes transverse (+ or -) and L - longitudinal helicity states (0). The q 
and p stand for the photons four- momenta: q = pi — p[, p = p2 — p^- The pi and 
P2 are the photon (with q and p four-momentum, respectively) density matrices. The 
o'tt,tl,lt,ll and ttt,tl denote the corresponding cross sections and interference terms 
(the first subscript corresponds to the photon with four- momentum q). The is the 
angle between two scattering planes of the scattered electrons in the 77 CM system. 

For large virtualities of both photons, = = 3> 4mg, eq. ^ can be 

written in the following form: 



, , da{ee eeX) 



E[E: 



Ltt(o'TT + eiO"LT + ^20'TL + ^1^20'LL + 



d^p'ld^p'2 

+ ^eie2i"LL COS 20 - 2^JeJYT~e{)^Je^(lTe2)TTL coS( 
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where 



a 



-'TT 



1/2 



and 



For ei close to 1, which is a typical value in present experiments, one can furthermore 
simplify the above formulae, and obtain: 

^, ^, daiee — > eeX) 1 - - 

E1E2 ^3^;^3^,^ = LTT{(yeff + -TllCOs2<P - 4rTLCOS0), 

where an effective cross section o"e// = o"rr + ctlt + o"rL + o"Lr was introduced. 

For the deep inelastic scattering, where > P^, the corresponding structure func- 
tions for a polarization-averaged photon target with virtuality can be introduced 



(29) 



47r2a /(g . p)2 _ g2p2 



with 



F2 =^2-^ 



and X = Q^/2pg. 



(30) 



(31) 



3.2.1 Theoretical description 

The structure function of the virtual photon can be obtained in the Parton Model 
assuming the production of the gg pairs in the 7*(Q^)7*(-P^) collision. In the mass 
range 



> > ml 



(32) 



it has the form [^: 



Ff (x, Q\ = <Q^> -x{[x^ + (1 - xY] In + 6x(l - x) - 2}, 



where 



1 

i=i 



(33) 



to be compared with the eq. dTHI). 

One can clearly see that the scale of the probe, Q^, has to differ from the P^ in order to 
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test the structure of the virtual photon in the leading logarithmic accuracy (LL). The 
corresponding quark density in the virtual photon defined in the LL approximation 
has the form: 



ILL 
\PM- 



Ztt 



2 r 2 I 



XI 



p2 



The QCD evolution equations for the virtual photon are analogous to those for the 
real photon with the inhomogeneous term given by the corresponding splitting function 
In the case of the virtual photon there is a hope that to solve evolution equations 



the initial conditions are not needed, since for ^qcd the nonperturbative 

effects should be absent (see ref.[^). Indeed, the virtual photon may play a unique 
role in testing the QCD. 

However, recently there appeared papers [p.62| , |17(]| , |163|| showing that the parton 
content of the virtual photon is not solely described by purely perturbative contribu- 
tions in this region, in contrast to the expectation from p2[ . 

Recently the LO ^7^, [l98l and the NLO QCD pl, |195|, |19|, |193 calcula- 
tions have appeared for the DISe7* ^^e jet production in the considered processes 
involving virtual photon(s) 7*7, 7*7* and 7*p. 

The existing parton parametrizations for the virtual photon are described in the 
Appendix. 



3.2.2 Measurements of F2 

For the DIS^* measurements, the quantity 

x,is = QViQ' + P' + w^,,) 

needs to be converted to the true xbj = Q'^/'^PQ variable 0. 

We start with the old data from PLUTO, the first "DIS"-type measurements for 
the virtual photon. 



■^^At finite a modified variable xbj, wliicli extends over tfie whole range between and 1 may 
be introduced. For small ratio P^/Q^ ^ 1 they coincide. 
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DATA 



•PLUTO 84 [§ (PETRA) 

The double-tag events were measured in e~^e~ collision, with the energy of the beam 
equal to 17.3 GeV. One of the virtual photons (the probe) had, on the average, virtu- 
ality < >=5 GeV^ and the other (the target) < >=0.35 GeV^. For the cross 
sections, the assumptions ari ~ c"lt and aLi ~ were made. The experiment was 
sensitive to the following combination of the virtual photon structure functions (the 
effective structure function): 

F,jf = F2 + {3/2)FL. 



Results for the extracted Fg// as a function of xbj are presented in fig. |143| together 
with the theoretical predictions. In fig. |144| the quantity Q^a^y/Air'^a'^ ^ Feff/a 
eraged over both xbj and Q^, is shown as a function of the measured (0.2 ■ 
GeV^), in comparison with QPM and VMD (for the target) contributions. 



av- 
0.8 



a 



Q6- 



PLUTO 



QPM lm:0) 



HOQCD 




Figure 143: The data for the effective structure function F^ff/a for the virtual photon. 
For the averaged = 5 GeV^, and = 0.35 GeV^ , the dependence on xbj is shown in 
comparison with the QPM (with massless and massive quarks) and QCD calculations, both 
with the VMD for the target contributions added (from ^). 
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Fi gure 144: The results for the -^^-^cr^^ for virtual photon as a function of P'^ (averaged 
over XBj and Cf'), compared with QPM and VMD contributions (from ^). 



•L3 2000 g2|] (LEP 1) 

The measurement of the structure function for the virtual photon (< P^ >= 3.7 
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GeV^) at \fs= 89-92 GeV was performed based on the 1991-1995 data. The large 
Q'^ region, between 40 and 500 GeV^ (< >=120 GeV^), was analysed. The real 
photon structure function was studied as well, see sec. |2.2.2| . 
In the analysis the effective structure function: 

-^e// = 4^('^TT + + C^TL + LL) 

was measured. 

The data were analysed using the JAMVG generator modelling the QPM with 
Nf=A, the PHOJET 1.05c with a cutoff p?^^"=2. 5 GeV and TWOGAM generating the 
three processes: QPM, VMD and QCD resolved photon contributions. The hadronic 
final state was investigated (see section |2.3| for details). 

The F^fj data are presented in table ^ and fig. |145| , where the results are compared 
with the predictions of the QPM and QCD calculations performed for the structure 
function of transverse photon only. The QCD contribution is needed to describe the 
XBj region below 0.5. 



Table 26: 



<Q2> 


XBj 


{stat. + syst.) 


120 


0.05-0.20 
0.20 - 0.40 
0.40 - 0.60 
0.60-0.80 
0.80-0.98 


0.42 ±0.16 ±0.05 
0.71 ±0.24 ±0.09 
0.72 ±0.34 ±0.09 
1.27 ±0.51 ±0.16 
1.48 ±0.66 ±0.19 



LL 











Q^ = 120 GeV^ 
• DATA 
— QPM 
-- GRS 



= 3.7 GeV^ 



L3 



0.4 0.6 



0.8 



Figure 145: F^^j^/a as a function of XBj for < >=i20 GeV^ and < >=3.7 QeV^. 
The solid line corresponds to the QPM prediction. The QCD prediction (only for transverse 
target photon states) using GRS LO parton paranietrization is given by the dashed line, 
(from H/j. 



The averaged values of F^j-^/a as a function of < > are presented in table ^ 
and in fig. |146 
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Table 27: 



< p2 > 

[GeV^] 


< F'^fJa > 
(stat. + syst.) 
= 0.05 - 0.98 





0.83 ±0.06 ±0.08 


2.0 


0.87 ±0.25 ±0.11 


3.9 


1.00 ±0.32 ±0.13 


6.4 


1.02 ±0.70 ±0.13 



LL 

V 
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Figure 146: The dependence of Fjjj/a averaged overxsj = 0.05—0.98 for the double-tag 
data, compared with the QPM prediction (sohd hne) (from p^/j. 



3.3 Measurements of resolved virtual photon(s) processes 

As in the case of the real photon, the large pt jets or particles may resolve the virtual 
photon(s). Provided the corresponding mass relation ~ pf. ^ (P2^) ^ ^qcd 
holds, one introduces the contribution due to a partonic content of 7* and a direct 
7* subprocesses, as for the real photon (eqs. (|20| , |2^) ). For the same reason as in sec. 
H we will limit ourselves to the jet(s) production. The exception is the hard forward 



production, where a single particle production is also considered, see sec. ^A . 

As we already discussed in sec. |2]^, in analyses of hadronic final state, accompany- 
ing the F2 measurements, jet production in 7*7 collision was studied, e.g. TOPAZ 
94, DELPHI 96b,conf, OPAL 97d,conf. In some of these analyses, e.g. TOPAZ 
94, there were events with the hard scale ~ p"^ much larger than the virtuality of 
7*, so 7* can be considered as being resolved. 

However, so far there are no data from dedicated experiments on the resolved 
virtual photon processes in e+e~ collisions at LEP and TRISTAN, whereas a large 
amount of data from the ep collider HERA have appeared recently. The single and 
double jets produced in the resolved virtual photon processes are being studied in 
detail at HERA. Three and four jets including remnant jets, and jet shapes were also 
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studied, however we will omit these results in our presentation. For the first time the 
effective parton density, based on polarization-averaged parton densities in the virtual 
photon, was extracted from the data (HI 2000b). 

The transition region between the interaction of an almost real photon and of a 
virtual photon with the proton is studied by the HI and ZEUS collaborations in the 
ep collisions at HERA. In such analysis various Monte Carlo generators (LEPTO, 
ARIADNE, MEPS, DISJET, PROJET, RAPGAP and PHOJET, HERWIG), used to 
describe the photoproduction and the DISep events as well as the rapidity gap events 
at HERA, are used. 

There are some discrepancies observed between the virtual photon cross sections 
for the dijets production at HERA and the NLO QCD |194|, |19|, |9|, [Tg^] calculations, 
see ZEUS 2000b. 

3.3.1 Jet production in 7*7 and 7*7* collisions 

The jet production in 7*7 collision were studied in DlSe^* experiments, e.g. ALEPH 



97a,conf, TOPAZ 94, see sec. p.Sj There are no data so far for the jet production 
in the double-tag events. 



3.3.2 Jet production in 7*^ collisions 



In the experimental analyses at HERA the flux of virtual photons is introduced (HI 
97b) in the form of the integral of the transversely polarized virtual photons over the 
relevant range of y, and over its positive squared mass (virtuality), P^: 



2 

^ - dy dP'f^./MP')- (34) 



' ymin m 



The kinematical variables used in the analysis of the single or double jet production 
are defined as for the real photon case (sees. p.4.2| and p.4.4|) . In the resolved virtual 



photon processes some kinematical variables are defined also in the 7*p CM system, 
and are denoted below by a star, e.g. E^. Note that here we use the notation 
for the virtuality of the photon although in the context of the DlSgp events (on the 
proton) at HERA it plays the role of the Q^. The hard scale is usually provided 
by the transverse energy or transverse momentum of jets. 
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DATA 



•HI 94 (HERA) 

The jet rate measurement in DISep events with ^ was reported. The JADE jet 
clustering algorithm was used. 



•HI 97b g (HERA) 

The single jet cross section for the events with 0.3 < y < 0.6 was studied in the 
transition between the photoproduction and the standard DlSgp regime at HERA. 
The three ranges of the squared mass of the virtual photon were studied: < 10~^ 
GeV^ (1994 data), 0.65 < < 20 GeV^ (the 1995 shifted vertex data) and 9 < < 
49 GeV^ (1994 data). The jets with > A GeV (for the photoproduction > b 
GeV only) for —2.5 < rj* < —0.5 were measured. The kx - clustering algorithm was 
used. For the photoproduction PHOJET 1.03, for DISep events - LEPTO 6.5 and 
ARIADNE 4.08 were used in simulation. RAPGAP and HERWIG 5.9 (with p??^" = 
1.5 GeV) were applied to model direct and resolved real and virtual photon processes. 
The GRV94 HO parton parametrization was taken for the proton. GRV HO with the 
Drees - Godbole parametrization (cj^ = 1 GeV^), denoted GRV HO/Drees-Godbole, 
and SaS2D parametrizations for the virtual photon were used. 

The measured da^p/dE^ as a function of the transverse energy of the jet for various 
P^ ranges integrated over y and rj* ranges is presented in fig. p.47| . It was found to 
be in agreement with the HERWIG model (using the GRV HO/Drees-Godbole parton 
parametrization). In fig. |148| the corresponding data for the rapidity distribution for 



jets with E^ > 5 GeV are shown in the various virtuality ranges. 
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Figure 147: The inclusive dcTep/dE^ jet cross section as a function of the transverse energy 
Ej^ for various initial photon virtuality (denoted as Q'^) ranges and for —2.5 < rj* < —0.5. 
The HERWIG GRV HO/Drees-Godbole prediction is denoted by the solid line, the dashed 
line corresponds to the direct contribution to this prediction (from 
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Figure 148: The inclusive daep/drj* jet cross section as a function of the rapidity rj* for 
various initial photon virtuality (denoted as j ranges for the > 5 GeV. The HER- 
WIG (with the GRV HO/Drees-Godhole parton parametrization) prediction is denoted by 
the solid line, the dashed line denotes the direct contribution (from j^/j. 



To study the dependence of the virtuahty of the photon, the flux (|3^ is introduced 



and the cross section calculated 



0". 



(35) 



Note that it is not certain that the implied factorization really holds for the whole 
range of kinematical variables. 

The results for a^*p agree with the HERWIG and RAPGAP predictions based on 
the GRV HO/Drees-Godbole and SaS2D parton parametrizations in the virtual pho- 
ton, as shown in fig. 149 . Similar comparison with LEPTO and ARIADNE was made 
(not shown) and "neither model can describe the data when < and the virtual 
photon can be resolved" . 
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Figure 149: The inclusive 7*p jet cross section as a function of (denoted as Q"^) for 
various ranges of the transverse jet energy E^, for —2.5 < rj* < —0.5. The HERWIG 
(with the GRV HO/Drees-Godbole distribution) prediction is denoted by the dashed hne, 
RAPGAP (with the GRV HO/Drees-Godbole distribution) - the dotted line, RAPGAP (with 
the SaS2D distribution) - the solid line, and the dot-dashed line corresponds to HERWIG 
with the GRV HO parametrization as for the real photon ("no P^ suppression") (from 



Comment: "The inclusive jet cross-section can therefore be understood if a partonic 
structure is ascribed to the virtual photon. 

The data are best described by the RAPGAP model using the SaS2D parametrization 
of the virtual photon. " 

•HI 2000b (HERA) 

The dijet event rates have been measured for 5 ^ ^ 100 GeV^ and 10""* ^ xbj ^ 
10~^ [i.e. for the DISep events at small xbj), and for jet transverse momenta squared 
Pt ~ (see also HI 94). The data collected in 1994 correspond to y > 0.05. The 
condition | Arj* \< 2 was imposed. The jet pt was assumed to be at least 5 GeV, with 
the following requirements: for the symmetric > 5 GeV, for the asymmetric 

case p^^ > 5 GeV and p^2 — GeV, and for the sum scenario p^^ + p^2 — 13 GeV. 
The cone algorithm with R=l was applied in the r]* - (p* plane. 

The standard DISep Monte Carlo generators: LEPTO and ARIADNE were used 
to describe in LO the direct 7*g contribution with two sets of parton parametrization 
for the proton: MRS-H and GRV94 HO. The additional mechanism based on the 
resolved virtual photon interaction neglecting the longitudinal 7* was introduced in 
the analysis. Its contribution was obtained using the RAPGAP Monte Carlo model 
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with the factorization and renormahzation scale = + p\. To describe this 
contribution also the NLO calculations (for partons) in JeTViP program with SaSlD 
parton parametrization for 7* were studied. For comparison the DISENT program was 
used also to calculate the NLO direct virtual photon contribution. The factorization 
and renormahzation scale = + 50 GeV^ was used in both NLO calculations. 

The transverse energy flow with respect to the jet axis was studied. There is a good 
agreement between the data and the events simulated with LEPTO and ARIADNE, 
except for the rj* distribution of the jets (not shown). 

The comparison of the fraction of dijet events in all DISep events with the RAPGAP 



and ARIADNE models prediction can be found in fig. |150| , where the rates as a 
function of (integrated over xsj) and xbj (integrated over P^) are shown. The 
same distributions are shown in fig. |151| for the NLO calculations obtained using 
DISENT (the direct photon contribution) and JeTViP, where also the resolved photon 
component is included. The x°^'^ distribution is shown in fig. |152. 



Comment: The agreement with the data at low xbj and moderate P^ is improved when 
beside the direct also contributions from the resolved virtual photon are included in the 
Monte Carlo model (RAPGAP) and in the NLO calculations. "The CDM model, as 
implemented in ARIADNE, is also able to describe the dijet rate well. " 
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Figure 150: The dijet rate as a function of (denoted as Q^) (a,c,e), integrated over xsj, 
and as a function of xbj (b,d,f), integrated over P^. The symmetric (a,b), asymmetric (c,d) 
and sum (e,f) cut scenarios on the of the two jets are presented. The comparison is 
made with RAPGAP and ARIADNE predictions for the direct and resolved virtual photon 
contributions (from J^^/ j. 
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Figure 151: The dijet rate as a function of (denoted as Q"^) (a,cj, integrated over xbj, 
and as a function of xbj (b,d), integrated over P^. The data (the same as in fig. llSQ c-f) 
for the asymmetric (a,b) and the sum (c,d) cut scenario are compared to different NLO 
calculations for the direct and resolved virtual photon contributions (from j^). 
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Figure 152: The uncorrected distribution of x° in three different P (denoted as Q ) bins 
(a). In b), c) and d) the data in different P^ bins are compared to the RAPGAP and 
ARIADNE models for the direct and resolved virtual photon contributions (from jp^. 
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•HI 2000c |]100| (HERA) 

The dijet cross sections were measured in 1996 for 1.6 < < 80 GeV^, 0.1 < y < 0.7, 
and 30 < < 300 GeV^. For jets the ranges \At]\ < 1, -3.0 < f] < -0.5, and 
2 \AE\/E < 0.25 were assumed. Two Monte Carlo generators were used, HERWIG 
and RAPGAP, with different (LO and NLO) virtual photon parton densities: GRV 
HO/Drees-Godbole (p^" = 3 GeV), GRV LO/Drees-Godbole (p^'^ = 2 GeV), and 
SaSl,2D (p?^^" = 2 GeV). For proton, the GRV HO and GRV LO parton parametriza- 
tions were used. The fcr - clustering algorithm for jets was applied. Only the trans- 
versely polarized photon flux was assumed in the analysis. 

The data were used to extract (for the first time) an effective LO parton density 
for the virtual photon. 

The transverse energy flow was studied as a function of 6(j) and compared with 
the Monte Carlo predictions with and without soft underlying event (HERWIG and 



RAPGAP, respectively), see fig. f53. 
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Figure 153: The transverse energy How as a function of S(j) for different rj bins. Comparison 
with predictions of HERWIG (HO) using the GRV HO/Drees-Godhole parton distribution 
(dashed hne) and RAPGAP (HO) using the SaS2D parton distribution (sohd line) is shown 
(from jjlOdi] ). 



The triple differential cross sections as a function of x^, 
presented in figs 



E^ and P2( 
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Q^) are 

155| and |156| , respectively, together with predictions of the 



HERWIG and RAPGAP Monte Carlo programs. 
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Figure 154: The differential cross section as a function of in four bins of the virtuahty 
(denoted as Q^) and four bins of E^. The comparison with the LO prediction of the 
HERWIG Monte Carlo based on the GRV LO/Drees-Godbole parametrization with uj=0.1 
(solid line) and lu=0.2 (dashed line) is shown. Also the direct contribution is displayed 
(dotted line) (from ^f). 
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Table 28: 



p2 

[GeV^] 






( <ifnf —1— fill 1 
\^vIajv. \^ oyov. 1 


2.4 
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E /GeV' 

Figure 155: The differential cross section as a function of Ej> in four bins of the virtuality 
(denoted as Q^j and five bins of x-y. The comparison with the LO prediction of the 
RAPGAP Monte Carlo using the SaSlD (dashed line), SaS2D (dash-dotted line) and GRV 
LO/Drees-Godbole (solid line) parton distributions in the photon is shown (from / |iOC|/ }. 
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Figure 156: The differential cross section as a function of virtuality (denoted as Q^) 
in four bins of Ej^ and five bins of tlie Xj. The comparison with the LO prediction of 
the RAPGAP Monte Carlo: using the GRV LO/Drees-Godbole parton distributions in the 
photon (sohd line) and the the direct contribution only (dashed line) is shown (from fiOOjl). 



The Single Effective Subprocess Approximation |5J] was assumed in the analysis 
and the effective (LO) parton density of the virtual photon 



(36) 



with q'^ = + (p)j was extracted for the first time. Here q^{(p) and are photon 
polarization-averaged parton densities. This effective density is presented in table |2^. 
In figs. p.57| , p.58| and p.59| it is displayed as a function of x^, pi (playing the role of 
the hard scale Q^) and P^, respectively. Only points where averaged pi is larger than 
< > are shown. 



In fig. |160| the obtained LO effective parton density as a function of the virtuality 



is shown for x^=0.425 and 0.6, at 
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Figure 157: The effective parton density as a function of x-y for various values of the 
parton transverse momentum squared (= pi) and virtuality (denoted as Q^). The 
comparison with predictions of the SaS2D (dotted hne), SaSlD (dashed hne) and Drees- 
Godbole (with GRV LO parametrization for the real photon and using uj = 0.1, solid line) 
parton distributions in the photon is shown (from / [J06|/ j. 
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Figure 158: The effective parton density as a function of (= pj^ ) for various values of x-y 
and virtuality (denoted as Q^). Comparison with predictions as in fig. |J57| (from lilOOjj ). 
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Figure 159: The effective parton density as a function of (denoted as Q'^) for various 
(= pI) and x^. Comparison with predictions as in fig. |J5?f (from |JOC|/ j. 
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Figure 160: The LO effective parton density as a function of (denoted as Q^j for (= 
p1) =85 GeV^ , x^=0.425 and 0.60. Comparison with predictions as in fig. |J57| (from fiOOi j). 

Comment: "Both simulations give a good description of the energy flow in the core of 
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the jets", neither Monte Carlo model is able to describe the pedestal for all ranges of 
Et and r] (fig. \T5^ ). 



•ZEUS 93 UTOIP (HERA) 

The evidence for events with two or more jets in DISep (neutral current interaction) is 
presented for > A GeV^ (1992 data). 



ZEUS 95d PU^ (HERA) 



The two-jet production in DISep for 160 < < 1280 GeV^ 0.01 < xbj < 0.1 and 
0.04 < y < 0.95 was studied. In reconstructing jets the JADE algorithm was used (in 
the HERA laboratory frame). The LEPTO 6.1 generator was used in the analysis. 
The comparison was made with the NLO calculations. 

• ZEUS 2000b gi03|] (HERA) 

Dijet production for photon virtualities from up to 4.5 GeV^ was measured, based 
on the data collected in 1995. Events corresponding to 0.2 < y < 0.55, -Eti,2 > 5.5 
GeV and —1.125 < rj < 2.2 were studied in three bins: ^ (a quasi-real 
photon sample obtained by antitagging condition, which corresponds to P^ < 1.0 
GeV^ with median 10^^ GeV^), 0.1 - 0.55 GeV^ and 1.5 - 4.5 GeV^. The kr algorithm 
and the HERWIG 5.9 program (with p??^"=2.5 GeV), with the GRV LO (photon) 
and MRSA (proton) parton parametrizations were used in the analysis of data. The 
possible multiple parton interaction ('underlying event') was simulated for the quasi- 
real photon sample (denoted MI). To estimate uncertainties due to modelling of the 
jet fragmentation the PYTHIA model was used in addition. The LEPTO model for 
the direct 7*g events and the NLO prediction for the resolved 7* processes (JeTViP) 
are compared with data. 

Uncorrected distributions are presented in fig. |161| for real and for virtual pho- 
tons (in two bins), and compared to the HERWIG predictions for the direct and 
resolved (LO) photon contributions. 

In the low region the data for quasi-real photons disagree with MC prediction. 
"Disagreement is also observed in the rj, ysj, and Et distributions (not shown) and 
can be attributed to the presence of underlying event effects or uncertainty in the 
PDFs of the photon." The inclusion of MI is not enough to reproduce the shape of the 
low data (fig. |161| a), so a reweighting has to be performed to include this fact in the 
correction of the data for migrations and acceptance. The obtained results (showing 
an agreement with data) are presented in fig. |161| a (dashed histogram), also the rj, 



yBj, and Et distributions (not shown) "agree well with data" 
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Figure 161: Uncorrected x'^'^ distributions for various photon virtualities: a) (denoted 
here as Q'^)=0, b) 0.1 < < 0.55 GeV^ , and c) 1.5 < P^ < 4.5 GeV^ . Events with Eti,2 > 
5.5 GeV, —1.125 < < 2.2 were included. The predictions of the HERWIG program (GRV 
LO) are shown for: a direct contribution (shaded histogram), a direct + resolved + MI (only 
for P^ = sample) contribution (solid histogram), in a) a reweighted prediction for a direct 
+ resolved + MI contribution (dashed histogram) (from |JOj|/ j. 



The dijet cross section da/dx^, corrected to the hadron level, is presented in fig. |162 
for a low Et sample, where, as above, Eti,2 > 5.5 GeV and —1.125 < rj < 2.2. The 
data are compared with prediction of the HERWIG program, with different parton 
parametrizations for the photon: GRV LO and WHIT2 (without MI and with MI, then 
assuming = 2 GeV), which are valid for a real photon, and SaSlD. In addition a 
comparison with prediction of the LEPTO program is shown for events with > 1.5 
GeV^. 
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Figure 162: The dcj/dx'^'^ cross section for various photon virtuahties: a) (denoted here 
as Q'^)=0, b) 0.1 < P2 < 0.55 GeV^, and c) 1.5 < < 4.5 GeV^, for low Et sample: 
Eti,2 > 5.5 GeV, and -1.125 < rj < 2.2. The predictions of the HERWIG program with 
GRV LO (SaSlD) are shown by dotted (solid) histograms. In a) the HERWIG predictions 
obtained using the WHIT2 parton paramerization with and without MI (dot-dashed and 
dashed histograms) are shown, in c) the LEPTO prediction (dot-dashed) is given. The 
shaded bound represents uncertainties in the simulation (from ^103(j). 



Similar analysis has been performed for the high Et sample, where Ej-i > 7.5 GeV, 
Et2 > 6.5 GeV and -1.125 < r] < 1.875, see fig. |16| for results. 
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Figure 163: The da/dx"''^ cross section for various photon virtuahties: a) (denoted here 
as Q'^)=0, b) 0.1 < < 0.55 GeV^, and c) 1.5 < P"^ < 4.5 GeV^ , for high Et sample: 
Eti > 7.5 GeV, Et2 > 6.5 GeV, and -1.125 < r? < 1.875. The predictions of the HERWIG 
program with GRV LO (SaSlD) are shown by dotted (solid) histograms. In c) the LEPTO 
prediction (dot-dashed histogram) is given. The shaded bound represents uncertainties of 
the simulation (from |JOj|/ ). 



The ratio a{x^ < 0.75)/ a{x^ > 0.75) as a function of the photon virtuahty is 
shown in fig. |164| for both sets of Et cuts. Here a comparison with the prediction of the 
HERWIG program with the GRV LO and SaSlD parton parametrizations is shown. 
The small region is not properly modelled. Also the LEPTO program predictions 
are shown for P^ > 1.5 GeV. The NLO QCD predictions obtained by using the JeTViP 
program with the SaSlD and GS96/Drees-Godbole parton parametrizations and the 
squared hard scale equal to P^ + p"^, are compared to high Et data. They "lie well 
below the data." 
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Figure 164: The ratio a{x'^^ < 0.75)/cj(x"^'* > 0.75) as a function of the photon virtuality 
(denoted as Q'^) for both Et cuts. The predictions of the HERWIG program with the 
GRV LO(SaSlD) parton parametrization are presented as dotted (sohd) histograms. The 
LEPTO program predictions are shown for P^ > 1.5 GeV (dot-dashed histograms). In b) 
the NLO QCD prediction of the JeTViP program with the SaSlD (GS96/Drees-Godbole) 
parton parametrization with the squared hard scale P^ +Pt shown as thick sohd (dashed) 



histogram (from ilOcj). 



Comment: "The shape of the dijet cross sections, da/dx^, is compared to the prediction 
of HERWIG MC for a variety of photon PDFs. None of these models is able to explain 
the data for both high- and Iow-Et cuts in all . " 

For the dependence of the ratio cr{x^ < 0.75)/ cr(x^ > 0.75) on the photon virtuality 
P^..." none of the LO model, or the NLO calculation examined here, gives a good 
description of the data across the full kinematic region. " 

Note, that data were corrected for acceptance, smearing and kinematical cut using 
HERWIG 5.9 with the GRV LO parton parametrization, valid for a real 7. 
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3.4 Forward jet (particle) production in collisions and the 
resolved virtual photon contribution 

The production of high Ef jets (particles) in the direction of the proton in events with 
large virtualities, i.e. DISep events, small xbj and ~ was studied by both the 
HI and ZEUS groups at the HERA collider. 

The contribution due the standard DGLAP evolution based on the strong pt or- 
dering of parton emission along the whole chain between the proton and the photon 
is expected to be strongly suppressed for such events. In contrast, signals from the 
non-strongly ordered in transverse momenta partons from the proton (the BFKL 



or CCFM evolution) are expected to show up pOO|| . At first the BFKL calcula- 
tion for the forward jet (particle) production was based on the LO approach, recently 
the modified LO BFKL approach has appeared including the kinematical constraint, 
called the consistencv constraint , on the gluon emission. This constraint presumably 
embodies a major part of the NLO 1/x corrections |p03[, see also discussion on this 



point in ||104|| . Besides the BFKL, the mechanism with the two DGLAP evolution 
chains corresponding to the resolved virtual photon-proton interaction may also lead 
to forward jets (or particles), as was pointed out in |p01|| . The NLO prediction of 



this contribution, based on the parton level generator JeTViP ||196|| , was performed in 



199|| . Several Monte Carlo programs were used to model the radiation characteristics 



for a particular evolution. The HERWIG and LEPTO models are based on DGLAP, 
while LDC (called also LDCMC) - on the CCFM equation. In the ARIADNE pro- 
gram, which implements the Color Dipole Model, there is no pt ordering, so it has the 
properties of the BFKL evolution. The RAPGAP generator (also JeTViP) contains 
the resolved virtual photon contribution with the DGLAP cascades. 

Both forward jet and forward single particle production are measured. The ad- 
vantage of studying single particles is the potential to reach a smaller angle than it is 
possible with jets. 

The scale for the deep inelastic scattering on the proton (DlSgp), Q^, plays the 
role of P^, if the contribution due to the resolved virtual photon is considered. For 
the sake of transparency, in this section we keep notation for the virtuality of the 
photon. 
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DATA 



HI 95b |1T05[| (HERA) 



The first analysis of the forward jets in DISep was done. 
•HI 97c (HERA) 

The pt distribution of charged particles in DlSgp events was studied. The hadronic 
final state was analysed in order to test the different schemes for parton evolution in 
the proton (DGLAP and BFKL). 

Comment: The large parton radiation between the current and the proton remnant 
system was found. 



•HI 99b |]T07[| (HERA) 

The production of forward jets and forward particles (the charged particles and the 
neutral pions) with large pt is analysed in the DlSgp events at small xbj (1994 data). 
The forward particles were measured for 5*^ < ^ < 25° (LAB), and for the energy ratio 
Xcp{n) = Ecp(TT)/Ep between 0.01 and 0.15 (for charged particles below 0.015 only). 
Two Pt thresholds were used, 1 and 2 GeV (only for tt^), with additional requirement 
0.0002 < XBj < 0.00235. Forward jets with the energy ratio Xjet = Ej^t/Ep > 0.035, 
and 0.5 < Px/Q'^ < 2, 7° < 6jet < 20° and pt > 3.5 and 5 GeV were studied 
(10~^ < XBj < 4 ■ 10""^). A simple cone algorithm with R=l was used for jets. 

The predictions of the LO BFKL calculation ||202|| and of the Monte Carlo gen- 
erators: LEPTO 6.5 (with renormalization scale Q^), LDC 1.0 and ARIADNE 4.08 
(scale Pt) were compared with data. The RAPGAP 2.06 generator with resolved 
photon processes (with SaS2D parton parametrization for the photon, assuming scale 

+ Pt) ^l^o used in comparison. To describe a proton the GRV HO parton 
densities were used in all the programs. The transverse energy flow as a function of 



6ri was measured for jets and results are presented in fig. 165 
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Figure 165: The transverse energy How as a function of Scf) and Srj for two pT cuts for jets 
(from jjlO^] ). 



The comparison of the measured XBj distribution with prediction of different types 



169 



of Monte Carlo models for single particle spectra is shown in fig. p.66| . Here also the 
comparison with the results of the analytical calculation based on the BFKL approach 
202|| , with the factor | to describe the data, is made. 




Figure 166: The single particle spectra as a function of XBj- aj for charged particles, b) for 
charged pions and n^'s; c), d) for tt^'s with different cuts on x-,^ and pj^. Comparison with 
the Monte Carlo predictions, and in panels (c, d) in addition with the LO BFKL calculation 
J20^ , is presented (from JJO?{ /j. 



In fig. |167| the forward jet data for pt > 3.5 GeV and pt > 5 GeV are shown and 
compared to the Monte Carlo predictions and analytical calculations, LO BFKL and 
the NLO direct contribution (implemented in the DISENT program). 
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Figure 167: The forward jet cross section as a function of xsj for px cuts: 3.5 GeV (a, 
b) and 5 GeV (c, d). a) and c) Comparison with the Monte Carlo predictions: RAPGAP, 
LDCMC, LEPTO and ARIADNE, b) and d) Comparison with the LO BFKL calculation 
20^] and NLO direct prediction (DISENT) (from ^]). 



Fig. |168| (taken from ||199|| ) shows the data together with the standard DISep NLO 



predictions, and the NLO calculation by Kramer and Potter [ 199|| , where the interac- 
tion due to partonic content of the virtual photon is included. Note that in this calcu- 
lation the hard scale Q"^ is assumed in the form Q"^ = kM"^ + P^, with =< Et >^ 
fixed to be 50 GeV^ and parameter k between | and 3. 
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Figure 168: Results of the NLO calculations of the forward jet cross section as a function 
of XBj, compared to the HI data for Et > 3.5 GeV (a, b) and Et > 5 GeV (c, d). The 
NLO description within the DISep approach is presented in a) and c), while b) and d) 
show the NLO description based on the direct (DIR) and the resolved (RES) virtual photon 
contributions. Three different renormalization scales are used: Q^(= /s^) = M'^ + (solid 
line), Q2(= ^ 3jvf2 _^ q2 ('c/agjjej j^^ej and Q'^{= fi"^) = + (dotted line), with 

M2 =< Et >^=50 GeV^ (from Jw^). 

The azimuthal correlation A(f) between the jet and the scattered electron was also 
studied (not shown). 



Comment: "The conclusions obtained using the forward jet and the single particle 
measurements are in agreement" . 

Predictions based on the DGLAP approach fail to describe the data, except for those 
which allow for resolved virtual photon contributions. The RAPGAP model gives a 
good description of data. The BFKL calculation is above the measurements for the 
data. ARIADNE model reproduces the rise towards low x. 

'Predictions of a CCFM evolution, which should smoothly interpolate between the 
DGLAP and BFKL regimes, give a poor description of all measured distributions. " 

•HI 99c JlOg (HERA) 

The forward production of large pt 7r°-mesons was measured in 1996 in DlSep events 
at small xbj- The events correspond to 2.0 < < 70 GeV^, 0.1 < y < 0.6 and 
small polar angle between 5° and 25° (LAB) . Events with forward 7r° production were 
collected for p*j? > 2.5 GeV (in the hadron CMS) and = E^/Ep > 0.01. 

The comparison of the data is made with the prediction of the modified LO BFKL 
calculation [E03I1 and with the Monte Carlo models: LEPTO 6.5 based on the DGLAP 



description of the direct 7*g scattering, and RAPGAP 2.06, which includes the resolved 
photon contributions and uses the DGLAP evolution for 7* (SaSlD with = Q'^+Pt)- 
For the proton all models use the CTEQ4M parton parametrization. 
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In fig. |169| a the xbj dependence is presented for > 2.5 GeV as a function of 
XBj in bins of Q^. The observed strong rise of the cross section towards small xbj in 
the forward 7r° production is identical to the rise of the inclusive ep cross section for 
all events, see fig. |169| b. Comparison with the modified LO BFKL calculation ||203|| , 
and with LEPTO and RAPGAP Monte Carlo models are also present in the figure. 
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Figure 169: The inclusive forward cross section for > 2.5 GeV as a function of 
XBj in bins of . Comparison with the modified LO BFKL calculation (solid line) |203|/ , 
with the LEPTO (dotted line) and RAPGAP (dashed line) Monte Carlo models, a) The 
XBj distributions for ir^. b) The ratio of the XBj distributions presented in a) and the 
corresponding ones for the full sample of ep events (from jlOSi] ). 



The analysis of the p^ and rj distributions in the bins was performed and the 



results are presented in fig. 170 



173 



Q. 



-D 200 



T3 



100 



200 



100 



200 



100 



/ GeV^ 



[2.0... 4.5] • H1 data 

p; ,>2.5 GeV 




[ 15.0 ... 70.0] 



a) 



■■■ LEPT0 6.5 
■ - RAPGAP 2.06 

— mod. LO BFKL 



J I I \ I I I \ I I I I I I L 



1.5 



2.5 



3 



> 

B 10' 



% 10 



• H1 data 

p;,,>2.5GeV 



T3 



10' 



10 



10' 



10 



/ GeV^ 



[2.0... 4.5] 

I I I I L 



[4.5... 15.0] 

I I I I L 



b) 



[ 15.0 ... 70.0] 

I I L 



3 4 5 6 7 8 



Figure 170: The inclusive forward vr^ cross section for > 2.5 GeV as a function of 7]^^ (a) 
and of pj^ (b), in bins of Q^. Comparison with model predictions as in fig. \16d{ (from fiOS^j). 
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In fig. |171| the dependence is presented for events witli > 2.5 GeV, in 
comparison witli tlie modified LO BFKL and Monte Carlo predictions. 
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Figure 171: The inclusive forward tt^ cross section for > 2.5 GeV as a function of 



Comparison with model predictions as in fig. |J6£| (from llOajj ) 
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For comparison also the Xbj and dependences are presented for the threshold 



equal to 3.5 GeV, see fig. |T72 
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Figure 172: The inclusive forward tt^ cross section for > 3.5 GeV as a function of (a) 
XBj cLnd (b) Q^. Comparison with model predictions as in fig. |i6S| (from l[103(l ). 



Comment: The measurement extends down to xbj = 5 ■ 10~^, lower than the previous 
values. It covers lower polar angles and larger p^ range, as compared with HI 99b. 

"Inclusion of processes in which the virtual photon is resolved improves the agree- 
ment with the data, but does not provide a satisfactory description in the full xbj and 

range. A calculation based on the BFKL formalism is in good agreement with the 
data, particularly for the shape description, but the absolute normalization remains 
strongly affected by the scale uncertainty. " 



•ZEUS 99c |]T09f| (HERA) 

In the e+p scattering the forward jets produced with y > 0.1, > 5 GeV and xbj 
between 4.5 ■ 10~^ and 4.5 ■ 10~^ are studied for the 0.5 < E^/Q'^ < 2. Data were taken 
in 1995. The motivation was the search for effects of the BFKL type parton shower 
evolution. The jet was defined according to the cone algorithm (R=l). The LO BFKL 
and few Monte Carlo programs (ARIADNE 4.08, LEPTO 6.5, HERWIG 5.9, LDC 1.0 
and MEPJET for the NLO direct contribution) predictions were compared with data. 

The parton-level predictions for forward jet cross sections as a function of xbj are 
shown in fig. |173| (left) and for the hadron level in fig. |173| (right). The comparison 



with the LO BFKL and LO BFKL 1** term [i.e. with one parton rung in the gluon 
ladder) and the Monte Carlo predictions is also given. 

The E'^/Q'^ distribution was studied as well (see ZEUS 2000c for the new results). 
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Parton level 




Figure 173: Left: The xbj distribution at the parton level obtained in various approaches 
(Monte Carlo and QCD analytical calculations). The shaded band corresponds to the range 
of results obtained with MEPJET using a renormalisation scale between 0.25A;j. and 2.0/cj,. 
Right: The xbj distribution at the hadron level: data points and various Monte Carlo 
predictions (from / [iO£|/ }. 



Comment: "Three regions are identified in the E^/Q"^ distribution: i) the standard 
DGLAP region with <^ Q^, where all Monte Carlo models are in agreement with 
the data; ii) the region of phase space where BFKL dynamics is expected to contribute 
significantly with E^ ^ Q^, where only the Colour Dipole model describes the data 
well; Hi) the region with E^ ^ Q"^ , where none of the models reproduces the data. " 
No contribution due to partonic content of the virtual photon was considered. 

•ZEUS 2000c iTToj (HERA) 

The forward jet production (1995 data) was studied based on in DISep events with 
> 10 GeV2, y >0.1. Events with Et >5 GeV, 0.00025 < XBj < 8.0 ■ 10^^ ^nd 
E"^ IQ"^ between 10~^ and 10^ were considered, in addition the cut Xjet=Pz,jet/Pproton > 
0.036 was used to select forward jets. The cone algorithm was applied to jets in LAB 
frame with R=l. 

The LO MC models, standard for the description of the DlSgp events: LEPTO 
6.5, HERWIG 5.9, LDC 1.0 (the linked-dipole chain as in CCFM), and ARIADNE 
4.08 were used. In addition the RAPGAP 2.06 Monte Carlo program and the JeTViP 
(NLO calculation) program were included in the analysis to describe the resolved vir- 
tual photon interaction, with SaS2D and SaSlD parton parametrization for the photon, 
respectively. Both introduce the renormalization/factorization scale = E^ + Q^. 
In fig. p.74| the E^/Q"^ distribution, and in fig. p.75| the x^ distribution is shown and 
compared with the Monte Carlo models, with and without the resolved photon con- 
tribution. 
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Figure 174: Forward jet production cross section as a function of Ej./Q'^. Comparison with 
predictions of the various MC models, left: ARIADNE 4.08, LEPTO 6.5, HERWIG 5.9, 
LDC 1.0, and right: RAPGAP 2.06, JeTViP (from ^TWj] ). 
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Figure 175: Forward jet production cross section as a function ofxsj: data compared to the 
various Monte Carlo predictions: LEPTO based on the standard DGLAP evolution (dashed 
histogram), RAPGAP with the direct virtual photon contribution (dotted histogram) and 
RAPGAP with the direct and the resolved photon contribution (solid histogram) (from 
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In fig. |176| (taken from | |199| |) the data are compared to the NLO calculation | |199| , 
with and without resolved virtual photon interaction (and with different renormaliza- 
tion scales). 




Figure 176: Results of the NLO calculations of the forward jet cross section for Et > 5 
GeV, without (left) and with (right) resolved virtual photon interaction, compared to the 
ZEUS data. Three different renormalization scales are used: /x^ = + (solid line), /j? 
= 3M2 + Q2 ('cJagjjej line) and //^ = + (dotted line), with =< Et >'^=50 

GeV^ (from p^j. 



Comment: "All models describe DIS regime. In the BFKL regime (Q^ ~ E\) ARI- 
ADNE reproduces the measured distributions well. Only model with resolved photon 
component, RAPGAP, describes data everywhere. " 
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4 Leptonic structure functions of the photon 



The lepton pair production by two photons in e"'"e~ scattering, e~^e~ — > e~^e~l~^l~, is 
one of the basic test of QED0. It serves also as a test of experimental methods applied 
to more complicated two-photon reactions involving hadrons, especially to the process 
e~^e~ —>■ hadrons, where the hadronic structure functions of the photon are measured. 
To study the reliability of unfolding and tagging applied in these measurement sQ, 
including the effect of the (small) virtuality of the quasi-real target photon, the 
lepton pair productions can be used as an ideal tool. 

Single-tag events with lepton pair in the final state can be described in the unpo- 
larized e^e~ collision in terms of leptonic structure functions , p^^'^^^^ [=p^^^^^'>'^_^ 
p-yiQED) p-yjQED) ^ analogy to the hadronic case. The leptonic structure functions 
were measured already in early eighties^. Although the production of final states with 
electrons, muons and even taons have been measured, only the muonic case allows for 
the extraction of the structure function from the data. Also nowadays experiments 
at LEP measure only muonic structure functions. The invariant mass of the 
pair, and hence x^j, can be determined very accurately and the measurements of the 
muonic structure function F2 , denoted below asF2^(«^^),are only statistically limited 
in contrast to the hadronic F2 . Recently also two other (muonic) structure functions, 
FJ^'?^^) and have been measured at LEP. 

Recently, for the first time, a measurement of the leptonic final state mediated 
by two highly virtual photons (double-tag events) was performed at LEP. For both 
photons with high virtualities the presence of large interference terms was established 
(OPAL 99b) even for events with very different virtualities ^ P^. Thus in this 
case the (effective) leptonic structure function for a virtual photon cannot be intro- 
duced, nevertheless we report these results together with leptonic structure functions 
of the real photon. 



4.1 Theoretical description 

The QED structure functions for photon: F^^^'^^\ piiQm ^ pi(QED) ^ 

etc., unlike in 

the hadronic case, can be reliably calculated in QED. Moreover all the particles in the 
final state can, in principle, be directly observed. 

The p2^'^^^^ structure function is much harder to measure than p^^^^^^ ^ because 
its contribution to the relevant cross section is weighted by the small factor y"^ (as 
for hadronic F2, see eq. (^). However, this longitudinal structure function is not 
the only structure function that contains additional information about the structure 
of the photon. It has been shown that there are azimuthal correlations in the final 
state particles from two-photon collisions which are sensitive to additional structure 
functions |[113| , |114|| . Below we discuss some of them for the real photon. 



Beside the cross section '^^^^.J^ (eqs. (|^), (^)), one can also measure in the two- 
body 7*7 process one final state particle a. Then additional structure functions Fj 
(resulting from interference between the amplitudes involving both longitudinal and 



^^See ref. |111| for the review of the early experiments. 

^•^See discussion in sees. 2.2 and 2.3 

34 



The earhest measurements are summed up in (see also 
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transverse photons) and (describing interference between amplitudes for transverse 
photons only) appear ||114|| : 

da[e, ^ eaX) 2.a^ 1 + (1 - y)^^^^^^^:^ ^ ^^^^^.^ 



dxBjdydVla/A'K xbjV 



-p{v)Fa cos 0a + ^(^iy)F^ cos 20„ 



(37) 



Here ^la describes the solid angle for the particle a in the 7*7 centre of mass frame, 
and (pa is its azimuthal angle around the 7*7 axis, relative to the electron (tag) plane. 
The corresponding polar angle relative to the 7*7 collision axis is denoted by 6. The 
functions e{y) = 2(1 - y)/{l + (1 - y)^) and p{y) = (2 - y)./T^/{l + (1 - y)^) are 
close to 1, as they are ~ (1 — 0{y'^)). The standard functions Fj, and F2 (and also 
F\^) are obtained from the corresponding Fj^ {Pab) integration over the solid 
angle Via- Note that the formula (|37|) holds for two leptons (or two partons) produced 
in the final state. 

The expression for F^^^^^^ and F^^^^^^ (or F^^^^^^^ = F^^^^^^ + 2x^,F^(<3^^)) 
in the lowest order QED is given by the Bethe-Heitler type formula (see eq. P). The 
functions Fj^'^^^^ and F^^^^^^ calculated in the same approximation and with the full 
dependence on the muon mass up to terms of order m'j^/W'^ |p05|| , have the following 
form: 

F;(<3^^)(x^„/?) = -XBj^XB,{l-XB,Kl - 2xBjmi + (1 - Pf^^^] (38) 
TT * i — ZXBj 

+?^^V^V^l^arccos(vT^)}, 

i — ZXBj 



^r''^\xB,.P) = -X|,(l - XB,m^ - (1 - 

vr Zx Bj 



2 X Bj ZX Bj 1 — P 



where 



P 



The collection of other useful formulae for the above structure functions in the 
lowest order QED approximation with the full dependence on the fermion mass can 
be found in p05|], see also ||2U6|| for the first calculation in the lowest order QED. 



The function f'b'^^^^ is, in a specific logarithmic approximation [|115|| and zero 



muon mass limit, accidentally equal to F]}'°^^^\ although it involves quite different 
photon helicity structures. Thus extracting F^^'^^^^ can give us indirectly information 
on F]}^^^^ (in this approximation). 

The above description concerns only the leptonic production in 7*7 collisions, so 
with one (almost) real photon. Similar consideration can be performed for 7*7* col- 
lisions, when the notion of leptonic structure functions of the virtual photon can in 
principle be applied if Q"^ ^ P^, in full analogy to hadronic structure functions for 
the virtual photon (sec. The extraction of such functions from the corresponding 
cross section, e.g. F^ (Q^^)^ jg feasible if the interference terms can be neglected. 
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4.2 Measurements of single and double-tag events 

A measurement of single-tag events allows to extract p^^'^^^^ ^ as well as pj^^^^^ and 



p^iQED) ^ The function pj^^^^^ is antisymmetric in cos^* j^. To extract it from the 



h-yiQED) 
'A 



data, the x^j-dependence of the term in the differential cross section proportional to 
pjiQED) j^^g studied separately for the ranges cos 6* < and cos 6** > in order 
to avoid cancellation. 

In the last few years measurements of the leptonic structure function p^^^^^^ have 
been performed at LEP 1 by the four collaborations; the ratios of structure functions 



A,B 



to f; 



liQED) 



have also been measured by all the groups. Only the L3 and 



OPAL groups provided also the results for individual structure functions F4 



liQED) 



and 



The Vermaseren, Galuga and BDK Monte Carlo generators, which have imple- 
mented full dependence on both photons' virtuahties, were used in analyses (for back- 
ground simulation also other generators were used). 

Results of all p^^'^^^^ measurements discussed here will be shown together in 



figs. |189| and |190| . A similar compilation is also presented for the ratios Fj^ / Fr 



■yiQED) I r^-i{QED) 



(fig. [191]) and for individual functions pf^'^°\ p^iQ^^) (fig. [T92D . 

The relevant cross section for the double-tag events, i. e. for the production of the 
leptonic (fi^fi^) final state in 7*7* collision, is as for the hadronic final state (eq. pHD , 
with the same definition of kinematical variables. Note that for these events the ad- 
ditional angle between two scattering planes for the two scattered electrons appears. 

DATA 



►ALEPH 97b,conf JTTSl (LEP 1) 



has been measured (data from 1994) for two samples: 0.6< <6.3 GeV^ 
(< g2 ^^2.8 GeV^) with < > equal to 0.15 and 3.0< <60.0 GeV^ (< >= 
14.6 GeV^) with < P"^ > = 0.22 GeV^. Results for these two samples are shown in 



fig. [T73- 
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Figure 177: The structure functions Fi]^'^^^^ /a as a function xsj for a) < Q'^ >=2.79 
GeV^, < P2 >=0.153 GeV^ and b) < >=14.65 GeV^, < >=0.225 GeV^. The soUd 
hnes are the QED expectations (from 
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Often variables referring to the 7*7 centre of mass frame are denoted by adding a star 
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In this experiment also the azimuthal angle distributions have been measured and 
ratios Fj('3^^)/i.7(Qi^^) p,mD) f p,iQED) ^^^^^^^^^ 



Results for 1994 Data 



Results for 1994 Data 





Figure 178: Results for a) pfQ^^D) ^ p^^iQED) 1j,i(QED)^j,^(qed) ^ ^ ^^^^^.^^ . 

at < >=8.8 GeV^. The solid lines are the QED expectations (from ^dW). 



•DELPHI 96a [|5|] (LEP 1) 

p^iQED) j^^^ been measured at < Q"^ >=12 GeV^, as a test for the unfolding and 
tagging methods in extraction of the hadronic F2 (see sec. 2.2.2| ). The effect of 
nonzero target virtuality has been studied. A satisfactory fit to the measured p^^^^^'> 
is obtained for the fixed value of P^=0.04 GeV^. 



•DELPHI 99,conf |]TT^ (LEP 1) 

The single-tag events with the muon pairs (data from 1992-1995) for between 
2.5 and 750 GeV^ are used to extract the leptonic photon structure functions. The 
SAT (1991-1993) and the STIC (after 1994) samples were studied corresponding to 
< >=12.5 GeV2, and the FEMC sample for < >=120 GeV^. 

The effect of nonzero target virtuality has been studied, a fit to F2 data gives for 
SAT (STIC) events = o.032 ± 0.007 (0.022 ± 0.007) GeV^. The results for various 
distributions for SAT and FEMC samples are presented in fig. |179. 
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Figure 179: Left: Distributions for STIC sample: a) Etag/Eheam, b) Q^, c) W, d) xsj, and a 
comparison with BDKRC and DIAG36 generators. Right: The same for the FEMC sample 
(from jjudi] ). 

The results for F2^'^^^^ /a as a function of xbj for the combined SAT and STIC 
data (< >=12.5 GeV^) and for FEMC sample (< >=120 GeV^), together with 
QED predictions are shown in fig. 
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Figure 180: The results for F2^^^^'^ /a as a function of xbj for the SAT and STIC combined 
data (< Q2 >=12.5 GeV^ ) and for FEMC sample (< >=120 GeV^ ), together with QED 
predictions (from 
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Azimuthal correlations have been studied to derive Fj^'^^^^ and p'^^^^\ in order 
to increase the observed correlation only events with 20° < 9* < 160° were taken 
into account. The azimuthal angle distributions (not shown) were fitted to obtain 
the factors multiplying the cos 0a and cos20a contributions. Next the extrapola- 
tion in these factors was made and the final results for ratios pj^^^^^ j p^^'^^^'> and 
p^mD)ip^{QED) obtained, see fig. m 
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Figure 181: The results for i7T;{Qi?^)/^7(Qi?i?) ^^^^j f;^(«^^)/2 ^2^(0^^) (bottom) as 
a function of xsj- The QED predictions are also shown (from Jii6| /j. 



•L3 98b jTTTI (LEP 1) 

The production of e^e~ and pairs for 1.4 GeV^ < < 7.6 GeV^ was studied, 

based on the 1991-1994 run. For a muonic final state, data on p^^^^^^ ^ and on the 
ratios pj(Q^^)/p^iQ^^) and ip^^(Q^^)/p^(Q^^)^ as well as on F^^^^^ and F^^^^^^ 
were obtained. 

The effect of the photon target virtuality has been studied and is clearly seen in 
all the three structure functions (see fig. |rS2]for F2^(^^^) only). 
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Figure 182: F^^'^^^\xBj,Q^,P^)/a forlA < < 7.6 GeV^ as a function of x Bj compared 
to QED predictions at P'^=0 (dashed line) and at < >fit= 0.033 GeV^ (solid line) (from 
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Figure 183: Results for Fg^^^^V" and F^^^^^V" measured at 1.4 < < 7.5 q^v 2. 



The curves are the QED predictions for P =0 (from illT]). 



•OPAL 93 [fTH (LEP 1) 

The QED structure function F^^^^^^ for < >=8.0 GeV^ was extracted from single- 
tag events at CM energy ~ Mz, for the first time at LEP. 

•OPAL 97f [m (LEP 1) 

Here the extraction of F'^^^^^^ was performed for 0.85< Q"^ <31 GeV^ (< > =5.2 
GeV2). 

•OPAL 99b []12g] (LEP 1) 

The complete OPAL dataset taken in years 1990-95 at e+e^ CM energies close to the 
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mass of the Z boson has been used to determine the QED structure functions p^'^^^^^ ^ 
pjiQED) p-^QED) quasi-real photon [i.e. single tagged events, with y < 0.5 and 
< 1.4 GeV^, with < >=0.05 GeV^) and, for the first time, the differential cross 
section for highly virtual photons with y < 0.5 and 1.5 < P^ < 20 GeV^. For 
single-tag events the different samples (SW, FD, EE) were collected, depending on 
the range (from 1.5 till 30 GeV^) and the calorimeter used to detect the scattered 
electron. The double-tag sample is denoted by DB and corresponds to between 1.5 
and 300 GeV^. 

The various distributions for the single and double-tag events were measured and 



compared with MC (Vermaseren generator) predictions. In fig. p.84| the results for DB 
sample are presented. Here the energy and the polar angle of the second electron, as 
well as the and P^ distributions, are shown. 
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Figure 184: Various distributions for double-tag events (from jJ2C|/j 



For single tagged events (with < P^ >= 0.05 GeV^) the differential cross section 
have been measured and the p^^^^^^ was extracted. Results for various samples 

are shown in fig. |185| {-^-) and [186| (P^'-*^^^''), together with the corresponding QED 
predictions. 
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Figure 185: Differential cross section da/dxBj for various < > samples of singly tagged 
events. The full line is the differential cross section as predicted by QED within Vermaseren 
Monte Carlo (from |l2^/). 
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Figure 186: Structure functions F2 /a for various < > samples of singly tagged 
events. The full line is the QED prediction for < >= 0.05 GeV^, the dashed line - for 
< P'^>=0 (from p^/j. 



For the first time the differential cross section has been measured for two highly 
virtual photons and the importance of the interference terms ttt and t^l (see eq. ^8[) 
has been shown. 

If, after integration of the differential cross section (eq. ^8]) over 0, the contributions 
of the terms ttt and ttl vanish, one can introduce an effective structure function of 
virtual photon (see sec. p.2.2| , PLUTO 84 and L3 2000). This is not the case in the 
present measurement, as can be seen in fig. |187| where Monte Carlo predictions also 
for the options with ttt=0 and ttt = ttl=0 are displayed and compared with the 
data. 

Based on azimuthal correlations, the QED structure functions pj^'^^^^ and p'^'^^^^ 
for quasi-real photons have been determined for < >=5.4 GeV^. First, ratios 
piQED) I p~f(QED) pf^QED) I p^iQED) ^^^^ obtained from a fit to the 0„ distribu- 
tion. Then from the measurements of F^^'^^^^ also the individual structure functions 
pjiQED) p'^QED) calculated. The obtained results are shown in fig. 
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Figure 187: Differential cross section for production witli two virtual photons, for two 

< >, < > samples, as a function of xbj- The full line is the differential cross section 
predicted by the Vermaseren Monte Carlo. The additional three histograms represent the 
cross section calculations with the GALUGA Monte Carlo for the full cross section (solid 
line), cross section obtained for ttt = (dot-dashed line) and that for ttt = ttl = 
(dashed line) (from p^j. 
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Figure 188: The measured structure functions F^^^^^^ /a and F^^''^^^' /a as a function of 
XBj for < > = 5.4 GeV^. The solid lines are the QED predictions for Q'^ = 5.4 GeV^ 
and e = p = 1 (from JJ20|/). 



Comment: For the single-tag events leptonic structure functions were extracted in the 
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wide range of . The effect of dependence was studied. 

For the first time the leptonic final state was studied for the deep inelastic scattering 
on the highly virtual photon. Due to large interference terms there is no clear relation 
between structure functions and cross sections. 

The definition of Fj^"^^^^ adopted here differs by sign and a factor 2 from that used 
by L3 98b. It is worth noting that the analytical formulae used here with the muon 
mass dependence, give noticeably different values from those with zero mass. 

"The measurements presented here supersede the earlier structure function results" 
^OPAL 93, 97{). 



For an overall comparison we present below collective figures (from ref. P|j) of 
the QED structure functions for the real photon together with corresponding QED 
predictions for P^ = 0. Data on p^^'^^^^ from old and new experiments are presented 
in fig. |189| as a function of xsj, for < > between 0.14 and 130 GeV^, and in fig. |190| 
as a function of for various < xbj > bins. Also the summary of structure functions 
ratios and and of the individual structure functions 

pi^ED) pf^ED) gj-^Q^j^ jj^ gg p^9i| and p.92| , respectively. 



191 



-M M M M 1 1 1 1 1 - 
: * TPC/2y(0.14-1.28) (a) I 


- 1 1 1 1 1 1 

: • OPAL(2.2) 


1 V 

(b): 


- 1 1 1 1 

: ■ L3(1.4-7.6) 


III- 

(cj 

>^ 




M M 1 M 




M 




: • oPAL(4.2) (d) : 


:*PLUTO(5.5) 


(e): 


: • OPAL(8.4) 


M 1 U 

(O: 














- 










-A 

7*, , , , 1 , , 









X 



- 1 1 1 1 1 1 1 1 1 - 

: iis CELLO(1.4-35) (g) = 


-III 1 

* DELPHI(4-30) 

- yi 

^' 1 ' 1 ' 1 




- 1 1 1 1 1 1 + 

• OPAL(12.4) (i) : 

1 ' 1 ' 1 '- 


; • OPAL(21.) 0) : 


: * PLUTO(40.) 


(k); 


: * DELPHI prl.(120.) : 


- • OPAL(130.) (m)P 

L A] 

- ^--^ ^ 





0.5 1 

X 

Figure 189: Summavy of existing Fi^^'^^^^ /a data as a function of xbj for broad < > 
range shown with QED predictions for P^=0 (from l^j). 



192 



z 

+ 

a 12 



o 
o 



10 



10 



10 



-— F2Vd(Q''P'=0) 


1 

^ T 


M| 1 1 1 1 
< X > 

■ 0.935 

— 0.85 


N_ 
10 

9 - 


_ y Ur^AL/ — ^ 




i-ir— 1 


. j 0.75 


8 


L| DELPHI ^ 






0.65 


7 - 


1 PLUTO ' 




^ 0.55 


6 


Jf CELLO ' " 








. )|^TPC 


^-1-5 — *f 


0.45 

0.35 


5 . 
4 . 




^ i -I ! \ 

- i i 




3 








I 0.25 
0.15 


2 - 








0.05 


1 - 






1 


..1 . . . . 





10 



10' 
[GeV^] 



Figure 190: Summary of existing F^^^^^^ /a data as a function of for fixed < xbj > 
bins shown with QED predictions for P^=0 (from ^j). 



193 



0.3 
0.2 
0.1 


-0.1 
-0.2 
-0.3 




+ OPAL 

+ L3 -- Q^ = 3.25GeV^ 

. j DELPHI pri.-.-.. = 12.5 GeV^ 



_J_ 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 



0.9 



b, 0.25 
S 0.2 
0.15 
0.1 
0.05 


-0.05 



i OPAL — = 5.4 GeV^ 

> L3 - - = 3.25 GeV^ 

I DELPHI prl. -- = 12.5 GeV^ 




_J_ 



(b) 




0.6 0.7 0.8 0.9 



Figure 191: Summary of existing data: a) Ff'^D)jp^(,QED) piQED) ^p^{QED) ^ 

function of xbj for < > = 3.25, 5.4, 12.5 GeV^, shown with QED predictions for P'^=0 
(from i/jj. 



B 0.2 
i?i 0.15 
0.1 
0.05 


-0.05 
-0.1 

-0.15 
-0.2 

-0.25 



L^^ OPAL 
14- L3 


— = 5.4 GeV^ 
-- Q^ = 3.25GeV^ 




(a) 1 














T 
I 








, 1 , 


, , , 1 , , , , 1 , , , 1 , , , , 5 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 



0.9 



5 0.2 
0.15 



0.1 
0.05 


-0.05 



: OPAL — 
L^^ L3 


= 5.4 GeV^ 


— ^ (b) ~. 


= 3.25 GeVi^-""""^ 








, , , 1 , , , , 1 , , , , 1 , , , ,- 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 



0.9 



Figure 192: Summary of existing data: a) FJ^^^^^ /a and b) F'^^^^^ /a as a function of 
XBj. Data for < >=3.25 and 5.4 GeV^ are shown with corresponding QED predictions 
for P^=0 (from 



194 



5 Selected results on the heavy quark content of 
the photon 

An additional information on the "structure" of the photon is coming from the produc- 
tion of heavy quarks in a photon-induced processes ( open heavy-flavor or heavy-Havor 
hadron production). In such processes heavy quark mass provides in a natural way 
the hard scale allowing to apply the perturbative QCD. However, the question how 
to describe within the perturbative QCD the heavy quark (hadron) production in 
hadronic processes is still open. Both the so called "massive" {no heavy parton) and 
"massless" (with heavy quark treated as a massless parton) schemes are being applied 
207 |-| 2T3| . For hard processes initiated by photons, where direct and resolved photon 



contributions are possible, this problem is particularly complex (see ref. ||214|| - 



In both approaches calculations for the heavy quark production in 77 and 7p colli- 
sions were performed up to the NLO accuracy. The fragmentation of heavy quarks into 
heavy-flavor hadrons is usually described by a simple Peterson function ||226|| (some- 



times with QCD evolution (DGLAP), see e.g. [P23|| ). The NLO calculation of the 
charm quark production in the DISe7 (contributions to F2 and F2) was done in | |165 



The progress in the theoretical description of heavy quark processes is necessary as 
large discrepancies in the description of the data involving heavy quarks are observed. 
This progress is needed not only for understanding of hadronic interaction. It is also 
crucial in searches of "new physics", e.g. for reliable estimation of the background. 

On the experimental side - new results from LEP and HERA colliders have ap- 
peared recently on the production of c-quark (and D*^ and Df mesons) as well as b 
-quark, in the DlSe-y and in the resolved photon processes in 77 and 7p collisions. 

As it was already mentioned in sec. ^ in early F2 experiments the heavy quark 
distributions, mainly c-quark, were subtracted from F2 data, see e.g. p[, AMY 95 



and TOPAZ 94 and references in |]T3|. The recent DIS-type measurement at LEP 
121|| led to the extraction, for the first time, of the charm contribution to F2 , denoted 
F2c- In fig- |193| the OPAL data in comparison with the NLO calculation ||165|| are 
shown. 
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Figure 193: OPAL data at < >=20 GeV^ in comparison with predictions of LO and 
NLO calculations and Monte Carlo models, a) The cross section for cc production as a 
function of xbj- b) The extracted structure function F2^{xBj, < >)/« (from / |J2J|/ j. 



The first hint that the charm-quark production in 77 colhsions requires the resolved 
photon processes can be found in AMY 1.5 94, ||122|| , ||123|| ; for other early results 
see ref. [p. 24 ] - \T^]. The new measurements of this type were performed recently at 
LEP by L3 group ||130|| , for the first time with 6-quarks production, see the collective 
figure, fig. |194| . The data arrive also from the c-quark (also D*^ and Df mesons) 

|132||. Some disagreement 



and 6-quark photoproduction measurements at HERA [|131| , 
with the Monte Carlo predictions and the NLO QCD calculations, especially for the 
6-quark production, is observed, see fig. |195| for ZEUS preliminary results. 
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Figure 194: Summary of charm and beauty production in 77 collision at e~^e~ collisions 
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Figure 195: ZEUS preliminary data. Left: Measured bb cross section for ep collision with 
p2(= q2) ^ I Qgy2^ o_2< y <0.8 and for Eti,2 > 7,6 GeV, |r/| <2.4, > 1.6 GeV,]!]"'] < 
1.1, compared to HERWIG and PYTHIA MC predictions. Right: Extrapolated b(b) cross 
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6 Related topics 



6.1 On the polarized photon structure functions 

The polarized photon structure function gj (another name is the spin dependent struc- 
ture function, i. e. not spin-averaged hke Fi and F2) both for a real and virtual photon 
can be studied in the DIS-type experiments with polarized electron beams and polar- 
ized photon target. The other spin dependent function, F3 , can be studied in the 
conventional experimental setup with unpolarized electron beams, see sec. 2.2. The 
structure function gj is of great importance, since its first moments (a sum rule for the 
"spin " of the photon) involve strong and electromagnetic anomalies, and it is deeply 



connected with the chiral properties of QCD |p3CI|| . 

The calculation of the structure function gj for polarized real and virtual photon 
has been done in NLO in ||227|| and ||228| , |229| | , respectively. The NLO calculations of 



the jet production in the polarized real 7-hadron collision can be found in |P31|] - ||234 



There is no data on the polarized photon structure functions , g{ nor on the 
polarized partonic content of the real and virtual photons from resolved photon pro- 
cesses. The future measurements in linear colliders, both in e^e~ and je or 77 options 



may be crucial for testing the structure of polarized photons, see e.g. ||235|, |133 



6.2 On the structure function of the electron 

In the DIS experiments in e^e~ collisions the inclusive hadron production can be 
ascribed not to the photon target, but rather to the parent electron or positron target. 
This in some cases may be more straightforward than pinning down the structure of 
the virtual photon, as was mentioned before. 

This topic is discussed in e.g. Ref. p|f, ||228|| . See also the results in ||236|| , where 



the structure function of the electron (in general - lepton) is related not only to the 
structure function of photon but also to the electroweak gauge bosons W and Z. The 
structure of weak bosons which appears in this approach was introduced and discussed 
in Ref. 1231 . 



6.3 On the photon content of the proton 

The photonic content of the proton relevant for the elastic and inelastic processes with 
an initial proton has been studied in | 238 | and [239, 24(]|| , respectively. 
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7 Summary and outlook 



The status of recent measurements (~ 1990 and later) of the "structure" of unpolarized 
photon in "DlSe^" experiments as well as in large pt jet production processes involving 
resolved photon(s) at e'^e~ and e^p colliders is presented. In updating our previous 
work p|i we have kept its form of the collection of selected data, presenting the main 
results of individual measurements. The final content of the survey has been dictated 
by the large amount of pioneering data that have appeared lately. 

For the real photon the data from LEP 1 and higher energy data from LEP 1.5 
and 2, with the higher statistics and improved unfolding methods are now published in 
final or preliminary form. Final results from TRISTAN appeared. The basic features 
predicted by the QCD are observed with higher than previous accuracy, as qualitative 
change has appeared. The hadronic structure functions for real photon have been 
measured in the 7*7 processes in a wide range of and down to very small xbj- 
On the other hand impressive progress has been obtained in pinning down the parton 
densities, in particular the gluonic content, in large pt resolved photon processes in 
ep collisions. The data are not only from the larger than before kinematical range 
but are also more precise, which allowed to pose more quantitative questions, e.g. on 
the rise of the F2 or gluon density at the small xbj- The prompt photon production 
has been measured in the resolved jp collisions at HERA (and also at TRISTAN 
in 77 collisions). This process can give an additional information on the partonic 
content of the photon as compared to the jet(s) production. The heavy (c and b) 
quark contribution to the photon structure functions has been measured in both e~^e~ 
and e^p collisions. 

The survey contains a large amount of data related to the "structure" of the virtual 
photon from the hard photon-proton collisions at HERA, where the effective parton 
distribution has been extracted from the dijet production. More than 15 years after 
the first measurement at PLUTO the corresponding new precise DlSg-y* measurement 
of the structure functions of virtual photon have been performed at LEP. 

The new measurements of the partonic content of the photon are accompanied by 
the impressive progress made in the NLO QCD calculations for the resolved real and 
for virtual photon processes. 

The discrepancies in describing the hadronic final states in the DlSe^ experiments 
as well as in the resolved photon processes in 77 and 7p collisions may have common 
origin. The need of additional pt in the standard QCD picture in the distribution of 
produced hadrons and jets may suggest an extra interaction involving the constituents 
of photon(s) - a multiple interaction. The newest data seem to give us a message 
that the partonic content of the photon is not properly described by existing parton 
parametrizations. 

The future high-energy linear e~^e~ colliders (LC) ||134|| as well as related 67 and 
77 colliders, based on the backward Compton scattering on the laser light, will offer 
a unique opportunity to measure the structure of photon in a new kinematical regime 



135|| . Moreover, at these colliders measurements of structure functions for photon 



with a definite polarization should become feasible with a good accuracy. To what 
extent the new generation of ep or 7J9 accelerators under discussion {e.g. THERA 



136 1) may extend our present knowledge, remains to be seen. 
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8 Appendix 



8.1 Parton parametrizations for the real photon 

Parton densities are given for fixed number of massless flavours (DO,LAC,GRSch,GIKO,IO) 
or for the number of flavours dependent on the scale (DG,GRV,GS,AGF,WHIT,SaS). 
Duke - Owens (DO) [|TM| 

A leading logarithmic parametrization of the parton distributions in an asymptotic 
form. Quarks with equal charges have the same distribution functions: /„ = /c, 

= (iV^ = 4). 
Drees - Grassie (DG) p5|] 

A parametrization of the full solution of the leading order evolution equations. The 
input parton distributions with free parameters are assumed at = 1 GeV^ and 
fitted to the only data on F2 existing at that time, at = 5.9 GeV^, from PLUTO. 
Field - Kapusta - Poggioli (FKP) [||] 

In this approach, there are two contributions to the nonperturbative ("hadronic") 
part, F2^^^, and the perturbative (pointlike) one, -F^^^. The F2^^ arises from the 
basic direct 7* qq coupling (and higher order QCD corrections), and it is included 
if the intrinsic quark pt is above a cutoff p^. If pt is smaller than p° then the qq 
pair creates a bound state leading to F2^^^ , usually taken from the VMD model. 
The perturbative -F^^'^ is calculated using the first order splitting functions and the 
one-loop as- 

Levy - Abramowicz - Charchula (LAC) 



A parametrization of the full solution of the leading order evolution equation fitted 
to all available in 1991 measurements of F2 for > Qq. The number of fiavours is 
fixed to Nf = 4. Three sets are provided with the different choices of an input scale 



Qq, and the a; ^ behaviour of a gluon distribution G{x), namely: 
• LACl: Q§ = 4 GeV2 

Qo = 4 GeV^, xG{x) const. 
Ql = 1 GeV2. 



LAC2 
LAC3 



Gordon - Storrow (GS) [IS? 



The LO and NLO parametrizations. The input structure function at the scale Qq = 5.3 
GeV^ [|1571| a and Ql = 3 GeV^ [|157|| b is chosen in the LO analysis as a sum of a hadronic 
part from the VMD model and of a pointlike part based on the Parton Model. Free 
parameters (also light quarks masses) are fitted to the data with Q^ > Q^. The input 
gluon distribution in LO [|157|] a is assumed in two different forms (set GSl and GS2). 
The NLO distributions in the MS scheme are obtained by matching of the F2 in the 
LO and the NLO approaches at the Qq scale. Nf = 3, 4, and 5 is used in [ p.57|] a. In 



TSTj b the number of fiavours equal to iV/ = 3 for Ql < Q'^ < 50 GeV^ and iV/ = 4 



for 50 GeV^ < was used 
Gliick - Reya - Vogt (GRV) [|15| 



The LO and NLO parametrizations of the parton distributions generated dynamically 
with the boundary conditions at Ql given by a VMD input. The physical photon is 
then assumed to be a coherent superposition of vector mesons, whose parton distri- 
butions are further assumed to be similar to those of a pion. The low initial scale 
Qq(LO) = 0.25, Ql{NLO) = 0.3 GeV^ is universal for structure functions of p, n, 7 



201 



etc. The DIS^ scheme is introduced to avoid the laige-x Bj instabihty problems. The 
one free parameter, which is the VMD input normahzation constant relative to tt, is 
fixed by the data. 

Aurenche - Chiapetta - Fontannaz - Guillet - Pilon (ACFGP) p9| 

A solution of the NLO evolution equation with the boundary condition taken at Ql 
= 0.25 GeV^. The input parton distributions were obtained from the VMD model at 

= 2 GeV^ and evaluated down to Q^. 
Aurenche - Guillet - Fontannaz (AGF) [ |16U|| 

The NLO parton distributions obtained with the input distributions (shown to be 
scheme-dependent) at Ql = 0.5 GeV^. The input distributions are based on the VMD 
model modified to agree with the MS scheme used in this analysis. 
Hagiwara - Izubuchi - Tanaka - Watanabe (WHIT) [[L61|| 
A set of six LO parametrizations obtained by fitting the input distributions to all 
available data for 4 GeV^ < < iqo GeV^. The parametrizations WHIT1-WHIT6 
are based on different input gluon distributions. A massive charm contribution is cal- 
culated from the Quark Parton Model below = 100 GeV^, above this scale from 
the massive quark evolution equations. 
Schuler - Sjostrand (SaS) p2l 

Four sets of the LO parametrizations. The nonperturbative input distributions at 
Qo = 0.6 GeV in the SaSlD, SaSlM and Qo = 2 GeV in the SaS2D, SaS2M sets are 
based on the VMD model (their normalization is fixed by the momentum sum rule, 
and the x - dependence is obtained from the fits to the data). The fully calculable 
pointlike contribution to F2 is expressed as an integral over the 7* —>■ qq virtuality fc^. 
The nonleading term is included into the F2 , leading to the MS distributions 
(SaSlM, SaS2M). = gives distributions in the DIS scheme (SaSlD, SaS2D). 
Gliick - Reya - Schienbein (GRSch) 

The LO and NLO (DIS^ scheme) parameter-free parton distributions obtained from 
the VMD input based on a coherent superposition of vector mesons applying the new 
Gliick - Reya - Schienbein pionic distributions |p.64|| . The low initial scale (universal 
for the proton, pion and photon distributions) Qo(LO) = 0.26, Qg(NLO) = 0.40 GeV^ 
is assumed. The parton distributions are calculated for light quarks, Nj = 3 (contri- 
bution of heavy quarks is taken into account in the structure function ^2^) . 
Laenen - Riemersma - Smith - Neerven (LRSN) ||165|| 

The NLO calculation of the structure function F2 and F2, including the heavy quark 
contribution. 

Gorski - loffe - Khodjamirian - Oganesian (GIKO) ||166 



The structure function F2 of the real photon, obtained from the results for the virtual 
photon, where beside PM also a gluon condensate contribution is included (see below). 
Only the ligh quarks contributions are considered {i.e. Nf=3). This approach is valid 
for Q'^ below 10 GeV^ due to the lack of evolution, and for intermediate xbj values: 
0.2 - 0.7. 

loffe - Oganesian (lO) ||167| 

Gluon distributions in the real photon calculated using the LO and NLO evolution 
equations, with the light quark densities input from GIKO. Two sets of gluon density 
are studied: 1) without the intrinsic gluons at the input scale Ql=l GeV^, 2) with the 
gluon density at the input scale equal to gluonic content of vector mesons p,uj,(f). The 
range of applicability: XBj below 0.6 down to 0.03-0.05, and 1 < < 100 GeV^. 
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Here also results for F2 of the real photon, with the LO and NLO evolutions in 
are given (with A^/=3), vahd for 0.2 < xsj < 0.7 and 1 < < 50 GeV^. 
H. Abramowicz, E. Gurvich, A. Levy (GAL) [|169|| 

The LO ||169|| a and NLO ||169|| b parton densities in the photon are extracted from the 
existing F2 measurements and the \ow-x proton structure function (related to the F2 
by Gribov factorization at low-x). 



8.2 Parton parametrizations for the virtual photon 

The notion of the partonic content of the virtual photon (P^ 7^ 0) can be applied 
if scale is larger than (typically ^ 0.2 Q^). The parton distributions in 
the virtual photon can be calculated in the perturbative QCD for ^ ^ 
without any input, for lower P^ extra assumptions are needed. The parametrizations 
described below are valid for < P^ and become the corresponding parametrizations 
for the real photon in the limit P^ — > 0. In the low-P^ region some of them are based 
on interpolation between parton distributions in the real photon (P^ = 0) and the 
QCD predictions in high-P^ region (A^ ^ P^). All the parametrizations presented 
below deal with the transverse polarized virtual photon, with one exception (Chyla 
parametrization) . 
Drees - Godbole |]f 

The parametrization of the parton distributions in the electron and in the virtual pho- 
ton. The parton distributions in the virtual photon are equal to the corresponding 
distributions in the real photon {e.g. DG, GRV) multiplied by a factor dependent 
on virtuality P^. In the limit P^ ^ ^ h? the quark distributions in the virtual 
photon approach the leading logarithmic (In^) parton model predictions arising from 
7* ((5^)7* (P^) — > qq subprocess. The gluon distributions in this limit approach the 
- behaviour (~ In^ ^2) following from the 7*(Q^)7*(P^) — > gqq cross section. The 
parametrization contains one free parameter P^?, a typical hadronic scale (in applica- 
tions denoted also as u"^). 
Gliick - Reya - Stratmann (GRS) [ |T70| 

The LO and NLO (DIS^ scheme) distributions obtained by solving the Q^-evolution 
equation with boundary conditions being a smooth interpolation between boundary 
conditions valid at high and low P^. For high P^ —>■ ^ the solutions correspond 
to the cross section for 7* ((5^)7* (P^) qq subprocess. For P^ = they are given by 
VMD input parton distributions at the scale Ql{LO) = 0.25 or Qg(NLO) = 0.3 GeV^ 
as for the real photon in | |158[ . The applicability of the virtual photon parametriza- 



tions is assumed for < P^ < 10 GeV^ and 5P^ ^ Q^. The number of active flavours 
is fixed to Nf = 3. 
Schuler - Sjostrand (SaS) [|T62] 

An extension of the SaS parton distributions in the real photon to the virtual photon 
case, valid for < P^ < Q^. In the pointlike contribution, the integral over the virtu- 
ality k"^ of the 7* — s> qq state, is modified by a factor (^)^ with the lower limit of the 
integration equal to a cutoff scale Pq (several choices for Pq have been considered, e.g. 

2 

Pq = max{P'^,Ql}). In the hadronic contribution a factor ( ^r^p2 Y is introduced, 
where my is a vector-meson mass. The notation SaSlD, SaS2D, SaSlM, and SaS2M 
is the same as in the real photon case. 
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Gliick - Reya - Schienbein (GRSch) 

The LO and NLO parton distributions valid for < ^ and Nf = 3. It is argued 
that whenever the virtual photon is probed at ^ it may be treated as the real 
one: the longitudinal polarization should be neglected and a partonic cross section 
involving 7(P^) should be calculated as if = (in this approach the virtuality 
P^ 7^ is present in the flux e —>■ e^* and in the virtual photon parton distribu- 
tions). The boundary conditions in the limit ^ ^ h? should be calculated 
from 7*(Q^)7(P^ = 0) — gg subprocess instead of 7*(<5^)7*(P^) — ^ qq] this implies 
changing the boundary conditions compared to the previous approach ||17CI|| . Now, at 



P^ = ^ the parton distributions vanish in both LO and NLO in DIS-y scheme. 
The input distributions at the initial scale = max{P'^, Q^} are assumed to be equal 



to the VMD parton distributions in the real photon (see sec. |8.1| ) multiplied by the 
factor [-^^^f, where nip = 0.59 GeV^, g^(LO) = 0.26 GeV^ and Qg(NLO) = 0.4 
GeV2. 

Gorski - loffe - Khodjamirian - Oganesian (GIKO) ||166 



The structure function P^ of the virtual photon, for P^ > 0.5 GeV^, calculated using 
the Operator Product Expansion as a series in 1/P^, where beside PM also a gluon 
condensate contribution is included. Using a dispersion relation the structure func- 
tion is then presented in terms of the contributions from physical states (p, u, 0) and 
continuum. The proper continuation to P^=0 can be performed, leading to prediction 
for the real photon (see sec. |HT^). Only the light quark contributions are considered. 
This approach is valid for intermediate xbj values: 0.2 - 0.7, and for larger than 
P^ but below 10 GeV^, due to the lack of evolution, 
loffe - Oganesian (lO) |T67| 

Gluon distributions in the (transverse) virtual photon calculated using the LO and 
NLO evolution equations, with the light quark densities input from GIKO. Two types 
of gluon density are studied: 1) without the intrinsic gluons at the input scale, Ql=l 
GeV^, 2) with the gluon density at the input scale, Ql=l GeV^, equal to the gluonic 
content of the vector mesons p,uj,(f). The range of applicability: xbj below 0.6 down 
to 0.03-0.05 and 1 < < 100 GeV^. 

Here also results for F2 of the (unpolarized and transverse) virtual photon, for 50 
GeV^ > P^ > 0.5 GeV^ with the LO and NLO evolutions in Q'^ are given (with 
Nf=3). Results for F2 for a real photon are also given (see above). 
Chyla [|l68f| 



The LO parametrization of the pointlike quark and gluon densities in the longitudinal 
virtual photon, 1 GeV^ < P^ < 0.2 Q^, obtained from the homogeneous evolution 
equation, without extra input. For P^ — »• the gluon distribution Gl in the longi- 
tudinal virtual photon vanishes, while the quark densities gf = q^^ approach the PM 
formula for 7* ((5^)7* (P^) — > qq- Parton densities satisfy the momentum sum rule. 
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